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abstract 


The adsorption and desorption characteristics of activated carbo 
have been experimentally and theoretically investigated with 
following overall goal, To determine the applicability . regene - 
able sorption to the control of airborne trace contaminants wi 
spacecraft cabins for long mission durations. Capacity correiations 
pure and mixed contaminants have been established, h theory has 
been derived and successfully tested for vacut. desorption rates 
from single particles and beds. Finally, a quantitative design 
methodology has been developed for the practical design of regener- 


able systems. 
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trace contaminant adsorption and sorbent regeneration 

By A. J. Robe 11, C. R. Arnold, A. Wheeler, G. J. Kersels, 

and R. P. Merrill* 

Lockheed Palo Alto Research Laboratory 


Section 1 
INTRODUCTION 


Xa any closed ecological system the buildup of trace contaminants from both metab 
processes and material degradation must be controlled by a suitable contaminant 
Umoval system. If a sorbent is used in such a system for extended manned space 
Missions, unduly large quantities of sorbent will be required unless a regenera iv 

Operation is employed. 

The overall goal of this work was to acquire sufficient information on adsorption/ 
desorption phenomena to enable design, optimization, and operation of a regenerable 

system. 

Sorption processes on activated carbon were studied with the goal of arriving at 
quantitative relationships enabling the correlation and estimation of sorbent capacRy 
interference effects of multiple contaminants (including water vapor), and desorption 
rates under spacecraft regeneration conditions. 

Despite the use of activated carbon as a sorbent for many years, little or no quantita- 
te work had been done hi the above areas at the time of inception o, this program. 

In fact, the state-of-the-art was as follows: 

• Activated carbon was the prime candidate for use as a general sorbent, and 
Bamebey-Cheney type BD was found, during performance of Contract NAS 
9-3415 (Design and Fabrication of a Trace Contaminant Removal System 
for Apollo) to be quite suitable for spacecraft use. 

^Assistant~Professor of Chemical Engineering, The University of California, Berkeley 
(CoMvdtantat Lockheed Palo Alto Research Laboratory). 
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• The potential theory had been used to correlate adsorption capacities of 
single contaminants. It had not been extended to mixtures. There was no 
other general theory applicable to the mixture problem, nor were there 
mixture data relative to conditions of interest for spacecraft operation. 

• Vacuum desorption rates were unknown. _ 

I 

• Concepts for system design were not based on a fundamental understanding 1 
of sorption phenomena; design was performed more or less by rule-of- thumb | 

A systematic approach was used to fill these gaps in order to arrive at sound design 
principles. It was felt that an understanding of sorption by a single particle was 
fundamental to an understanding of how a collection of particles functions in a bed. 

Thus, a microbalance system was used to study single particles in a parallel effort 
to a flow system for beds. 

Theoretical extention of the potential theory to very low coverages was undertaken to 
attain a better understanding of the empirical parameters of the correlation and in its 
extension to estimations of some of the necessary kinetic parameters. This approach 
was merged with the satisfactory theoretical analysis of mixed adsorption to give a 
generalized capacity correlation. The mathematics of transport of sorbate from a 
bed of particles during vacuum desorption was developed. All of these theoretical 
advances were combined with the empirical data produced during the contract to 
indicate a rational design procedure. 
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Section 2 

ADSORPTION OF PURE CONTAMINANTS 


The length of time that a trace contaminant system will be usable without regeneration 
is directly related to the capacity of the adsorbent material for the various contam- 
inants to be encountered. In order to gain a better understanding of the relation of 
this capacity to contaminant levels, it was necessary to determine extensive adsorp- 
tion isotherms, which show the relation of adsorptivity to contaminant concentration 
or partial pressure. These were determined for a few selected components which 
were felt to represent a typical spectrum of contaminants that are amenable to re 
moval by adsorption on activated carbon. 

Since the method of choice for trace contaminant removal is a fixed adsorbent bed 
in a flowing stream of space cabin air (actually a 50/50 mixture of nitrogen and 
oxygen at one-half of atmospheric pressure), a number of experiments were conducted 
using a flow system incorporating a fixed bed scaled down to laboratory dimensions. 
Using this method the performance of the sorbent bed could be determined under 
conditions similar to those encountered in a spacecabin atmosphere, but using only 
one contaminant at a time at a concentration near its maximum allowable concentration 
(MAC) value. It was also necessary to use this experimental approach to obtain sorbent 
beds saturated with typical contaminants at these concentrations in order that these 
beds might then be tested for their vacuum desorption characteristics, as described 

in Section 7 . 

Due to the time-consuming nature of the flow experiments, a faster method of determin- 
ing more extensive adsorption isotherms of the pure components was needed. This 
was accomplished by use of the gravimetric system described below, in which an 
isotherm that might take weeks to determine by flow methods could be determined in 
1 or 2 days. The core of this system was the microbalance unit, which provided 
instantaneous weight readings as adsorption partial pressure was varied at will. The 
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weights were then translated directly into adsorptivity or sorbent capacity values. 

By having these extensive isotherms for pure contaminants, it was then possible 
to test various theoretical models (Section 4) that would enable prediction of sorbent 
bed performance under real conditions in which one or more contaminants would be 
present at one time. These isotherms were also of value in correlating desorption 
(regeneration) data in Sections 6 and 7. 

2.1 EXPERIMENTAL TECHNIQUES 

2.1.1 Flow System 

This system was designed so that a stream of gas contaminated with a known concen- 
tration of adsorbate or adsorbates could be passed through a cylindrical packed bed 
of sorbent at subatmospheric pressure while effluent from the bed was monitored 
for the contaminant. A schematic diagram of the system is shov-'x in Fig. 2-1. 

The sorbent bed consisted of a 60-cm stainless steel tube, of 1/4-, 3/8-, or 3/4-in. OD 
with stainless steel screen spacers used to hold both ends of the sorbent bed in place. 

(A slightly different bed design was used in early runs in the F2- and F3-series.) 

That portion of the tube upstream from (below) the sorbent bed was filled with quartz 
chips 1 to 2 mm in diameter in order to preheat the gases entering the sorbent bed in 
experiments at elevated temperatures. In flow run F13 and thereafter, the quartz 
chips were eliminated because it was of interest to obtain the upstream and down- 
stream pressure of the bed during vacuum desorption. The quartz chip packing would 
result in misleading upstream pressure readings. A drawing of the flow bed is shown 
in Fig. 2-2. 


Those portions of the flow system upstream from the sorbent bed were constructed 
from l/4-in. OD copper and stainless steel tubing, in order to minimize dead volume, 
and were connected with Swagelok fittings. Portions downstream were made of 
1-in. OD copper tubing with sweat soldered fittings in order to provide high 
conductance for vacuum desorption of the bed. 
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Fig. 2-2 Adsorbent Chamber 
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Gases were introduced into the system through manifolds, as illustrated m Fig. 2-1. 
Manifolds 1 and 2 allowed a premixed carrier gas without contaminant to be instanta- 
neously substituted for a similar gas with contaminant, without any significant change 
in flow rate. Manifold 3 provided for addition of carrier gas components in order to 
adjust their concentrations or those of the contaminants. The gas cylinders were 
equipped with standard pressure regulators and, except for those cylinders containing 
premixed contaminated gases, with molecular sieve 13X traps for removal of traces 
of extraneous contaminants. 

Effluent gas samples were taken automatically, on a preset time cycle, with a gas 
sampling valve. Such valves have the advantage of giving very reproducible injection 
volumes, and they have low maintenance needs. 

A vacuum pump placed downstream from the gas sampling valve allowed the system 
to be operated at subatmospheric pressure. This pressure was maintained at approxi- 
mately 0. 5 atm by judicious adjustment of the needle valves located on each inlet 
manifold and at the outlet to this system vacuum pump. 

Valves used in the system were of three types: toggle valves for on-off control of 
upstream gas flows; needle valves for precise control of flow rates; and screw-type 
bellows valves for bypass and vacuum portions of the system. Specific locations of 
these valves are shown in Fig. 2-1. 

Since it was found early in this work that rotameters provide marginal measurement 
reproducibility, total system flow was monitored with a wet test meter located on the 
exhaust of the system pump. Average flow rates were determined by dividing the 
difference in wet test meter readings by elapsed time. 

The sorbent bed was heated by a Hevi-Duty tube furnace in which it was vertically 
mounted. Furnace temperature was kept constant by two West JP controllers, while 
the skin temperature of the bed was measured with a chromel-alumel thermocouple 
attached near the midpoint of the bed. This thermocouple was connected to a tem T 
perature recorder, which provided graphic temperature records. (See Fig. 2-2.) 
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t„h ;rrr ure was measured with waua ° e and -odd ** ™ ^ ri 0 n- 

tabe differential pressure gages. A model FA129 gage, or in later experiments a 
Barocel capaeitanee manometer, was used to measure the pressure immediately 

upstream from the sample injection valve to determine the pressure of the standard 
vo ume o gas injected into the chromatograph. 

7o“ o“fr T‘ ySi t 01 ^ e " 1Uent gaS6S WaS prOTided bF “ F « model 
810 gas chromatograph equipped with a flame ionisation detector. Current 

output of this detector was converted by an electrometer to a voltage that was dis- 
played on a strip chart recorder. This recorder was equipped with a disc chart 
integrator for evaluation of the total ionisation charge generated by contaminant 
ructions coming out of the chromatographic columns. Separation of effluent gas 
components for analysis by the detector was effected by the use of appropriate 
columns, flow rates, and temperatures. 

Raw sorbent materials obtained from outside sources were crushed and sieved to 
sises differing by a faetor of two. Specifications for these sorbents are given in 

Table 2-1. Ix>cksorb-2A is a Lockheed proprietary sorbent prepared by impregnating 
30 x 40 BD with a Group 1 metal. 


Table 2-1 

SORBENTS STUDIED 



-j 


Manufacturer 

Type 

Sieve 

Size 

(U.S. No.) 

Approx. 

Diam. 

(cm) 

Surface 

Area 

(m2/ g ) 

Textual 

Designation 

Barnebey- Cheney 

Union Carbide, 
Linde Division 

Lockheed 

BD, activated 
coconut shell 
charcoal 

Molecular 
Sieve 13X 
(Zeolite) 

Locksorb-2A 

8 x 12 
16 x 20 
30 x 40 
50 x 70 

16 x 20 
30 x 40 

0.20 
0.10 
0.05 
! 0.025 

0.10 

0.05 

1100 

1100 

1100 

1100 

725 ± 75 

8 x 12 BD 
16 x 20 BD 
30 x 40 BD 
50 x 70 BD 

16 x 20 
LMS 13X 

LS-2A 


A mixture of 50% oxygen and 50% nitrogen was used as the carrier gas in flow experi- 
ments. Contaminants to be studied were obtained premixed in the above carrier or in 
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oxygen. These premixed gases were supplied by the Matheson Company with their 
analyses of the contents. 


For initial sorbent outgassing, samples were weighed out in air and placed in the 
60— cm tube described above. The tube was then inserted into the system and sub- 
jected to a flow of helium. Bed temperature was then raised to 400-450° C, while 
system pressure was maintained at approximately 0. 5 atm. After 4 hr the temperature 
was allowed to come to ambient or to the temperature at which the experiment was to be 
run. In a few cases it was found that after the heat pretreament small air leaks de- 
veloped around the O-ring fittings which hold the bed in place. (See Fig. 2-2. ) To 
assure that no leaks were present after the pretreatment, the bed was evacuated to a few 
microns and isolated. The increase in pressure in the bed tube was measured as a 
function of time and the leak rate , if any , was determined. 

In some cases this pretreatment was performed at a system pressure of 1 atm, but this 
did not appear to affect the results as demonstrated in subsequent flow runs. 


Adsorption experiments on the flow system consisted of five steps: 

(1) Calibration of the chromatograph detector for each contaminant to be 
studied in that experiment 

(2) Determination of the inlet concentration of each contaminant, by allowing 
the flow to bypass the sorbent chamber 

(3) The actual adsorption, during which effluent from the bed was monitored 

for each contaminant 

(4) Redetermination of inlet concentration by again bypassing the sorbent bed 

(5) Re calibration of the detector for each contaminant 

Vacuum desorptions, if any, were next performed. Finally, flow desorption was 
performed by following steps (1), (3), and (5) of the above procedure. 


Calibration was performed by measuring detector response to a known volume of 
premixed gas of known concentration. After calibration for each contaminant, all 
premixed gases were then allowed to flow, along with nitrogen or oxygen-nitrogen 
mix needed to adjust contaminant concentrations and keep the oxygen-nitrogen ratio 
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at or near unity. The adsorption step was performed by merely diverting the flow 
through the sorbent chamber and observing chromatograph detector response to the 

contaminant in question until the effluent concentration appeared to have risen back 
to the inlet concentration. 

Flow desorption was performed, after calibration of the detector, by passing the oxygen- 
mtrogen mixture through the sorbent bed and monitoring the effluent as in the adsorp- 
tion step. Strongly adsorbed contaminants were sometimes removed in a similar 
manner at higher temperatures, using helium instead of oxygen- nitrogen. 

Bed effluent contaminant concentrations were determined by measuring the area of 
peaks on strip chart recordings corresponding to detector output for each contaminant. 
Since holdup time for each contaminant under study had been determined before begin- 
ning the experiment, location of the peak on the recording time axis was relatively 
simple. Peak area was determined simultaneously by means of an integrator built 
into the recorder. Records were kept of peak area, electrometer range and attenua- 
tion settings, injection pressure, wet test meter readings, and time. A sample 
calculation of the amount adsorbed appears in Appendix A. 

Since detector response to the contaminant varied slightly during the course of each 
experiment, it was evaluated at two times, as indicated by steps (1) and (5) on page 9 
Average response values before and after the run were then interpolated in a linear 
fashion and used to determine effluent concentrations during the course of the experi- 
ment. Sample calculations of electrometer response and of effluent concentration 
also appear in Appendix A. 

Linearity of detector response to contaminants in the concentration range of interest 
was determined by injecting premixed propylene (50 ppm) in oxygen-nitrogen into the 
chromatograph and observing detector output. Variation of the amount injected from 
one-tenth to ten times the volume of interest resulted in a log-log plot with slope 
equal to unity, as shown in Fig. 2-3 , thus indicating completely linear response over 
this range. Injection of larger volumes of gas tended to blow the detector flame out, 
while injection of smaller volumes led to poor reproducibility. 
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Fig. 2-3 Calibration Curve for Propylene in F&M 700 Flame 
Ionization Detector 
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2.1.2 Gravimetric System 


The sorbent materials used were the same as those described in section 2. 1. 1. 
Gravimetric experiments were performed with a Cahn model RG microbalance. 

This microbalance has an ultimate sensitivity of 0.2 and can measure weight 
changes as large as 200 mg on its most sensitive beam position. Its sensitivity is 
1 pg per chart division while measuring weight changes as large as 10 mg. A vacuum 
chamber accommodating the microbalance was constructed, as shown in Fig. 2-4. 

This gravimetric system consisted of a number of functional parts, most of which 
are typical of any vacuum system. The jar in which the weighing mechanism was 
mounted was equipped with provision for three hangdown tubes, one for each beam 
position. In the present work, only the more sensitive position was used, leaving 
the hangdown tube port for the less sensitive position available for attachment of a 
Barocel capacitance- manometer pressure-transducer system, which was used to 
obtain the pressure data. In earlier experiments this port was used for attachment 
of a McLeod gage, which was used to calibrate a Pirani gage used for pressure data. 
This calibration appears in Fig. 2-5. 

An ionization gage was provided for checking ultimate vacuum during desorption 

steps, and for determining apparent system leak rates, which were generally of the 
~ '5 

order of 10 to 10 torr per minute for a system volume estimated at 5 liters. 
Pumping was provided by a 4- in. water-cooled oil diffusion pump with a baffled liquid 
nitrogen trap which was connected to the system by a 4-in. Viton-sealed gate valve. 
Adjustment of range and sensitivity was provided by the microbalance control unit, 
which was connected to the balance mechanism by electronic cables and appropriate 
vacuum-type connectors. Readout was provided by a strip chart recorder attached 
to the output of the control unit. 

Introduction of vapors into the leak valve was effected by connection of a manifold, 
which is shown schematically in the system line drawing in Fig, 2-6, This manifold 
had provision for introduction of both liquid and gaseous adsorbates. Connections 
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Fig. 2-4 Gravimetric Adsorption System 











Fig. 2 5 Calibration of Pirani Gage in n-Butane 
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Fig. 2-6 Gravimetric System Line Drawing 





to a vacuum pump were provided so that the manifold could be quickly flushed with 

different vapors. Pirani-gage tubes were provided in each section of the manifold 
so that leak rates could be estimated. 


Heating was provided by two identical cylindrical tube furnaces, which were located 
by means of adjustable clamps around each hangdown tube in the vicinity of the sample 
and counterweight pans, as shown in Fig. 2-7. Temperature control was provided 
by a West JP temperature controller on the sample furnace and by Powerstat auto- 
transformers in series with windings of both furnaces. The temperature controller 
sensor was an iron-constantan thermocouple placed in the sample tube furnace. 
Temperature readout was provided by another iron-constantan thermocouple located 
directly opposite the sample pan at the hangdown-tube surface, and connected to an 
accurate dc vacuum-tube voltmeter (VTVM). 


The microbalance was calibrated periodically to read from 45 to 55 mg with ± 1 M g 
precision. Charcoal sample weights were kept as close to 50 mg as possible. Since 
adsorbed atmospheric contaminants were generally about 5 percent by weight, initial 
sample weights of about 52. 5 mg were weighed out in air with the hangdown tube off. 
After replacing the hangdown tube and putting the furnaces in place around each tube, 
the system was rough-pumped for 5 to 10 min and then opened to the oil diffusion pump 
through the gate valve. After the sample weight rate of decrease was sufficiently 
low (10-20 pg/hr), both hangdown tube furnaces were turned on, first to low heat, 
which was gradually increased to about 400° C, where it was maintained for a few 
hours to get constant weight. The furnaces were turned off and allowed to cool for 
at least an hour before adjusting the temperature to the desired value. 

To determine isotherms, first, the introduction manifold (Fig. 2-6) was flushed two 
or three times with the vapor to be used by first admitting the vapor and then pump- 
ing it out with the auxiliary pump. Next, the system gate valve (Fig. 2-4) was 
closed and the ionization gage (P 3 ) filament turned off, after noting that the system 
leak rate was reasonable (less than 10~ 5 torr/min). Then the vapor was admitted 
to the system by again admitting vapor to the manifold and opening the leak valve 
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sufficiently to allow system pressure to be increased to the lowest desired value 
for an isotherm point. After equilibrium had been obtained, as evidenced by constant 
weight and pressure (see strip-chart records in Appendix B), the weight and pressure 
values were noted and pressure was increased to the next desired value. In some 

cases, pressures were approached in the opposite direction by cracking the gate valve 
slightly in order to pump out only a part of the vapor. 

Temperature was maintained by judicious adjustment of the temperature controller. 

In experiments at 27°C, the furnace wrapping was kept damp with water to provide 
a "wet-bulb" cooling effect, which was then counteracted by heating at very low voltage 
settings, which were judiciously changed to maintain temperature within ±0. 2°C. 

After sufficient points had been obtained to determine an isotherm, the system was 
pumped out completely by opening the gate valve fully. In cases where more than 
one point was obtained at a given pressure, a slight weight increase was noted at lower 
temperatures and with lighter contaminants, such as n-butane and propane. This was 
assumed to be due to a background of "heavy" contaminants, and was extrapolated 
over the period of the isotherm experiment as a corrective addition to the initial 

weight used in calculating the amount adsorbed. In most cases this correction was 
very minor. 

2.2 RESULTS 


Isotherm data were obtained for a number of contaminants on the various sizes of 
BD charcoal listed in Table 2-1. These data were obtained on both the flow system and 
the gravimetric system. Summaries of data obtained on these systems appear in 
Appendixes C and D, respectively. 

Selected gravimetric data for n-butane and toluene were plotted according to the 
Freundlich relation, 


q 


a 



(2-1) 
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where 


q & = amount adsorbed (ml STP/g) 

k_, = a constant [g/ml STP (torr) ] 

F 

p = partial pressure of contaminant (torr) 

n = a constant 


An advantage to plotting the data according to the Freundlich isotherm is that increased 
accuracy is obtained in reading low pressure and amount of adsorbed values. Such 
accuracy is required for calculation of thermodynamic state functions. 

Early data from gravimetric experiments with n-butane on 0. 5-mm charcoal particles 
based oh Pirani-gage readings at various temperatures are plotted in Fig. 2-8. 

Later data based on Barocel readings for 2. 0-mm particles at 27 and 67 °C and for 
0.5- and 1. 0-mm particles at 27 and 107° C are shown in Fig. 2-9. 

Data for toluene on 1. 0-mm charcoal particles at 27, 67, and 107 °C and on 0. 5- and 
,1.0-mm particles at 27 °C are shown in Fig. 2-10. Some additional isotherm data 
were determined for various contaminants in the course of desorption experiments de- 
scribed in section 5. 1. 1 and are tabulated in Appendix E. 

Selected isotherm data, which were obtained in the flow adsorption experiments 
described in section 2.1.1 and summarized in Appendix C, are also tabulated separately 
in Appendix F. In addition, some adsorption data for n-butane at much higher pressures 
which were obtained in other work are tabulated in Appendix G. 

The data in Appendixes E and F that were obtained under conditions similar to those 
of the isotherms in Figs. 2-9 (n-butane) and 2-10 (toluene) are plotted along with these 
isotherms. Open symbols are used to indicate gravimetric data, while flow data are 
shown as solid symbols. 

Isotherm data were obtained for water on 30 x 40 BD at pressures near the vapor 
pressure of water at 25° C: 5, 10, 15, and 20 torr. These data are plotted on linear 
coordinates in Fig. 2-11. 
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Fig. 2-8 Isotherms of n-Butane on 30 x 40 BD Activated Carbon 
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AMOUNT ADSORBED (ml STP/g) 




AMOUNT ADSORBED (ml STP/g) 








Section 3 

ADSORPTION OF MIXED CONTAMINANTS 


While adsorption isotherms such as those determined in Section 2 give one an excellent 
feeling for the performance of sorbent beds under conditions in which only one 
contaminant is present in the air , they do not necessarily present sufficient informa- 
tion as to how the presence of other contaminants will affect adsorptivity. The next 
logical extension of the experimental work then was to perform experiments in which 
there were two contaminants present. To this end, and also with thought to the vacuum 
desorption experiments on fixed beds described in Section 7, a number of flow experi- 
ments, similar to those described in Section 2 but with two and in one case three 
contaminants, were performed. Pairs of contaminants with both similar and widely 
different adsorptivities were selected for these experiments. Since water is not 
normally considered a contaminant, experiments incorporating it as the second 
component are described separately in a later section (Section 9). 

Due to the time-consuming nature of the flow experiments, which were exaggerated 
greatly by using two contaminants of widely varying adsorptivities, a faster method 
was needed to get extensive data on adsorptivity of two contaminants at one time. The 
microbalance apparatus described in Section 2, this time modified by inclusion of a 
mass spectrometer for determination of partial pressures of multiple component 
vapors, was again found to provide adsorption isotherm data in much shorter 
experiments. Two typical contaminants of quite different adsorptivities in the pure 
state were chosen for these experiments. These experiments provided valuable data 
for testing a general theory for prediction of sorbent performance in the presence of 
more than one adsorbing component, as described in Section 4. 
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3.1 EXPERIMENTAL TECHNIQUE 


3.1.1 Flow System 

The same system and procedures were used in multicontaminant flow adsorption as 
in single contaminant work. There is only one variation in the experimental procedure. 

In blockage and displacement experiments, the inlet concentration of one contaminant 
was kept at zero until the other contaminant had been adsorbed on the bed. Oxygen 
or oxygen-nitrogen was substituted for that contaminant from just prior to the beginning 
of the adsorption step until some convenient time after adsorption of the first contaminant 
was complete. In flow experiment F17-3 (see Table 3-1), simultaneous adsorption of 
n-butane and toluene was studied; therefore, the oxygen-nitrogen mode described above 
was not used. 

3.1.2 Gravimetric System 

The apparatus used to determine mixed- contaminant adsorption isotherms was the 
same as that described in section 2. 1. 2 for single- contaminant isotherms with the 
addition of a mass spectrometer section, as shown in Fig. 2-4. This mass spec- 
trometer, a Varian Quadrupole Residual Gas Analyzer (QRGA), was included in a 
sampling manifold with a 15-liter/sec getter- ion pump at the outlet and with molecular- 
flow conductance limitations at both the inlet and outlet so that the relative concentra- 
tions of each species were the same in the mass spectrometer as in the microbalance 
jar. (See Appendix H.) Mass spectrometer readout was obtained by connecting a 
strip-chart recorder to the picoammeter (electrometer) used to amplify ion currents 
that were detected as various mass-to-charge ratios were electronically swept over 
the range of interest. The strip-chart recorder, dual-pen version, was also connected 
to the pressure transducer readout to obtain total adsorbate pressure as a function of 
time. 

Vapor introduction procedures were the same as those described in section 2. 1. 2, 
except that in no case was the microbalance system pumped by cracking the gate valve. 
This limitation was necessary in order to make certain that all the first species 
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Table 3-1 


Run No. 
F9-1 

F9-2 

FI 1-4 

Fll-7 

FI 1-6 

F6-1 

F6-2 

F7-2 

F7-3 

F7-7 

F12-1 

F17-3 

F17-4 


SUMMARY OF FLOW SYSTEM MIXED ADSORPTION 
DATA - BLOCKAGE EFFECTS 


Contaminant 

Temperature 
(°C) 

Concentration 

(ppm) 

Blockage 

(%) 

n-Butane 

25 

20.4 

V 

Propylene 


0.54 

86 ± 2 

n-Butane 

25 

1.08 


Propylene 


13.2 

29 ±7 

n-Butane 

27 

0.9 


Propylene 


1580 

15 ± 8 

n-Butane 

25 

19.5 


Propylene 


13.4 

69 ±3 

Propylene 

25 

13.2 


n-Butane 


20.2 

0 

Propylene 

25 

21.0 


Vinylidene Chloride 


18.3 

0 

Propylene 

46 

20.2 


Vinylidene Chloride 


25.8 

0 

Vinylidene Chloride 

102 

25.1 


Toluene 


29.5 

6 ±9 

Toluene 

102 

24.8 


Vinylidene Chloride 


25.2 

91 

Toluene 

102 

6.1 


Vinlyidene Chloride 


23.5 

93 

Propylene 

25 

13.8 


Ethylene 


27.2 

0 

n-Butane 

107 

26.9 

30 

Toluene 


26.8 

Toluene 

67 

26.5 


n-Butane 


46.0 

96 

Freon- 12 


38.5 
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•* produced was still in the system, except for an insignificant amount that was pumped 
out for analysis through the mass spectrometer section. 

Un determining these mixed-adsorption isotherms, the more strongly adsorbed con 
4aminant was admitted first until the weight gain of the charcoal corresponded to a 
predetermined amount adsorbed, usually 1,2,5, 10, or 20 ml STP/g. As equilibrium 
tfas being attained, both adsorbate pressure and characteristic mass spectral peak 
Wight were noted. Next, the less strongly adsorbed species was admitted, after 
{lushing and pumping the vapor admission manifold with the new vapor. After each 
Admission of the less strongly adsorbed contaminant, total pressure and the mass 
Spectral peak heights of both species were noted. 

-The peak heights corresponding to the first, more strongly adsorbed contaminant were 
Checked to see if they were remaining constant in order to assure that this contaminant 
Vvas not being displaced by the less strongly adsorbed one. Partial pressure of the 
Second, less strongly adsorbed contaminant was obtained by subtracting the initial 
pressure of the first, more strongly adsorbed one from the total at each point in the 
Isotherm determination, while weight gain corresponding to the second contaminant 
was obtained by subtracting the weight after adsorption of the first one. Both calcu- 
lations were based on the assumption of no displacement as determined above. 

Examples of mass spectra and pressure records and of weight records for n-butane on 
charcoal with preadsorbed toluene are given in Appendix B. 

1.2 RESULTS 

3.2.1 Flow System 

Results of multiple -component flow adsorption experiments are included in the experi- 
mental summary in Appendix C and are summarized in Table 3-1. This table is 
broken into five sections corresponding to the five contaminant pairs studied: n-butane 
propylene; propylene-vinylidene chloride; vinylidene chloride -toluene , propylene 
ethylene ; and n-butane -toluene. The last group also includes an experiment in which 
toluene was preadsorbed, followed by concurrent adsorption of both n-butane and 
Freon-12. 
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3.2.2 Gravimetric System 


Isotherm plots showing adsorption of n-butane at 27° and 67 °C with various amounts : .. 
of preadsorbed toluene are presented in Figs. 3-1 and 3-2, respectively. Data are 
for 5, 10, and 50 ml STP/g of preadsorbed toluene. In the 50 ml STP/g cases, some 
displacement of toluene was evidenced by increases in mass spectral peak intensities 
as n-butane pressure was increased; thus , the butane adsorption shown in these 
isotherms may be misleadingly low. Virtually no displacement of toluene was experi- 
enced at 5 or 10 ml STP/g toluene preadsorbed at either temperature. 

In Table 3-2 blockage of n-butane by toluene when n-butane is at MAC pressure 
(0. 02 torr) is shown as a function of amount of preadsorbed toluene. Also shown is 
the partial pressure of toluene needed to obtain this amount adsorbed. As can be 
seen, the equilibrium pressures of toluene are far below the MAC value, 0.02 torr , 
for blockages of n-butane of well over 50 percent. At room temperature (27° C), 
adsorption of 50 ml STP/g of toluene at only one-third of MAC partial pressure results 
in 93 percent blockage of n-butane. Thus, a sorbent bed would have to be approximately 
14 times as large as predicted for MAC n-butane if toluene were also present, even at 
this reduced level. 


Table 3-2 

DATA FOR BLOCKAGE OF n-BUTANE BY TOLUENE ON 8 x 12 BD CHARCOAL 


Amount of 
Toluene 

Toluene Equilibrium 
Pressure 

Blockage of n-Butane 
’ at MAC (0.02 torr) 

Preadsorbed 

27 °C 

67° C 

27° C 

67° C 

(ml STP/g) 

( torr ) 

( torr ) 

<%) 

(%) 

5 

< 1 x 10“ 5 

~ 8 x 10' 5 

29 

47 

10 

< 1 x 10" 4 

5.6 x 10“ 4 

51 

65 

50 

7.2 x 10" 3 

9.1 x 10 -2 

93 (a) 

ioo (a) 


(a) Possibly somewhat less, due to displacement of toluene by n-butane. 
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AMOUNT ADSORBED (ml STP/g) 



PRESSURE (torr) 


Fig. 3-1 Isotherms of n-Butane on 8 x 12 BD Charcoal at 27°C: Blockage by Toluene 
at Various Amounts Preadsorbed 
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AMOUNT ADSORBED (ml STP/g) 



PRESSURE (torr) 


Fig. 3-2 Isotherms of n-Butane on 8 x 12 BD Charcoal at 67°C: Blockage by Toluene 
at Various Amounts Preadsorbed 
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order to compare these gravimetric data with the flow system data in Table 3-1, 
blockage effects at n-butane partial pressures equivalent to 26.3 ppm at 0. 5 atm 
-total pressure are listed in Table 3-3. Thus, these are the blockage effects that 
Would be expected in a flow experiment in which toluene and n-butane were present m 

•tVe concentrations indicated. 


Table 3-3 


GRAVIMETRIC SYSTEM MIXED ADSORPTION BLOCKAGE DATA 
EXPRESSED IN CONCENTRATION TERMS 


Run No. 


Temperature 
Contaminant (° c ) 


Concentration Blockage 

(ppm) (%) 


/YI20-22 to -30 
-74 to -81 
-42 to -49 


-65 to -73 


-31 to -38 


-50 to -56 


Toluene 

n-Butane 

Toluene 

n-Butane 

Toluene 

n-Butane 

Toluene 

n-Butane 

Toluene 

n-Butane 

Toluene 
n- Butane 


27 

27 

27 

67 

67 

67 


0.01 

26.3 

39 

0.1 

26.3 

58 

19 

26.3 

97 

~0.2 

26.3 

57 

1.5 

26.3 

73 

240 

26.3 

>99 
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Section 4 

DISCUSSION AND THEORETICAL TREATMENT OF ADSORPTION DATA 


It is essential, because of the wide variety of compounds known even now to be present 
as trace contaminants in spacecraft environments, to provide some rational basis for 
the estimation of expected levels of sorbent loading. While it may not be possible 
to produce a correlation that is quantitatively accurate and of high precision over the 
five or six orders of magnitude that are of interest, it is necessary that the correlation 
be quantitative enough for design purposes and that the physical-chemical data required 
to apply the correlation be readily available. The so-called "potential theory" of 
adsorption seems to meet these rather stringent requirements suprisingly well. In 
this section it is demonstrated that the amount adsorbed as a function of pressure and 
temperature may be estimated within engineering accuracy merely from molar volume 
of the liquid at its boiling point and vapor pressure of the pure liquid. The theory is 
also generalized to the case of multicomponent adsorption. 

4. 1 POTENTIAL PLOT CORRELATION 

The potential theory of adsorption was originally proposed by Polanyi (Refs. 1, 2, 3) 
and later modified by others (Refs. 4 and 5). It was first used to correlate adsorption 
data collected under spacecraft conditions by Robell et al. (Ref. 6). The correlation, 
called a potential plot, enables estimation of q , the amount adsorbed, from the 
adsorption potential: 


A = 




(4.1) 


where 

A = adsorption potential (mol °K/ml liq) 
T = temperature (°K) 
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V = contaminant molar volume at normal boiling point (ml liq/mol) 
m 

P° = contaminant vapor pressure (atm) 
p = contaminant partial pressure (atm) 

The adsorption potential can be calculated from two basic physical properties - molar 
volume and vapor pressure. These properties are listed, for a few selected trace 
contaminants considered in this work, along with estimated values for molecular area 
and monolayer coverage, in Table 4-1. 

Table 4-1 

SELECTED PHYSICAL PROPERTIES OF TRACE CONTAMINANTS 



Molar 

Vapor 

Molecular 

Monolayer Coverage 


Volume 

Pressure 

Area 

<1™ 


at 27°C 




L'OllWIlilIiw'i* 

Vjn 

o i 8 


(ml liq/g) 


(ml liq/mol) 

P° (atm) 

(A ^/molecule) 

(ml STP/g) 

Toluene 

118.2 

0.042 

37.4 

109.4 

0.578 

n-Butane 

96.4 

2.55 

31.9 

128.3 

.552 

Vinylidene 

Chloride 

79.7 

.84 

28.6 

143.1 

.510 

Ethanol 

62. 1 

.084 

25.7 

159.2 

.441 

Propane 

74.5 

9.8 

27.1 

151.0 

.502 

Propylene 

66.6 

12.0 

25.4 

161. 1 

.479 

Ethylene 

49.4 

42.5 

20.2 

203.0 

.447 

Dichlorodi- 
fluoromethane 
(Freon- 12) 

80.7 

6.9 

28.9 

141.6 

.510 

Water 

18.80 

.032 

13.1 

311 

.262 


Molecular area, /3 , has been shown by Hill (Ref. 7) to be equivalent to the two 
dimensional van der Waals constant, which can be evaluated by the equation: 


b = 6.354 x 10" 16 (T c /P c ) 2/3 (4-2) 


33 




where 


b -j two-dimensional van der Waals constant (cm 2 /molecule) 

T c = critical temperature (°K) 
p c = critical pressure (atm) 

Values of b were calculated from critical temperature and pressure data of Robell 
(Ref. 8). These values are to be used only as approximations in surface coverage 
calculations, since they are quite possibly smaller than actual effective molecular 
areas. For instance, Davis (Ref. 9) and Emmett (Ref. 10) have shown that n-butane 
molecular area ranges from 43 to 52 A 2 /molecule depending on the sorbent. Since 
similar values for other contaminants studied are not readily available, the b-values 
provide a consistent means of estimating relative coverages of these contaminants. 
Values of q m were calculated from the equation 

q n, * = S s v m /W (4.3) 

monolayer coverage (ml liq/g) 

sorbent surface area, assumed = 1 . 1 x 10^ (cm 2 /g) 
standard state molar volume (ml STP/mol) 

Avogadro's number = 6. 02 x 10 23 (molecules/mol) * 

The derivation of the potential plot correlation is semiempirical and has previously 
been discussed thoroughly (Ref. 5, 6, 11). 

A wealth of adsorption data on Barnebey-Che-ney type BD activated carbon have been 
obtained in this and other programs. 


where 


q m = 

S s = 

V = 
o 

N = 


Some of these data are plotted in Fig. 4-1 as a potential plot. The correlation is quite 
good and extends over four orders of magnitude in amounts adsorbed. It is striking 
how consistent the results are, even though the data were obtained on the three 
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ADSORPTION POTENTIAL (mol°K/ml liq) 


Fig. 4-1 Potential Theory Plot for BD Activated Carbon 


separate experimental systems indicated in the figure legend. The correlation is by 
no means perfect; there are deviations of the order of 30%. However, the potential 

plot is still the best correlation found to date. It is recommended that Fig. 4-1 be 
used for design purposes. 

The equation shown in Fig. 4-1 delineates the straight-line portion of the plot which, 
for this work, is the region of interest. All potential plots appearing elsewhere in 
this report include the line from Fig. 4-1 for reference purposes. 

A more detailed exposition of the potential plot correlation for one adsorbate is 
shown in Fig. 4-2 and 4-3. The dotted line represents the potential plot as determined 
m Fig. 4-1. The n-butane data shown in Fig. 4-2 fit the correlation admirably, 
especially considering the temperature range studied and the fact that the data were 
obtained from three different experimental systems. (Run numbers beginning with M 
refer to microbalance gravimetric experiments; those beginning with F refer to flow 
system experiments. ) 

The toluene data in Fig. 4-3 are quite consistently above the potential plot line. There 
is no explanation for this variance at present, although the correlation still holds to 
about 30% as mentioned earlier, and this is adequate for design purposes. 

These potential plots were made for single contaminants present in pure form (micro- 
balance runs) and in nitrogen-oxygen mixtures (flow runs). There is no apparent 
difference m adsorptive capacity between these two cases. However, in the case of 

mixtures or when humidity is present, adsorptive capacity is drastically affected, as 
discussed later (Section 9). 

A new derivation of the Polanyi potential plot is presented in the following paragraphs. 
This derivation starts from fundamental principles and goes through to the convenient 
calculation of thermodynamic quantities, such as heat of adsorption, and of kinetics 
quantities, such as desorption rate order. 
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AMOUNT ADSORBED (ml liq/g) 


Temp. 

Particle Size 


(°C) 

(Mesh) dp(cm) 

Run No 


U 27 

8 x 12 

0.20 

M20-1, -19. -22, -42, -4 

67 

8 x 12 

0.20 

M20-31, -50, -65; F17-4' 

k 107 

8 x 12 

0.20 

F17-1 

□ 27 

16 x 20 

0.10 

M18-28 to -37, -69 to -74 

O 67 

16 x 20 

0.10 

M18-38 to -47, -61 to -68 

0 107 

16 x 20 

0.10 

M18-49 to -59 

O 27 

30 x 40 

0.05 

M9-1 ; M10-2 ; M19-25 

O 87 

30 x 40 

0.05 

M13-1 to -3 

^ 102 

30 x 40 

0.05 

F7-3, -7 

Vl07 

30 x 40 

0.05 

M19-16 to -22 


1 ml liq/g = 189.5 ml STP/g 


1. 

. 732 Monolayers 




ADSORPTION POTENTIAL (mol°K/ml liq) 


Fig. 4-3 Potential Plot for Adsorption of Toluene on BD Activated Carbon 




The general form of the Polanyi isotherm can 


be derived from the Gibbs equation: 


P 

7 T = (y - y s ) = RT J rd In f 


(4.4) 


where 

2 

7 T = spreading pressure (cal/cm ) 

y = surface free energy of adsorbed layer per unit area (cal/cm ) 

v = clean surface free energy per unit area (cal/cm ) 
s 

R = gas constant (cal/mol °K) 

2 

p = pressure (cal/cm ) 

r = surface excess concentration (mol /cm ) 

3 

f = fugacity of adsorbate (cal/cm ) 

and two assumptions that are quite reasonable for very low coverages. When this is 
done the derivation gives exactly the linear form of the isotherm, which is so often 
observed at very low coverages. Firstly, it is presumed that the surface is composed 
of portions that are essentially bare and portions covered with absorbate. Secondly, 
it is assumed that the bare portion has a surface free energy per unit area identical 
to the clean surface free energy, Y g , and that the surface free energy per unit area 
of the portion of the surface covered with adsorbate is a constant, y A , independent 
of coverage for portions of the surface covered only by a monolayer of adsorbate. 

The surface free energy of the adsorbed phase, y , may then be written: 

y = ey A + (1 - 0)Y S (4 ‘ 5) 


where 6 = fractional coverage. Thus, 


tt = 0 ( y A 


Y s> 


(4. 6) 


Equation 4.4 may now be differentiated to give: 

d7r = (y A - y g )d0 = FRTdlnf 
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but 


e = pr 


(4.8) 


where 0 - monolayer area in cm 2 /molecule. 


Equation 4.7 may be rearranged to give: 


d In T 


HT 

p(y A - yj 


In f 


(4.9) 


Hie quantity <v A - y g ) in Eq. 4. 9 is the excess surface free energy per unit area 
associated with the adsorbate. There is no a priori reason to expect (y A - y ) to be 
constant in the multilayer region; nevertheless, Eq. 4.8 maybe, in principle? 
integrated even for multilayer adsorption to give the isotherm; 


in r/r° = rt 


i 

/ 


din f 

0 (y A - r s ) 


(4.10) 


where r and f are limits suitably chosen for the integration. 

At low coverages, as mentioned earlier, it is assumed that y and y Q are constant 
and independent of coverage. Then the denominator in Eq. 4. 10 is constant, and 
integration is straightforward. 


Following the original idea of Polanyi (Ref. 1) or taking guidance from the success 
of regular solution theory, it is more likely that the surface free energy per unit 
volume would be independent of the adsorbate on a particular adsorbent than the free 
energy per molecule. Thus, it might be reasonable to expect that, for a given class 
of adsorbates (say a homologous series of hydrocarbons), the quantity (y A - y )/6 
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where 6 m is the molecular diameter, should be relatively constant. Thus, Eq. 4.9 

may be rewritten: 


in r/r° 


RT 


} 

po 


<5 d In f 
m 


(4. 11) 


and may be integrated to give: 


, r 
ln-^ 


RT 

V m 


m 


(yn ~ y*> 


In f° /f 


(4. 12) 


which is identical to the empirical Coverage limit of the Poltbiyi isotherm Il f 
is tahen as the fugacity of samrated adsorbate liquid at the temperature of the adsorption. 
Of course, I~ has no real physical significance under such circumstances, and 

Eq 4. 11 departs from the empirical isotherm as soon as y A begins to change With 

coverage. Nevertheless, a rather simple physical interpretation results from t e 
slope of the linear portion of the Polanyi isotherm. It is the reciprocal of the excess 
surface free energy per unit volume of adsorbate. 

A number of thermodynamic properties can be derived simply from Eq. 4. 12. The 
isosteric heat, q g ^. , is given by 


*st 


= - R 



(4. 13) 


but the fugacity usually can be approximated by the pressure at low coverages where 
Eq. 4. 12 is appropriate. Thus: 


q st - AH V - i* <? A - ln r/r " 


(4. 14) 


m 
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This has the same tactical form as the relation that Halsey and Taylor (Ret 121 hav 

derived for the Freundlich isotherm and indeed Eo 4 „ „ ' ’ aVe 

’ maeed - E q. 4.11 can be rearranged in the 


Freundlich form : 


if we set: 


and 


r = k f p 


1/n 


(4.15) 


k = r° [p° ] ~ 1//n 


(4.16) 


n . VA - V 

6 RT 
m 


(4.17) 


While Halsey and Taylor (Ref. 12) derived their 


results for a heterogeneous surface, 


this development shows that the Freundlich isotherm and the logarithmic dependence 
o e isosteric heat(Eq.4. 14) can come as well from a homogeneous surface. 


The integral enthalpy of adsorption, AH 


a , is given (Ref. 13) as: 


AH = H - H = q + -l/chr) 

a a g 4 st r \at/. 


where 


H^ adsorbed-phase enthalpy (cal/mol°K) 
H g = gas-phase (cal/mol°K) 

using Eq. 4.6 and 4.14 


(4.18) 


- - 1st + T ^[|<r A - y] r 

and under the assumptions consistent with Eq. 4.12, Eq. 4.19 becomes: 


(4.19) 


AH = - q + TV 

a 4 st m vat/r 
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(4.20) 




R 

<P 

<P 


negative slope of Polanyi plot (ml liq/mol °K) 

T A 7 s excess f ree energy per unit volume (cal/ml liq) 
6 

m 


The integral entropy change is AH & /T . In fact, the isotherm results in a very 
SWsple equation of state for the adsorbed layer: 


which any of the desired thermodynamic properties may be calculated by well' 
J^jown methods. 


Usually <p is independent of temperature or nearly so in which case 

q . « AH = TAS 
H st a a 


(4.22) 


^Kere AS g = entropy of adsorption in cal/mol °K. 

^firie temperature dependence of <p is always quite small when Eq. 4. 12 is appropri 
it all, and reducing experimental data according to Eq. 4.12, and/or Eq. 4.21 sub- 
stantially facilitates the extraction of thermodynamic properties in contrast to other 
Approaches. These ideas can be applied to the nonlinear portion of the isotherm if 
*4- is recognized that outside the linear range one may use for (p 



(4.23) 


VJhere 

— Polanyi adsorption potential (ml liq/mol ( K) 

A = f^ln f°/f 
v m 


(4. 23a) 
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4.2 CALCULATION OF HEAT AND ENTROPY OF ADSORPTION FOR BUTANE AND 
TOLUENE 

A good deal of information about the thermodynamics of an adsorbed gas can be 
obtained from a set of isotherms such as those shown in section 2.2. These 
thermodynamic properties can then be compared with those calculated for pertinent 
idealized models of the adsorbed layer. It is shown in the following paragraphs that 
the picture which emerges is consistent with the desorption mechanism deduced from 
desorption rate data. For example, it is shown how entropy data as a function of 
coverage establish the relative surface mobility of adsorbates and how enthalpy data 
contribute to characterizing surface heterogeneity. In this section these methods 
are applied in detail to the adsorption of n-butane and toluene on BD charcoal. 

Several heats of adsorption may be defined, but the one most easily identified with the 
adsorption isotherm is the isosteric heat of adsorption, , which may be defined 
(Ref. 14) as the heat released in a constant temperature calorimeter when a differential 
amount of gas is adsorbed at constant pressure. 

1st = < 3<5/3,I s>T,p < 4 - 3 °) 


where 

q gt = isosteric heat of adsorption (cal/mol) 

Q = reversible heat change of the system (cal) 

n = number of moles of gas adsorbed (mol) 
s 

T = absolute temperature (°K) 

3 

p = partial pressure of the adsorbate in the gas phase (cal/cm ) 

The isosteric heat may be calculated from the adsorption isotherms by applying the 
following equation (Ref. 14): 


q st = - R[ainp/9(1/T)] 0 


(4.31) 


W here 


R = gas constant = 1. 98 (cal/mol K) 

9 = fraction of the surface covered with adsorbed molecules 


T* Fig 4-4 logarithm of the partial pressure of n-butane is plotted as a function of 
reciprocal temperature for several ualues - 

3* Fig. 4-5 data for more recent experiments with 0 5 cat* P ^ ^ 

vOhile data for 1- and 2-mm carbon appear in Fig. 4-6 ’ P 

o( these isosteres is , S ,/R . The values of %t as deduced from Fig. 

tabulated in Table 4-2. 

A similar plot for toluene adsorption on 1-mm particles appears in Fig. 4-8 and values 
of q are tabulated in Table 4-3. 

S v 

isosteric heat can be shown (Ref. 13) to be related to the integral molar enthalpy 
by the following equation: 


AH == H - H = - q t + ( T ^/ 0 H a,r / 8T) 0 

si a £ 


(4.32) 


Where 

AH = integral molar enthalpy of adsorption (cal/mol) 

H a = enthalpy of the adsorbate in the gas phase (cal/mol) 

H S = enthalpy of the adsorbate in the adsorbed phase (cal/mol) 

a 2 
/3 = molecular area (cm ) 

7 T = spreading pressure (cal/cm ) 


It is shown in the previous section that systems 
theory have a very simple equation of state; 1 . e 


which follow the Polanyi potential 


n 


v < pe 

m 


P 


(4.33) 
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Fig. 4 4 Estimation of Isosteric Heat of Adsorption for n-Butane on BD 
Charcoal at Various Coverages 
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PRESSURE (torr) 



Fie 4-5 Estimation of Isosteric Heat of Adsorption of n-Butane on 30 x 40 
BD Activated Carbon (d p = 0. 05 cm) at Various Amounts Adsorbed 
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reciprocal TEMPERATURE (°K x ) 


4-7 Estimation of Isosteric Heat of Adsorption of n-Butane on 8 x 12 BD 
Activated Carbon (d = 0.20 cm) at Various Amounts Adsorbed 









Table 4-2 


HEAT OF ADSORPTION OF n-BUTANE ON BD ACTIVATED CARBON 


Particle Size 
(cm) 

Amount Adsorbed 
(ml STP/g) 

Fractional Coverage 

e 

Isosteric Heat 
of Adsorption 
Qst 

(kc^l/mol) 

0.05 (a) 

8. 0 

0.0624 

12.3^ 


4.0 

.0312 

12. 0^ 


2.0 

.0156 

12.0 (a) 


1.0 

.0078 

12.5 (a) 


0.5 

.0039 

13.4 (a) 

0,05 

10.0 

0.078 

13.2 


5.0 

.039 

13.9 


2.0 

.0156 

14.2 

0.1 

10.0 

0.078 

12.5 


5.0 

.039 

13.6 


2.0 

.0156 

13.8 

0.2 

50.0 

0.39 

10.8 


20.0 

.156 

11.8 


10.0 

.078 

11.6 


5.0 

.039 

11.2 


2.0 

.0156 

11.4 


1.0 

.0078 

10.6 


(a) Early data based on Pirani-gage pressure readings 


Table 4-3 

THERMODYNAMIC PROPERTIES OF TOLUENE ADSORBED ON 1. 0-nrai BD 
ACTIVATED CARBON (16 * 20 MESH) AT 27 °C (300 °K) 


Amount 

Adsorbed 

q 

(ml STP/g) 

Fractional 

Coverage 

e 

Isosteric Heat 
of Adsorption 

(kcal/mol) 

50 

0.457 

15, 1 

20 

.183 

15.7 

10 

.0914 

17 

5 

. 0457 ' 

17 
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where (p is the excess free energy per unit volume of adsorbate. This excess free 
energy is given by: 


4 > 


R 


(W 


-1 


(4.34) 


When the potential plot is linear, 0 is a constant. Even when <p varies with coverage, 
it is usually nearly independent of coverage so that 



« 0 


and then 


’at " AH a - TAS a 


(4. 35) 


For the n butane and toluene adsorption, however, there is some curvature in the 
potential plot (see Fig. 4-2 and 4-3). The temperature dependence of <p is small 
enough to permit the use of Eq. 4.35. 

In Fig. 4-9 the values of AH^ for butane are plotted as a function of- coverage along 
with data reported by de Boer and Kruyer (Ref. 15). Although there is a good deal of 
scatter in the data, it is clear that the isosteric heats decrease more or less linearly 
with coverage, and that the data obtained for the very low coverage of interest for 
contaminant control are in essential agreement with previously reported data at much 
higher coverages. The decrease in the value of q t results from the energetic 
heterogeneity of the carbon surface. 

When a system obeys the Polanyi equation (Eq. 4.1), Eq. 4.31 has a very simple 
form: 


q 


st 


AH - V <p In « AH 
v nrr q a 


(4. 36) 
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ADSORPTION ENTHALPY, AH a « q . (kcal/mol) 



Experimental (q . ) Data of de Boer and Kruyer 
Sl (Ref. 15) for Charcoal 

O 30 x 40 BD at 353°K 0 

380°K O 

□ 16 x 20 BD 408°K A 

436°K ^ 

q 8 x 12 BD 464 °K [7 

-•-Calculated From Eq. 4.36 




FRACTIONAL COVERAGE 


Fig. 4-9 


omparison of Experimental n-Butane Adsorption Enthalpy With Literature 
a W 44 -K uoi.ioo r.al mil atsd From " Pntfiuti a 1-Plot Correlc 


o Correlation 


of Data for BD Activated Carbon 



The solid line in Fig. 4 9 is the adsorption enthalpy calculated from Eq. 4.36 and 
the potential plot of the n-butane data. It fits the scattered data rather well and is 
consistent with the de Boer and Kruyer data. This method for deriving the thermo- 
dyanmic properties of adsorbed layers seems to be much preferred over the direct 
application of Eq. 4. 31 on the individual isotherms. 

The entropies of adsorption consistent with the solid line in Fig. 4-9 are plotted in 
Fig. 4-10. It is instructive to compare theoretical entropies for idealized models with 
these experimentally derived values. From such a comparison it is often possible to 
infer something about the state of the adsorbed molecules on the surface. 

DeBoer and Kruyer (Ref. 16) have outlined the calculation of standard- state adsorptive 
entropies for both mobile and immobile adsorbed layers. In the case of mobile layers, 
they select the following for standard states: one atmosphere in the gas phase and an 
adsorbed layer with an average molecular separation equal to that in the gas -phase 

standard state. Thus, if only the entropy associated with translation is considered, 
then 


where 


s g = s g - Rlnp 


Sg - entropy of the gas at adsorption partial pressure (cal/mol°K) 
= standard-state entropy of the gas (cal/mol°K) 
p = adsorption partial pressure (atm) 

and 


(4. 37) 


S = 
a 



(4. 38) 


where 

~ entropy of adsorbate at actual coverage (cal/mol °K) 
= standard- state entropy of adsorbate (cal/mol °K) 
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ADSORPTION ENTROPY, -AS a (cal/moI°K) 





Experimental Values, 
(q^From Eq. 4.36) 

Idealized Model (Eq. 4.41) 



Fig 4-10 Comparison of n-Butane Adsorption Entropy Derived From a Smoothed 
S ' Potential-Plot Correlation of Data With an Idealized Model for the 

Adsorbed Layer at 300 °K 



°o ~ standard-state molecular area 
= 4. 08 T x 10 16 (cm 2 /molecule) 

= monola y er molecular area (see Table 4-1) (cm 2 /molecule) 
0 = fractional coverage 

P S and s ° are given in Ref. 16 as- 

5 

S g = Rln (M 3/2 T 5/2 ) - 2.30 


and 


S 


o 

a 


§S° 
3 S 



+ 5 R1 "2705 


Thus, the adsorption entropy may be calculated by 


(4. 40) 


AS a S a " S g 3 R ln (M 3y/2 T 5j/2 ) - 2.28 + 0.662 In T - R In ~ 


P 


+ R ln P + R In 0 


(4.41) 


The theoretical entropies for a completely mobile adsorbed layer of n-butane calculated 
from Eq. 4.41 are plotted as the solid line in Fig. 4-10. These calculations are 
given in Appendix I. They compare quite well with the measured values, suggesting 
that the absorbed layer is indeed very mobile. It will be shown later that this result 
is in agreement with the kinetics of desorption from single particles. 


The toluene adsorption data may be treated in an analogous manner. The isosteric 
heats (integral enthalpies) are plotted in Fig. 4-11. These calculations are also given 
m Appendix I. Again the heats drop with coverage, and the values derived from the 
potential plot shown by the solid line in Fig. 4-11 are in agreement with the values 
obtained directly from the isotherms. 
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ADSORPTION ENTHALPY, AH * q . (kcal/mol) 



FRACTIONAL COVERAGE, 0 

Fig. 4-11 Comparison of Experimental Toluene Adsorption Enthalpy With Values 
Calculated From a Smooth Potential- Plot Correlation of Data for 
8 x 12 BD Activated Carbon (d = 0. 20 cm) 



In Fig. 4-12 the entropies are plotted and compared with Eq. 4.41. Again the 
agreement is nearly quantitative and it may also be concluded that absorbed toluene, 
like n— butane, is mobile. It is difficult to say with certainty whether all contaminants 
will be mobile, but it is certainly not unreasonable to assume that all components 
lighter than n-butane and toluene are mobile except perhaps for materials which are 
highly polar, exhibit marked hydrogen bonding, or for other reasons have very strong 
molecule-molecule interactions. 

4.3 FIT OF TOLUENE AND BUTANE ISOTHERMS TO EMPIRICAL EQUATIONS 

It is of some interest to fit the adsorption isotherms of toluene and of butane obtained 
under this contract on BD activated carbon to empirical equations. The toluene 
isotherms of Fig. 2-10 are fitted by the equation 

16 300 - 57.6 q 

log 10 P = 2.303 RT + 5.81 + log 1() q & (4.42) 


where 

p = partial pressure of adsorbate (torr) 

R = gas constant = 1.987 (cal/mol°K) 
q a = amount adsorbed (ml STP/g) 

The above fit was obtained from the empirical observation that the isotherms when 
plotted as RT log p/q & versus q & gave reasonably straight lines, and the plots for 
different temperatures were roughly parallel. The quantity 

16 300 - 57.6 q a (4.43) 

is the isosteric heat of adsorption derived from the above plots by the Clausius- 
Clapeyron equation, and this value agrees well with the values derived in section 4. 2. 
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ADSORPTION ENTROPY, -AS o (cal/mol°K) 

a 



FRACTIONAL COVERAGE, 6 


Fig, 4-12 Comparison of Toluene Adsorption Entropy Derived From a Smoothed Potential- 
Plot Correlation of Data With an Idealized Model for the Adsorbed Layer at 300° K 



The older butane isotherms shown in Fig. 2-8 are fitted by: 


12 500 - 100 q , 

log 10 P ~ 2.303 RT + 8 . 0 - 2. 6 q & (4. 44) 

Here 12 500 - 100 q & is the isosteric heat of adsorption (cal/mol), which is a smoothe- 
fit to the data given in section 4. 2. The entropy term 8. 0—2. 6 q -1//5 was thei derived 

3 . 

to give a good fit to the isotherm data. 

4.4 CORRELATION AND PREDICTION OF MIXED ISOTHERM DATA 


An important factor in the proper design of a space-cabin contaminant adsorption bed 
is the prediction and correlation of adsorber performance in the presence of many 
contaminants. In Section 3 the experimental adsorption isotherms for several con- 
taminant mixtures have been presented. In this section it is shown how these may be 
predicted and/or correlated with some accuracy. 

Various theories of mixed adsorption have been proposed and the work up to 1962 has 
been summarized by Young and Crowell (Ref. 14, pp. 371-91). More recent work 

includes that of Myers and Prausnitz (Ref. 17), Grant and Manes (Ref. 18), and Loven 
(Ref. 19). 

The main objective of a theory of mixed adsorption is to predict the adsorption 
isotherms of the mixture from the isotherms of the pure components. The equation 
which appears to be most general and empirically accurate is that proposed by Grant 
and Manes and by Loven: 


P ml = X lP 0 l (q t> 


(4.45) 


where 


P m ^ - equilibrium partial pressure of component 1 in the gas mixture (torr) 
= mole fraction of component 1 in the adsorbed phase 
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p (, ) . pressure of pure component 1 (at the same temperature as the 
01 ‘ mixture) which would be in equilibrium with a quantity of adsorbe 

gas q T equal to the sum total of all the adsorbed gases in the mix- 

ture (torr) 

Equation 4 45 sufficies to compute the adsorption isotherm for a mixture of many 
Ration 4 45^ of fte adsorbed phas e is fixed at an amount q. 

Components. Thus lUhe c p ^ , ften these are summed to give 

J ? - component f in the adsorbed phase is computed as 
ratio of q : to the sum of the q's (q in ml STP/g). 

F to, - isodierm of pure component I 
Ayv amount q of component 1 is read off. 

ponent 1 , p^ , in equilibrium with this presumed adsorbed mixture (q, • V 
| 5 calculated. 

f ~ t io and 4-14 were computed. The agreement 

clarify the method of pr^lrbed. 5.8 ml STP/g of 

a is at 27° in which, with 5 ml STP/g ol toluene pie . 

pig. 4-13 at 27 m w -2 To CO mpute this point, 

butane is adsorbed at a butane pressure of 1.78 1 

. , that q (amount of butane adsorbed) is 5.8 ml STP/g T 2 

-the assumption is that q x famoui =58 + 5.0 = 

(amount of toluene adsorbed) is 5 ml STP/g, so that q t q i q 2 * = 

.0 8 ml STP/g and X, - mole fraction of butane in the adsorbed phase - 5.8/10.8 

to. 8 ml blP/g ana 1 10 8 ml STP/g is now 

0.537. The pressure of pure butane in equi i^ rium ^ us / the desire d butane 

read from the pure butane iso erm as S TP/g can be immediately 

pressure in equilibrium with q x - 5.8, q 2 

computed from Eq. 4.45. 

p = 0.537 X 2.8 x 10 2 = 1.5 x 10 torr 

F ml 
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AMOUNT ADSORBED (ml STP/g) 



Fig. 4-13 Blockage of n-Butane by Toluene on 8 x 12 BD Charcoal at 27°C: 
Comparison With Blockage Predicted From n-Butane Isotherm 
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n-BUTANE PRESSURE (torr) 


Fig. 4-14 


Blockage of 
Comparison 


n- Butane by Toluene on 8 x 12 BD Charcoal at 67 °C 
with Blockage Predicted From n- Butane Isotherm 



which thus locates the desired point and is within probable experimental error equal 
to the measured value , 1.78 x lo” 2 torr. 


The origin and meaning of Eq. 4.45 are of some interest. Grant and Manes, and 
Loven appear to have intuitively assumed it as a logical extension of the Polanyi 
potential theory. In Appendix J, it is shown that 

• Equation 4.45 may be derived by a simple kinetic argument. 

• The Langmuir theory of mixed adsorption, and an equation due to Glueckauf 
(Ref. 20) for mixed Freundlich-type adsorption on a more uniform surface 
are special cases of Eq. 4.45. 

The basic assumption of Eq. 4.43 is that the heat of adsorption depends only on the 
total adsorbed coverage, q T , and not on the exact composition of this adsorbed 
layer. This may be shown by differentiating both sides of Eq. 4. 43 (after taking 
logarithms) with respect to temperature to obtain the usual Clausius- Clapeyron-type 
of equation for q gt , the isosteric heat of adsorption. A second assumption implicit 
in Eq. 4.43 appears to be that the mixed adsorbed layer forms a perfect solution 
so that the partial molar entropy of mixing for component 1 is 


AS m " Rtax i 


(4.46) 


where 

A ^m = P ar tial molar entropy (cal/mol °K) 

R = gas constant = 1. 987 (cal/mol °K) 

Whatever its origin, not only this research but also the work of Loven (Ref. 19) and 

that of Grant and Manes (Ref. 18) show that it predicts and correlates a large number 
of mixed adsorption data. 

As shown above, Eq. 4.45 allows rapid computation of the mixed adsorption isotherms 
when the composition of the adsorbed phase is assumed known, and the unknowns are 
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the partial pressures of gases in equilibrium with this adsorbed composition. The 
converse case, when the partial pressures are known (and the desired quantities 
are the unknown amounts adsorbed) , becomes more unwieldy mathematically. Here 
there will be an equation of the same type as Eq. 4.45 for each component, and this 
equation must be solved simultaneously for the unknown q's, after fitting the isotherms 
for the pure components to a mathematical expression. 

Finally, it is easy to show that Eq. 4. 45 also predicts the correct blockage of pro- 
pylene by butane, as shown in Table 3-1. For runs F9-1, F9-2, and Fll-4 the 
adsorbed mole fraction of propylene may be predicted from Eq. 4.45, inverted to 
read 


X 


1 


ml 


Pol (q t> 


(4.47) 


and these computed values are compared with the experimental ones in Table 4 4. 

Table 4-4 

PREDICTED BLOCKAGE OF PROPYLENE BY n-BUTANE AT 27° C 


Run No. 

X-l- 

(Calculated) 

X^ (Experimental) 

F9-1 

0.54 

1000 

= 5.4 x 10 -4 

5.7 x 10~ 4 

F9-2 

13 

150 

= 8.7 x 10" 2 

1. 14 x 10 _1 

Fll-4 

1580 

1800 

= 8.8 x 10 _1 

8.7 x 10 _1 


The agreement is within die experimental error. In computing the predicted values 
of , an approximate propylene isotherm was derived from the points of Runs 
F9-3, Fll-5 , F6-1 (see Table 3-1). This isotherm can be represented by 
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(4.48) 


q = 2.13 x 10 -4 c 0,72 

d 

where 

q a = amount adsorbed at 27° C (ml liq/g) 
c = propylene gas phase concentration at a total pressure of 0. 5 atm (ppm) 

The mixed adsorption equation (Eq. 4. 43) can be incorporated in adsorption potential 
theory parameters as follows: 

In a mixture , the adsorption potential A T which correlates with the total amount 
adsorbed q^ is 

a t - t- lo « rfb < 4 - 49 > 

m *ol t' 

Substituting Eq. 4.47, 

T P ° X 1 

A T = — log — (4. 50) 

m F ml 

or 

A T = A + V" log X 1 ( 4 - 51 ) 

m 

where 

T P° 

A = V- log — (4, 52) 

m p ml 
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Thus, Eq. 4.48 or 4. 51 can be used to correlate mixed adsorption capacity data on 
a potential plot. 

This procedure is used to correlate the humid air data (see Section 9) and is also 
applied to several flow runs involving three different contaminant pairs , propane-n- 
butane, ethylene -propane , and n -butane -toluene. These data are summarized in 
Table 3-1 and plotted as suggested by Eq. 4. 51 in Fig. 4-15. All the data correlate 
well with the slope of the generalized potential plot (dotted line in Fig. 4-15) but 
all are somewhat higher than the generalized correlation. The deviations, however, 
are not appreciably greater than the scatter in the data used to establish the generalized 
correlation (see Fig. 4-3) , nor are they greater than the consistent deviation exhibited 
by the extensive toluene adsorption data shown in Fig. 2-10. Within the accuracy and 
precision of the generalized plot, this simple theory of mixed adsorption seems to be 
quite adequate and it is thus recommended for design purposes. 
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Fig. 4-15 Potential Plot for Mixed Adsorption 


TOTAL AMOUNT ADSORBED (ml liq/g) 



O Propane and n-Butane 
A Ethylene and Propane 






Section 5 

VACUUM DESORPTION FROM SINGLE PARTICLES — EXPERIMENTAL 

AND RESULTS 


Regeneration of sorbents after saturation with contaminants has been one of the major 
goals of this effort. In order to be able to make any sort of meaningful attempts at 
design of a spacecraft sorbent regeneration system involving the use of space vacuum, 
it is necessary to have extensive information on the effects of vacuum on the individual 
sorbent particles. Once the effects of temperature, particle size, and contaminant 
species have been established, it should be possible to relate these to the performance 
of packed beds under vacuum desorption conditions. 

With this end in mind, a series of experiments was performed on the gravimetric 
system described in Section 2. hi this series the above parameters were varied. 

In addition, experiments were performed, perhaps for the first time, in which more 
than one contaminant was monitored instantaneously during simultaneous vacuum 
desorption. This was made possible by using the gravimetric system with mass spec- 
trometer described in Sections 2 and 3. A few experiments were also performed on 
a volumetric apparatus in order to supplement the gravimetric experiments. 

5.1 EXPERIMENTAL TECHNIQUES 

5.1.1 Single-Component Experiments 

The apparatus used in vacuum desorption experiments consisted of the same vacuum 
microbalance system as described in section 2.1.2. Extremely high pumping speed 
was provided by the 400-liter /sec oil diffusion pump and baffled liquid -nitrogen trap 
located next to the 10-cm gate valve on the microbalance jar. 
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The procedure for adsorption of contaminants was the same as that described in 
section 2.1.2, except that the amounts adsorbed were controlled to certain nominal 
values by allowing weight to increase by a predetermined amount in each experiment. 

For the desorption portion of the experiments, the strip-chart recorder speed was set 
sufficiently high (2 or 4 in. /min) so that desorption time data could be taken to the 

nearest 0.5 sec. At time zero, the gate valve was quickly opened, and the ionization 
gage was turned on after about 5 sec. 

5.1.2 Multicomponent Experiments 

The apparatus used in multicomponent desorption experiments consisted of the system 
as described in sections 2.1.2, 3.1.2, and 5.1.1. The adsorption procedure was the 
same as for the first point on a multicomponent adsorption isotherm, with the equilibrium 
amount of the second, less strongly adsorbed contaminant somewhat greater than that 
of the first one. Desorption was initiated in the manner described in section 5. 1. 1 but 
with the mass spectrum /pressure strip-chart recorder running at higher speed. The 
ratio of desorption rates at any given instant was assumed to be equal to the ratio of 
partial pressures in the mass spectrometer as determined by the mass spectral peak 
heights corresponding to each species. Since the total amount of each species adsorbed 
was known, and assuming that only the more strongly adsorbed contaminant remained 
adsorbed in significant quantity after the desorption was 90 percent complete, it was 

possible to get instantaneous desorption rates for each species, as explained in 
section 5.2.2. 

5.1.3 Volumetric Method 

To supplement the gravimetric method for measuring vacuum desorption rates from 
single activated carbon granules, a volumetric desorption system was designed. In 
this method, the desorbed contaminant was caught in a cold trap (which also acted as 
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the pump) and was measured volumetrically after warming up the cold trap. The 
method has several advantages over the microbalance technique: 

(1) Freedom from thermogravimetric and other spurious effects, such as 
were encountered when gas was suddenly admitted to or evacuated from 
the microbalance 

(2) Great simplicity, since a large stopcock (or better, a solenoid valve), a 
suitable vacuum gage, and a vacuum pump are the only equipment needed 

(3) The capability of inserting a thermocouple directly into an activated carbon 
granule, so that adiabatic cooling effects can be directly measured 
simultaneously with desorption rate 

(4) The capability of varying pumping speed by simply lowering the liquid 
nitrogen level on the cold trap, which acts as the pump 

The sensitivities of the two methods are about equal, since a microgram of n-butane 
in 300 cm 3 volume results in a pressure of about 1(T 3 torr, which is easily detectable. 
A disadvantage of this method is that only one point on the desorption curve is measured 
in a given run, since the desorption process is interrupted to measure the amount of 
butane desorbed. However, each point takes only about 15 min, so that in a few hours 
the whole curve (amount desorbed versus time) can be established. 

The all-glass apparatus that was finally used after a little experimentation is shown in 
Fig 5-1. The activated carbon granules were supported on a fine mesh nickel screen 
in the vertical tube, T . The charcoal could be isolated from the rest of the system 
by the very large hollow-plug, single-hole, high-vacuum-type stopcock, S . The choice 
of S is critical, as an ordinary stopcock or valve will reduce the pumping speed too 
severely. Likewise, the diameter of the tube, T t , must be kept large (about 2 cm) 
to provide adequate pumping speed to the charcoal. In this apparatus the pumping 
speed was limited by the 7.0-cm length of 20-mm OD tubing which connected the 
charcoal to the larger stopcock. The pumping speed of this tube was 10 liters/sec 
(as read from Fig. 2-1, p. 92, of Ref. 21). When the large stopcock was open, the 
carbon was connected to a U-tube cold trap, T,,, 20 mm OD and 15 cm long, which 
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Fig. 5-1 Volumetric Apparatus for Determining Rate of Vacuum Desorption 
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Could be immersed in liquid nitrogen. A Consolidated Vacuum Corporation Pirani gage 
Vas connected on a side seal, and a second stopcock isolated the system from the 
V-acuum pumps and from a McLeod gage used to calibrate the Pirani gage. Obvious 
1 mprovements would be to eliminate stopcock grease by using a metal solenoid valve 
pr bakeable gate valve in place of S , and to use a disphragm-type pressure transducer 

] n place of the Pirani gage . 

Carbon granules were introduced into by blowing a hole in the top of the tube, 

Which was later sealed. The optimum weight of carbon used was a small problem. 

-foo large a weight would cause results to be affected by pumping speed, and too small 
* weight would cause sampling errors and errors due to contaminants desorbed from 
ihe apparatus walls. In these experiments 20-mg samples were used, although 10-mg 
Samples were also tried in initial experiments. 

A typical experiment was carried out as follows: The charcoal was first degassed in 
kigh vacuum for 4 hr at 300- 400°C by slipping an electric furnace over the tube, . 
During bakeout the large stopcock, S , was turned at frequent intervals to degas the 
stopcock grease. (Bakeout temperature was varied from 250°C to 400°C with no signifi- 
cant change in desorption rate of butane at room temperature.) The carbon was next 
cooled to room temperature and the stopcock, S , was closed. A desired pressure of 
butane was then admitted to the volume containing the Pirani gage and large stopcock. 
This volume was 366 cm 3 for the case described here. Thus, from the perfect gas 
law, 0.364 torr in 366 cm 3 gave an initial loading of 8 ml STP/g on a 20 mg sample, 
corresponding to the initial loading in the microbalance experiments. This pressure 
of butane was admitted to the carbon at 26°C, and the rate of adsorption was determined 
by measuring the amount adsorbed as a function of time. This rate was initially too 
fast to determine, due to time lags in the Pirani-gage system. Within about 10 sec 
the pressure dropped by a factor of 3 to 10, and adsorption equilibrium was established 
in 10-20 min. (A point on the adsorption isotherm was thus measured in each experi- 
ment. ) The stopcock, S , was then closed and the cold trap, T 1 , was completely 
immersed in liquid nitrogen. The desorption rate was then measured by opening the 
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Fig. 5 8 Gravimetric Desorption of n-Butane From BD Activated Carbon: Variation 
with Particle Size 
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Fig. 5-9 Gravimetric Desorption of n-Butane From BD Activated Carbon at 87 C (360 K) 


Data from experiments on 1.0-mm charcoal at 87°C using both the shorter and the 
longer furnace designs are compared in Fig. 5-10. Initial desorption rates were 
somewhat higher using the longer furnaces but approached the same rate after about 
30 percent of the n-butane had been desorbed. Data are shown from the two sets of 
longer furnaces, in which only the heater winding material was different. While the 
furnaces with magnetic windings resulted in very unstable weight readings, desorption 
rate data do not appear to have been appreciably affected. 

Data for an additional set of experiments in which n-butane initial loading was varied 
are plotted semilogarithmically in Fig. 5-11 and 5-12. In the former plot the actual 
amount adsorbed is shown as a function of desorption time, while in the latter the 
amounts adsorbed are normalized by dividing out the initial amount adsorbed. 


Data for the desorption of toluene from 1.0-mm particles at various temperatures are 
shown logarithmically in Fig. 5-13 and 5-14. The initial amount adsorbed was 50 ml 
STP/g. In addition, the results of some preliminary experiments with toluene on both 
BD charcoal and sorbent LS2A are shown in Fig. 5-15, which unlike previous plots is 
on linear coordinates. 

Propane desorption data for two particle sizes at 27°C are shown semilogarithmically 
in Fig. 5-16. The initial amount adsorbed was 1.1 ml STP/g. 

5.2.2 Multiple -Component Desorption 

The measurement of the desorption rate of two components from charcoal involves a 
new factor, since the rate of total weight change as read from the microbalance record 
must be supplemented by analysis of the off-gas, so that the individual rates may be 
computed. In the present work a quadrupole mass spectrometer was used for the 
instantaneous gas analyses. 


86 



FRACTION REMAINING ADSORBED 


1 . 0 * 


Hangdown Tube Furnace: 


0.5 


0.2b 


ID Length 
(cm) (cm) 


o 

7.8 

12.7 

□ 

7.8 

12.7 

o 

7.8 

12.7 

A 

4.5 

30.5 

V 

4.5 

30.5 


4.4 

30.5 


4.4 

30.5 


Winding 

Mat'l 

Ni/Cr 

Ni/Cr 

Ni/Cr 

Fe/Cr 

Fe/Cr 

Ni/Cr 

Ni/Cr 


Preadjustment for Thermo- 
gravimetric Effects 


Run 

No. 


(Done Before Run No. M13-1) 
(Done Before Run No. M13-1) 
Height Adj. w/Ar at 20 /xHg 
C-tube Volt. Adj. at 20 pHg 
C-tube Volt. Adj. at 20 ^Hg 
None, Checked w/Ar at 20 pHg 
None, Checked w/Ar at 20 p Hg 


M18-1 

M18-2 

M18-4 

M18-16 

M18-18 

M18-19 

M18-20 


Amount Adsorbed in All Cases — 7.7 ml STP/g 



0.05 b 


0 L 
0 


-4— 1 U 

100 200 

DESORPTION TIME (sec) 


300 


Fig. 5-10 Gravimetric Desorption of n-Butane From 16 x 20 BD Activated Carbon: 
Effect of Various Heating Arrangements at q 7°C 
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Fig. 5-11 Desorption of n-Butane From 16 x 20 BD Activated Carbon at 87°C and 
Various Amounts Previously Adsorbed 
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Fig. 5 13 Gravimetric Desorption of Tolune From 16 x 20 BD Activated Carbon 
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It follows that an equation is needed which allows the individual desorption rates to be 
computed from the total weight change and the instantaneous gas analyses. This 
equation is 


where 


N, = 


w 


t M t + ayx/x^ 


(5.1) 


N f = 


w = 


M f % - 

x f »b - 


instantaneous desorption rate of toluene (mol/sec) 
total desorption rate (g/sec) 

molecular weight of toluene and butane, respectively (g/mol) 
mole fraction of toluene and butane, respectively, in the off-gas 
at that time 


Equation 5. i follows directly from the equation of mass balance: 


M t N t + “iA “ * 


(5.2) 


plus the equation 




(5.3) 


which expresses the basic fact that the mole ratio in the off-gas at any instant is equal 
to the relative rates of desorption. An equation exactly the same as Eq. 5.3 with 
subscripts interchanged holds for butane. 


Thus, the total weight of toluene desorbed up to any time can be computed by integrating 
Eq. 5.1 graphically: 

t 

N t “ J M t + M fe (X t /X b ) dt ( 5 * 4 ) 

where = number of moles of toluene desorbed up to time t 
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[n determining the off-gas composition, the mass spectrometer peak at M/e 51 was 
chosen for toluene and at M/e = 43 for n-butane. If the peak heights at M/e = 51 and 
13 are, respectively, h gl and h^, then the mole ratio is computed from 


5t = S 43 h 51 (5.5) 

*b " S 51 h 43 


where S and S are empirical sensitivity factors, related to ionization probabilities. 

The sensitivity ratio in Eq. 5.5 can be determined either empirically from a known gas 
mixture or alternatively from Eq. 5.4 since only the correct ratio will give the correct 
number of moles of toluene desorbed at infinite time. 

For this experiment at 27°C, toluene was first adsorbed at 8 x 10 torr equilibrium 
pressure which gave an adsorption of 31. 6 ml STP/g. Butane was then admitted at 
0.30 torr pressure which gave an equilibrium adsorption of 1. 125 ml STP/g. This 
relatively high butane pressure was necessary to get a significant adsorption. 

Evacuation was then started, with recording of both the total weight versus time, and 
the mass spectrometer peaks versus time. 


The total weight versus time is shown in Fig. 5-17, while Fig. 5 18 displays the mole 
ratio versus time, using a sensitivity ratio S^/S^ = 0.04. This ratio was determined 
from the fact that it gave a good mass balance when Eq. 5.4 was, in effect, integrated 
out to infinite time. The increase in this ratio during later stages of the desorption is 
probably due to n-butane being gradually reevaporated from the liquid nitrogen trap 
used in pumping the system. The ratio dropped to a much lower value when the trap 
was filled at about 1550 sec desorption time. 


The rate of toluene desorption versus time is shown in Fig. 5 19 and was computed 
from Eq. 5.1 using the sensitivity ratio of 0.04. 
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Initial Weight, w° = 50.618 mg 

Weight With Adsorbed Toluene, wt = 57.206 mg 
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Fig. 5-17 Desorption of n-Butane and Toluene From 8 x 12 BD Activated Carbon at 27°C 




Fig. 5-18 Mole Ratio of n-Butane to Toluene During Vacuum Desorption From 8 x 12 BD Activated Carbon at 27°C 
as Determined From Mass Spectrometric Analysis of Gases Being Desorbed 


DESORPTION RATE (fi mol/g-sec) 




The integrated version of Fig. 5-17, converted to the fraction remaining adsorbed 
versus time, is shown in Fig. 5-20. Also shown in Fig. 5-20 is a curve representing 
previous desorption data for toluene on 0.1-cm particles, one-half the size of the 
present ones. This curve is from a previous experiment in which somewhat more 
toluene (50 ml STP/g) was initially adsorbed than in the present experiment (31.6 ml 
STP/g). This curve has been replotted for comparison purposes with zero time 
assigned to the point where 31.6 ml STP/g remained adsorbed. It is predicted that 
desorption time for this smaller particle size would be decreased by a factor of about 
1.8, as was experienced with n-butane desorptions on various particle sizes as reported 
in section 5.2.1. The ratio here is somewhat less - about 1.2; however , it would be 
best to compare the present data with pure toluene desorption data obtained on the same 
particles at the very same initial loading, since there is some variation in the charcoal 
and probably as increased cooling effect in the case of the data from the previous work 
due to the larger amount initially adsorbed. It thus appears that the desorption rate of 
toluene is not greatly affected by the presence of n-butane. 

5.2.3 Volumetric Apparatus 

Experimental data were obtained using the volumetric system for desorption from two 
particle sizes, 16 x 20 mesh and 8 x 12 mesh (corresponding to average diameters of 
0.1 and 0.2 cm, respectively). Experiments were performed at 26°C on two carbon 
sample weights, 10 and 20 mg. Within the experimental error of reading the Piram 
gage (2 to 5 percent), all desorption rates were essentially second order; that is, a 
plot of the inverse of the fraction remaining versus time was a straight line over the 

whole range of time. 

Data for the experiments with 10-mg samples were not very reproducible, especially 
in the case of the 0.2-cm particles, 10 mg of which consists of only 2 or 3 particles. 
In one case a sample of 10 mg of 0. 2-cm particles actually desorbed faster than a 
similar sample of 0.1-cm particles, resulting in an experimental scatter of more 

than a factor of 2. 
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Desorption of Toluene From 
Desorption of n-Butane 


8 x 12 BD Activated Carbon at 27°C During Concurrent 



Data for the experiments using 20-mg samples were considerably more reproducible 
than those using 10 mg. The 20-mg data are given in Fig. 5-21, in which time for half- 
desorption is plotted against initial loading of n-butane. Straight lines of unit slope 
have been drawn through the data, since second-order desorption kinetics seem to be 
well established. Both the upper line for 0.2-cm particles and the lower line for 
0.1-cm particles are within 10— 15 percent of the values obtained in gravimetric 
experiments. Experiments were reproducible to 5 percent or better for a given set 
of charcoal granules at a given loading, but changing to a new sample or to a new 
loading of n-butane introduced a scatter of up to 20 percent. 

Some additional details of the volumetric data are presented in Appendix K. 
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Section 6 

VACUUM DESORPTION FROM SINGLE PARTICLES — DISCUSSION 


Once extensive vacuum desorption data have been obtained for single activated carbon 
particles, it should be possible to gain considerable insight into the physical processes 
that are taking place when such particles are subjected individually to the approximate 
equivalent of space vacuum. In this section, data from Section 5 are discussed m 
relation to various physical models for desorption. 

The kinetics of desorption (desorption rate, rate order, rate constants, and desorption 
activation energies) are discussed with respect to variations in temperature and particle 
size for a few selected contaminants in order to have some basis for predicting the 
effect of such variations on the performance of packed beds during vacuum regeneration. 
Various mechanisms for single-particle desorption are examined for similar reasons. 
Cooling effects are treated in order to determine the magnitude of the cooling which 
naturally results when single particles with no thermally conductive path to the 
surroundings are subjected to fast vacuum desorption. This is needed in order to 
make valid comparisons with vacuum desorption data from packed beds, in which 
thermal conductivity of closely packed particles is much better than in the 
single -particle case. 


6. 1 FIT TO APPARENT RATE ORDER 

A means of correlating the gravimetric kinetic data is simply to seek the empirical 
order of the desorption rate. That is, it is assumed that the desorption rate follows 

simple power-law kinetics: 
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where 


(ml STP/g-sec) 

[(ml STP/g)n-l/ sec ] 
(ml STP/g) 


I'd = desorption rate 
kj = desorption rate constant 
q r = amount remaining adsorbed 
n = desorption rate order 

Desorption rates which obey such kinetics give rise to a Freundlich isotherm since 
the adsorption rate may be defined by 


where 

(ml STP/g-sec) 
(torr-ml STP/g-sec) 
(torr) 

The net rate of adsorption is 


r a = adsorption rate 
k a = adsorption rate constant 
p = pressure 


( 6 . 2 ) 


r 

n 



i 

k,q 

dr 


At equilibrium, this net rate is equal to zero, so that 


% = (kAj 1 /n P 1/n = k p 1/n 


(6.3) 


(6.4) 


where 

kp - Freundlich constant [ (tor r) n ml STP/g J 

Furthermore, a linear potential plot (see section 4.1) is consistent with Eq. 6.3, and 
hence it is expected that the simple power-law kinetics of Eq. 6.1 are appropriate. 
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The most direct way to determine desorption rate order is by the method outlined in 
Appendix L. There it is shown that a plot of log q versus log r (t being a characteristic 
time) is linear. The slope at any point during the desorption is equal to 1 - n, where 
n equals the rate order. Also, the rate constant may be calculated from the intercepts 
of such a plot, as shown in Appendix L. 

Using this method, rate orders and constants were graphically determined for the 
n-butane desorption results given in section 5.2.1. Rate orders were found to be 
essentially constant after about half the n-butane had been desorbed. For desorption 
from 0.5-mm particles, the slopes of the lines thus determined are all only slightly 
greater than -1, as shown in Fig. 6-1. Rate orders in this case are then equal to 
slightly more than 2. For 0.20-cm particles, the rate order increases from 2.3 at 
26°C to n = 2.8 at 87°C, as shown in Fig. 6-2. Comparisons of rate orders for different 
particle diameters at 27 and 87°C are given in Fig. 6-3 and 6-4, respectively. Rate 
orders for 0.025- and 0.05-cm particles are near 2.0, while those for 0.10- and 
0.20-cm particles are somewhat higher. Rate orders for desorptions from 0.10-cm 
particles at 87°C with initial loadings other than the usual 7.7 ml STP/g are shown in 
Fig. 6-5, and are approximately equal to 2.0 except for very small initial loadings 
where the value goes as high as 2.5. 

In Fig. 6-6, the data from a second set of runs on 0.05-cm particles are treated in the 
same way as in Fig. 6-1. The temperature range is a bit greater for these data but 
still the rate order, n , varies between 2. 1 and 2.4. If the data of both these runs are 
taken together, there is no recognizible correlation of the order, n , with temperature. 
Furthermore, the data depicted in Fig. 6-6 are analyzed for n by graphical differ- 
entiation (to be described more completely below) in Fig. 6-7. Here the orders vary 
from 2.3 to 2.5. From all of this, it may be concluded that the order is known to 
be between 2.0 and 2.5 or so far the small particles with no easily recognizible trend 
with temperature. For the longest particles (0.20-cm), there appears to be a slight 
but significant increase in the order, perhaps up to as high as 2.8. Rate constants for 
each of these plots determined from the intercepts (see Appendix L) are indicated on 
the respective rate order curves. 
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Fig. 6-1 Determination of Rate Order and Constants for Gravimetric Desorption 
of n-Butane From 30 x 40 BD Activated Carbon (d =0.05 cm) 
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Temp. 

Run 

Order 

Constant 

<°C) 

No. 

n 

[(ml STP/g) 1 " r 

26 

M15-2 

2.3 

4.9 x 10" 4 

47 

M15-5 

2.4 

8.8 x 1(T 4 

67 

M15-4 

2.6 

1.3 x 10 -3 

87 

M15-3 

2.8 

1.9 x IQ" 3 


Initial Amount Adsorbed = 7. 7 ml STP/g 


AMOUNT REMAINING ADSORBED (ml STP/g) 

Fig. 6-2 Determination of Rate Order and Constants for Gravimetric Desorption 
of n-Butane From 8 x 12 BD Activated Carbon (d = 0.20 cm) 
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Fig. 6 3 Determination of Rate Order and Constants for Gravimetric Desorption of 
n-Butane From BD Activated Carbon at 27°C: Variation With Particle Size 
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Fig. 6-5 Determination of Rate Order and Constants for Desorption of n-Butane 
From 16 x 20 BD Activated Carbon: Variation With Amount Previously 
Adsorbed 
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Fig. 




6-6 Determination of Rate Order and Constants for Gravimetric Desorption of 
n-Butane From 30 x 40 BD Activated Carbon (d p = 0.05 cm) 
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n alternate method of determining rate order and constants, the one more commonly 
sed, requires graphic differentiation of the desorption curves at various points to 
etermine instantaneous desorption rates. These rates are then plotted against the 
mount remaining adsorbed on log-log coordinates, as shown in Fig. 6-8 for the same 
xperiments evaluated in Fig. 6-5 for various amounts initially adsorbed on 0.1-cm 
arbon. Since the axis for the amount remaining adsorbed in Fig. 6-8 (and also in 
'ig. 6-7) is negative, the rate order is the negative of the slope of the lines. This 
lethod is less precise because the slopes obtained are greater by one than those 
btained in the method of Appendix L. Thus, the deviation of rate order from n — 2 is 
ess obvious in Fig. 6-7 and 6-8 than in plots similar to Fig. 6-4. 


late constants are somewhat easier to obtain from this type of plot, the value being 
letermined directly from the intercept with the q = 1 axis, where q r = amount 


■emaining adsorbed. Units for the rate constants are (ml STP/g) 1 "sec" 1 . 

Dnce desorption rate order has been established, it is possible to determine rate 
:onstants by yet another method. Since the rate expression is 


- dq/dt = k d q” 


(6.5) 


it may be integrated to 


I 

Thus, a plot of q 1_n versus t will have a slope of kj . Such plots for four of the 
experiments whose rate order was determined in Fig. 6-1 are shown in Fig. 6-9, in 
which the resulting values of the rate constant are indicated. The rate order was 
n = 2 in all cases. A similar plot for rate orders determined in Fig. 6-2 is given in 
Fig. 6-10. In this case, rate orders varied from 2.3 to 2.8; thus, it was necessary 
to calculate q* " n for each point and then to divide the slope by n - 1 in order to get 
each rate constant, which is shown next to its respective curve in Fig. 6-10. 
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Fig. 6-9 Second-Order Rate Plot for Desorption of n-Butane From BD Activated Carbon 











DESORPTION TIME (sec) 

Fig. 6-10 Determination of Rate Constants for Desorption of n-Butane From 0. 2 cm BD Activated Carbon Particles 



Toluene desorption rate orders were determined for the one set of experiments with 
1-mm BD charcoal, in which temperature was varied from 27 to 107°C. Those 
determined by the method of Appendix L are shown in Fig. 6-11, which is plotted on 
log-log coordinates with unequal axes in order to bring the slopes of the lines nearer 
to unity for accuracy. In Fig. 6-12, desorption rates determined from the weight 
records have been plotted versus the amount remaining adsorbed on log-log coordinates 
with the latter axis reversed as in Fig. 6-8. In both cases, the rate order varies 
from 3.5 at 27°C down to 2.6 at 107°C. 

The resultant rate orders and rate constants along with the reciprocal slopes of the 
Freundlich-type isotherms given in section 2.2 are tabulated in Table 6-1 for n-butane 
and Table 6-2 for toluene. 

It has already been shown that the apparent rate order is related to the excess free 
energy of adsorption, (p : 


n 


V <P 
m 

kT 


(4. 24) 


The quantity <p may be obtained from the slope of the potential plot (see Fig. 4-3). 
For toluene, <p is constant over a large range of amounts adsorbed and independent 
of temperature. Thus, Eq. 4.24 gives explicitly the magnitude and temperature 
dependence of the reaction rate order. In Fig. 6-13 the observed rate orders are 
plotted as a function of reciprocal temperature and compared to Eq. 4. 24 with a 
value of <p = 19.6cal/cm calculated from the toluene potential plot. The points 
are very close to the line, indicating that Eq. 4. 24 predicts quantitatively both the 
value of n and its temperature dependence. 

Even though there is no easily recognized temperature dependence in the orders 
for butane desorption, it is instructive to see how those values compare with the 
prediction of Eq. 4.24 for n-butane. In this case, there is a slight curvature in 
the potential plot over the range of interest, but <p has been evaluated for the 
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Fig. 6-11 Determination of Rate Order a 
From 16 x 20 BD Activated Ca 


U 27°C (300°K) Run No. M18-25D 
O 67°C (340°K) Run No. M18-48D 
O 107°C (380°K) Run No. M18-60D 

Amoun t Adsorb ed = 50 ml STP/g 







Table 6-1 

SUMMARY OF n-BUTANE DESORPTION RATE EXPERIMENTS 

(q a = 7. 7 ml STP/g) 


Rate Order 
(Values of n) 

- 

Nominal Temperature — ► 

27°C 

47°C 

67°C 

87°C 

107°C 

50 x 70 Mesh, d p = 0.025 cm 

2.0 

, „ 




30 x 40 Mesh, dp = 0.05 cm 

12.0 

2.2 

2.3 

2.2 



\2.5 

— 

2.3 


2.3 

16 x 20 Mesh, d p = 0.10 cm 

2.3 

- 

- 

2.0 

_ 

8 x 12 Mesh, d p = 0.20 cm 

2.3 

2.4 

2.6 

2.8 

- 

Rate Constants 






l(ml STP/g) 1_n /sec] 






Nominal Temperature — ► 


47°C 

67°C 

87°C 

107°C 

50 x 70 Mesh, d p = 0.025 cm 

3.8 x 10 -3 






30 x 40 Mesh, d p = 0.050 cm 

11.7 x 10-3 

3.5 x 10“ 3 

6.6 x 10-3 

1.2 x 10~ 2 



19.0 x 10" 4 

— 

4.6 x 10”3 

— 

1.8 x 10 -2 

16 x 20 Mesh, d p = 0. 10 cm 

8.5 x io~ 4 ( a ) 

— 

— 

5.4 x 10" 3 


8 x 12 Mesh, d p = 0.20 cm 

4.9 x 10 _4(b) 

8.8 x 10 -4 

1.3 x 10 -3 

1.9,x 10“ 3 

- 

Reciprocal Slopes of Freundlich 

-Type Isotherms 




at q a - 3.85 ml STP/g (Half-Desorption) 





Nominal Temperature — 

27°C 

47°C 

67°C 

87°C 

107°C 

50 x 70 Mesh, d p = 0.25 mm 

— 

— 





30 x 40 Mesh, d p = 0.50 mm 

|l . 84 

2.27 



1.59 


11 . 93 





16 x 20 Mesh, dp = 1.0 mm 

1.93 




1.59 

8 x 12 Mesh, dp » 2.0 mm 
fa\ A 1 v in-4 — 

1.79 


1.76 




(b) 7.0x10 4 = best fit for n = 2. 
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Table 6-2 

DATA FOR TOLUENE SORPTION ON 16 x 20 MESH (0.10 cm) BD CHARCOAL 


Temperature 

T 

(°C) 

Freundlich 
Isotherm Data 
at 25 ml STP/g* ! 

Overall 
Desorption 
Rate Order 
n 

Overall 
Desorption 
Rate Constant 
kd 

[(ml STP/g) 1-n sec _1 ] 

Half 

Desorption 

Time 

*1/2 

(sec) 

Toluene 

Vapor 

Pressure 

P° 

(torr) 

Adsorption Pressure 
for 50 ml STP/g 

Absolute 

P 

(torr) 

Relative 

p/P° 

Slope 

1 /Slope. 

27 

0.27 

3.7 

3.5 

1.42 x 10-7 

630 

32 

0.0060 

1.87 x 10 -4 

67 

0.27 

3.7 

3.0 

5.5 x 10“ 6 

95 

177 

0.091 

5.1 x 10" 4 

107 

0.27 

3.7 

2.6 

7.8 x 10“ 5 

23 

660 

1.070 

1.62 x 10 -3 


♦From Fig. 2-10. 




loading equal to one-half the amount initially adsorbed. Also the large particle size 
(0. 20-cm)data have been omitted, since as mentioned earlier there seems to be 
an increase in the order for these particles relative to the smaller ones. The com- 
parison between these data and values calculated from Eq. 4. 24 is also made in Fig. 
6-13. Within the experimental scatter of these data, which admittedly is quite large, 
Eq. 4-24 is again quite adequate using 0 = 14.6cal/cm3, which was derived as 
indicated from the appropriate point on the n-butane potential plot. 


These quantitative comparisons between rate parameters and equilibrium quantities 
are confiriming evidence of the generality of the potential-plot correlations and the 
validity of using them to predict some features of the rate phenomena. 

Because of the importance of the effect of particle size on the desorption rate in 
understanding the details of the desorption mechanism, one series of runs using 
n-butane was performed specifically to evaluate this dependence. The data for these 
runs are shown in Fig. 5-6 and 5-7. From these figures, the half-times as a function 
of particle size can be determined. These are shown in Fig. 6-14. The data included 
are quite linear with respect to particle size throughout the whole range tested. 

In section 6.3, it is demonstrated that this implies that the rate-limiting step is 
desorption from the external surface of the particles, accompanied by rapid surface 
diffusion within the particles. 

6.2 CALCULATION OF APPARENT ACTIVATION ENERGIES 

The most direct way to calculate isosteric activation energies for desorption is to 
plot the logarithm of absolute desorption rate at a given coverage versus reciprocal 
temperature. The activation energy can then be determined by multiplying the slope 

of such a plot by the molar gas constant (and by In 10 to convert the slope to natural 
logarithms), as shown by a version of the Arrhenius equation 

E d=™[^| (6.7, 
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3.5 


2.5 



DESORPTION RATE ORDER, n 
Fig. 6-13 Temperature Dependence of Desorption Rate Order 


DESORPTION HALF-TIME, t A _ (sec 



where 


Ej = activation energy of adsorption (cal/mol) 

R = gas constant = 1.987 cal/moT’K 
T(j = desorption rate (ml STP/g-sec) 
t = time (sec) 

T = absolute temperature (°K) 
q = amount adsorbed (ml STP/g) 

This has been done for n-butane on 0.5-mm carbon in Figs. 6-15 and 6-16 for toluene 
on 1.0-mm carbon in Fig. 6-16. Desorption rate data were taken from Figs. 6-8 
and 6rl2 , respectively. 

The apparent activation energies calculated from Figs. 6-15 and 6-16 are tabulated 
in those figures. The values are only approximate since data for the three tempera- 
tures used, 27, 67, and 107°C do not fall exactly in a straight line. These values 
are somewhat less than the heats of adsorption reported in section 4.2.1 and they 
increase slightly with decreasing coverage. Adiabatic cooling effects could be 
responsible for the low apparent activation energies, since at the higher temperatures 
the individual granules are actually cooled below the nominal furnace temperature in 
order to supply the energy necessary to desorb the n-butane. Otherwise, it would be 
expected that the activation energies would always be slightly greater than the 
isosteric heats, approaching them quantitatively only when the activation energy for 
adsorption is zero. 

6.3 THEORY AND MECHANISMS 

In order to interpret the experimental results of the preceding sections, the various 
possible mechanisms of desorption must be examined. In the desorption of a contaminant 
from a charcoal granule, there are four possible rate-determining processes: 

(1) Pumping speed 

(2) Rate of pore diffusion from within the charcoal granule to the external 
surface 
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DESORPTION RATE (ml STP/g-sec) 
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(3) Rapid gas -phase diffusion within the pores, and slower desorption of a 
molecule from the pore wall into the pore void space 

(4) Rate of evaporation from the external granule surface when this is supplied 
by a very fast surface diffusion along the pore walls 

Fortunately , all of these mechanisms but the last one may be eliminated, as described 
in the following paragraphs. 


6.3.1 Pumping Speed Rate-Determining 

When pumping speed is controlling, the concentration at the charcoal granule will be 
related to the equilibrium (isotherm) pressure. Since the total rate of removal equals 

the pumping speed times the pressure, the rate per unit volume of bed may be expressed 
as 


= S c (S) 
W p e ' 


- adsorbent concentration per unit bed volume (mol/cm 3 ) 

time , 

(sec) 

particle density (g/cm 3 ) 

pumping speed (cm3/ S ec) 

gas phase concentration in equilibrium with the adsorbed (mol/cm 3 ) 
concentration S 

weight of carbon used ^ 


= time 


c e (S) 


The above equation shows that the half-time depends on the weight of charcoal used 
(i.e. , on the number of granules used in this particular experiment). Rate will be 
independent of particle size (for a given weight used) and will depend strongly on tem- 
perature, since the equilibrium pressure depends strongly on temperature. The 
apparent activation energy will be exactly the heat of adsorption. The apparent 
"order of reaction" will be the order of the Freundlich isotherm. (That is, if the 
adsorption obeys the Freundlich isotherm with n = 2, desorption will be second order.) 
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Since in the present work desorption rate depends on particle size (at constant charcoa 
height), it appears certain that in the experiments on the microbalance the pumping 
speed is not rate-determining, and this is cheeked by experiments in the volumetric 
system (see section 5.2.3) in which the pumping speed was varied. This possi e 

mechanism can thus be eliminated. 

£,.3.2 Pore Diffusion Rate-Determining 

When diffusion through the pores to the particle surface is controlling, then the time 
^ half desorption can be shown (see Appendix M) to be given by 


l l /2 


d 

= k-^ — 

^ n c 


(6.9) 


Where 

t|y 2 = time for half-desorption 
d = granule diameter (presumed spherical) 

f Ir 

D = diffusion coefficient 
S Q = initial adsorbed phase concentration 
c o = initial equilibrium gas phase concentration 
K = constant which depends on the shape of the 
adsorption isotherm and which increases 
as the isotherm becomes more nonlinear 
(e.g. , K = 0.008 for the linear isotherm) 


(sec) 

(cm) 

(cm 2 /sec) 
(mol/cm 3 ) 
(mol/cm 3 ) 


It follows from Eq. 6.9 that for this mechanism the desorption rate is inversely 
proportional to the square of particle size. The apparent activation energy is exactly 
equal to the heat of adsorption for the case of gas diffusion through its pores, but for 
two-dimensional surface migration through the pores the exact activation energy 
beomces somewhat hazy due to uncertainties in formulating the diffusion equation, 
although it is certain to be less than the heat of adsorption. 
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This possible mechanism may be eliminated, since the present work shows the time 
for half-desorption to be proportional to the first power, not the square, of granule 
size. 

A rate -controlling mechanism of slow surface diffusion in the adsorbed layer to the 
external particle surface would result in the same dependence on R as above. The 
expression for t^^ would also be the same, except that D would refer to surface 
diffusion rather than pore diffusion. The apparent activation energy would be less 
than AH . Clearly, this mechanism is ruled out by the same argument as for pore 
diffusion. 

6.3.3 Desorption From Pore Wall Rate-Determining 

When the actual surface desorption step from the pore wall into the pore voids is 
rate-determining, the desorption rate will be independent of granule size, since the 
charcoal will behave as if its entire surface area were exposed as a plane surface. 
(That is, by hypothesis, pore diffusion is much faster than desorption.) This 
mechanism can thus be eliminated since, in the present results, rate depends on 
granule size. 

6.3.4 Desorption From External Granule Surface Rate-Determining 

If two-dimensional surface diffusion along the pore wall is very rapid, then a possible 
rate-determining step is the rate of desorption from the external granule surface, into 
the gas phase, supported by the very rapid surface diffusion. In Appendix N, it is 
shown that the rate for this process may be formulated as 

- jja = (1.98xl0‘ 3 )^^ap e (6.10) 

where 

q = adsorbed concentration (ml STP/g) 

T = temperature (°K) 
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V 

d 

P 

P P 

a 


p, 

The importance of Eq. 6. 10 is that it predicts the correct dependence of rate on granule 
size, and apparently this is the only mechanism that gives the inverse first power of 
particle size. The only unknown in Eq. 6.10 is a , the sticking coefficient, so that 
a may be calculated from the experimental rate data. Thus, a severe test of this 
mechanism is that the absolute value of a must be reasonable (i.e. , less than unity 
but not too small. Sticking coefficients for vapor molecules onto liquid surfaces are 
known to be in the range 0.02 to 1.0. (See, for example, Wyllie, Ref. 22, p. 389.) 


average molecular velocity 
particle diameter 
bulk density 

sticking coefficient for adsorption (i.e., the probability 
that a molecule striking the charcoal surface will be 
adsorbed) 

pressure of adsorbate in equilibrium with the adsorbed 
concentration q 


(cm/sec) 

(cm) 

3 

(g/cm ) 


(torr) 


To apply Eq. 6.10 conveniently to the experimental data, if the adsorbate obeys 
approximately the Freundlich isotherms 


then the equation may be rewritten 




( 6 . 12 ) 


where 


k 


1.98 x 10 


-3 273 v 


T d p 
P P 


(6.13) 


The time for 50% desorption may then be formulated from standard kinetics as 

500 (2 n_1 - 1) d p P p q o T f6. 14i 

V2 = (n-1) v ap° 273 

e 
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where 


P e pressure of adsorbate in equilibrium with the initial amount adsorbed 
% < Le *> amoun t adsorbed at zero time) (torr) 

From Eq. 6. 14, the sticking coefficient a may be calculated for butane at 27°C, 

using the experimental values for L /n . Since n = 2n-nsfi<r/ 3 - „ 4 

1/2 uoe n - p - 0.86g/cm , v=3,3xl0 

cm/sec, and for the data of Fig. 5-6, ^ = 7.7 ml STP/g, p e = 9. 6 X 1(T 3 torr, 
d p -0.1 cm, t 1/2 =95 sec, a may be computed at 0.0115. Since this is a very 
reasonable value, this may be taken as confirmatory evidence that the above mechanism 
is a reasonable one. A similar calculation for toluene desorption at room temperature 
gives a - 0,053, which is also very reasonable. 

If similar calculations are made for a at more elevated temperatures, the values turn 
out lower. For example, at 107’C for toluene the computed value is a factor of almost 
10 lower (0.0060). While it is true that energy accommodation is more difficult at 
higher gas temperatures, it is believed that this decrease in a at the higher tempera- 
tures is associated not with energy accommodation but with uncertainties in the 
temperature brought about by the adiabatic cooling which accompanies the rapid 
desorption rates observed at higher temperatures. 

Further evidence for this particular rate-controlling step is found from the analysis 
of the entropies of adsorption (see Section 4). The experimentally determined entropies 
are consistent with a freely mobile adsorbed layer. Such mobility, of course, is 
necessary for the rapid surface diffusion which the kinetics of desorption suggests. 

A summary of factors relating to these various mechanisms for vacuum desorption 
from single particles is presented in Table 6-3. 

6.4 COOLING EFFECTS 

All the possible mechanisms discussed in the last section predict an apparent activation 
energy equal to the heat of adsorption, * essentially because for each mechanism the 


Actually, the last mechanism discussed, desorption from the 
particle, predicts an activation energy somewhat higher than 
likely a function of temperature. 


external surface of the 
q st , since a is most 
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Table 6-3 

KINETICS OF VACUUM DESORPTION FROM SINGLE GRANULES 


Possible rate-determining 
mechanism ( a ) 

Dependence of rate on 
particle size, d^ 

Time for 50% 
desorption, t^y 2 

Pumping speed 

Independent, but depends 

i s 


on weight of particles 

W _J_ o 
n Sc 


used 

P P O 
d 2 S 

Gas diffusion in pores 

1/d 2 

P 

K» — 

D K C o 

Evaporation from internal 



pore wall with very fast 
gas diffusion in pores 

Independent of particle 
size 

1/k d 

Evaporation from granule 


. d S 

surface fed by fast surface 
diffusion 

1/d 

P 

1 p o 

6k. e 
A o 


(a) For these mechanisms , apparent activation energy = AH a = enthalpy of 
adsorption. 

(b) K* is a numerical constant depending on the shape of the adsorption 
isotherm. 

rate is proportional to the equilibrium pressure, at a given loading of adsorbate. 
Since at constant loading of adsorbate, the change in equilibrium pressure with 
temperature is given by 


d In p _ _ %t 
d(l/T) “ R 


(6.15) 


it follows that the desorption rate will also have q g ^ as its apparent activation energy. 

Here, q is the isosteric heat of adsorption. Physically, this corresponds to the 
st 
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fact that the adsorbed molecule must acquire its heat of adsorption before it can 
desorb, so that the activation energy should be at least this large. 

Experimentally, however, the apparent activation energy is much less than the heat of 
adsorption, e.g. , for n-butane the initial apparent activation energy is 10 kcal/mol or 
less, while the heat of adsorption is 12 — 14 kcal/mol. (See sections 4.2 and 6.2.) 


This dilemma may be resolved if one assumes that at elevated temperatures the 
granules are not at the furnace temperature but at some lower temperature. Using 
the furnace temperature to compute activation energies in such cases would result 
in apparent activation energies somewhat lower than the correct values. 

The charcoal granules can be colder than the furnace temperature for either of two 
reasons. (1) poor furnace design, so that the granules are not completely surrounded 
by the furnace, and (2) adiabatic cooling during desorption. The former can be 
eliminated from consideration, since the same furnaces are used to measure the 
adsorption isotherms, so that too low an isosteric heat of adsorption would also be 
measured. Also two different sets of furnaces have been used. The most likely 
explanation is that there exists an adiabatic cooling of the granules due to the loss 
of the heat of desorption as desorption proceeds. If the charcoal granule were 
completely thermally isolated from its surroundings, the cooling effect due to 
evaporation of adsorbate would be, from a simple heat balance, 


where 


(AT) 


<V q o 


ad 22.4 xC 


(6.16) 


< AT >ad 

qst 

C s 

% 


adiabatic temperature drop due to evaporation of adsorbate (°C) 
heat of adsorption (kcal/mol) 

heat capacity of the charcoal (cal/g°C) 

initial amount of adsorbate adsorbed (ml STP/g) 
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For butane with q Q - 7. 7 ml STP/g, C g = 0. 17 cal/g°C, q gt = 12 kcal/mol, the 
adiabatic cooling effect is 24°C. For toluene with the larger loading, q Q - 50 ml STP/g 
and the larger heat of adsorption (about 16.5 kcal/mol), the adiabatic temperature drop 

would be 217°C. 

In a given experiment, only part of these large adiabatic drops will be realized due to 
heating of the granule during the desorption by radiation and conduction. If the time 
constant for desorption (e.g. , time for 50% desorption) is much shorter than the time 
constant for radiative and conductive heating, then most of the adiabatic drop will be 
observed. Conversely, if the desorption time constant is much longer than the 
radiative heating time constant, then very little temperature drop will occur. 


The equation of heat balance which determines this is 


AH 


- P (ird 3 /6)C (dT/dt) = f W’/«) 22 ^ (- dq/dt) - (*d p )K(T Q - T) (6.17) 


where 


Pp = sorbent particle density -0.86 

(g/cm 3 ) 

= effective sorbent particle diameter 

(cm) 

P 

dT/dt = cooling rate 

(°C/sec) 

-dq/dt = desorption rate 

(ml STP/g-sec) 

K = rate of heat input to sorbent from surroundings 

per unit area and per unit temperature difference 
due to radiation and conduction 

T 0 = temperature of surroundings 

(°C) 

C D = heat capacity of sorbent 

(cal/g°C) 

s 

AHj = heat of desorption 

(kcal/mol) 


The left side of Eq. 6. 17 is the rate of heat loss from a single granule from the stand- 
pointof heat capacity and temperature drop. The right side is the net difference 
between the energy loss due to evaporation and the energy gain due to radiation an 
conduction from the surroundings. 
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At pressures of 10 torr and lower, which are typical during vacuum desorption 
calculation shows that gaseous heat conduction is negligible as compared with radiation 

" In ‘” B ^ Va ' Ue ° f K to * he equation. Thus, K is determined by ’ 

the relation J 


^( T 0 T) = eA(T Q - T 4 ) k 4eA T q(T q - T) 

where 

e = emissivity of sorbent surface (assume = 1) 

X = Stefan-Boltzmann constant = 5.67 x 10 -5 g/cm-sec 3 °K 4 
K = as above, but in cgs units (erg/cm 2 sec°K) 

Equation 6.17 may be rearranged to 

_ dT _ / dq\ 

' <Jt r ar/ ■ D < T o - T > 

where 


(6.18) 


(6.19) 


A 


^d 
22.4 C s 


( 6 . 20 ) 


24 eXT3 

° = X 2 * 39 x 10 ( 6 . 21 ) 

P p s ' 

In Eq. 6.21 D is the inverse of the time constant for radiative heating of a sorbent 
granule. For 0.1-cm particles with the emissivity e taken as 1.0, this time constant 
works out to be about 33 sec. Since the time for 50% desorption at temperatures above 
room temperature is of this order or smaller for both butane and toluene, it appears 
that a significant adiabatic cooling effect can be expected. The hill temperature history 
of a granule during desorption can be obtained by integrating Eq. 6. 19. Since this is 
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linear and of first order, it may be immediately integrated. It is eonvenient to eliminate 
time in favor of q by the relation 


dT _ dT dq 
dt ” dq dt 

The solution of Eq. 6. 18 then becomes. 

= «p [d j - ) 1 *i] / e *p [- D J - © d( i] d (^) (6 - 22) 

where (AT) ad is the adiabatic temperature drop equal to Aq Q which would be 
observed when D = 0 (i.e. , for a pellet completely thermally isolated from its 
surroundings). Numerical integration of Eq. 6.22 can be carried out using the 
experimental data for - dq/dt. 

For first-order desorption kinetics, i.e. , 

= V (6 ' 23) 

Equation 6.22 may be integrated in closed form to give the simple result 

T o- T 1 [/ jl\ d ' _ -a.1 (6 . 24) 

(AT) ad 1 - DAdlkJ 

Equation 6.24 confirms the intuitive guess that the main factor m determining the 
adiabatic cooling effects is the ratio of time constants, D/kj • In Fig. 6-17, Eq. 6.24 
is plotted, each curve corresponding to a different value of DA d • 
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Figure 6-17 indicates that significant adiabatic cooling effects can be expected even 
when the radiative time constant is four times larger than the reaction rate time 
constant. This curve (for DA d = 4) rises to only about 15% of the full adiabatic 

cooling effect, but for toluene this amounts to over 30°C of cooling since (AT) 
is over 200°C. 


The value of E>A d has been calculated from the actual data representative of those 

runs used for evaluation of the activation energy. For n-butane, DA d is about 

3.5 (Run No. M19 27). This results in a maximum temperature drop of around 4°C. 
This seems hardly enough to explain the large discrepancy between the desorption 
activation energy and the isosteric heats. Nevertheless, it is obvious that the adiabatic 
cooling even for n-butane is not entirely negligible. 

In the case of toluene, DA d is about 6. 3 (Run No. M18-48). The maximum tem- 
perature drop in this case, however, is about 20°C, since the adiabatic temperature 
drop is much higher for toluene (200-C). Such an error in the temperature would most 
assuredly affect the activation energy. Thus, for toluene the adiabatic cooling is 
obviously sufficient to account for the low apparent activation energies. Before 
applying Eq. 6. 21 in detail to the experimental data however, it may be apparent to 
perforin a relatively simple experiment in which a thermocouple is inserted directly 
into a charcoal granule and the cooling effect during desorption is actually measured. 
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Section 7 


VACUUM DESORPTION FROM ACTIVATED CARBON BEDS 

Sections 5 and 6 dealt with the effects of vacuum upon single particles of activated car- 
bon under various conditions of temperature, particle size, and contaminant species. 

In order to extrapolate these effects from the single particle case to the packed bed 
situation, it was necessary to obtain some typical data on desorption from such beds. 
To this end, a number of experiments were conducted in which, after saturation of the 
bed with one or more contaminants in a flowing stream of air, the bed was subjected 
to vacuum for varying lengths of time. In early experiments, the effectiveness of this 
procedure was measured at only one point, after a given vacuum desorption time, 
by desorbing the remaining contaminants in a gas stream, usually at elevated tempera- 
tures. Later refinements of the experimental method permitted determination of the 
amount desorbed for more than one contaminant at a number of points in the vacuum 
desorption, thus condensing the work of many long, tedious experiments into single 
ones. In addition, a number of experiments on diffusion of gases through such beds 
were performed in order to have necessary data to test various desorption models 
suggested below. 

Using the data presented in Section 5 and 6 and in this section, it was then possible to 
formulate a model for desorption of contaminants from packed beds. Using this model, 
it should then be possible to extrapolate desorption data from a few key experiments 

to the case of practical -sized packed beds for use in spacecraft contaminant removal 
systems. 

7.1 EXPERIMENTAL TECHNIQUES 

The initial vacuum desorption system used in Phase 1 is shown in Fig. 7-1. The sys- 
tem between the bed and diffusion pump was constructed of 5/8-in. OD copper tubing 
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Fig. 7-1 High Vacuum System Used in Phase 1 (see Fig. 2-1) 



with sweat-soldered fittings in order to provide high conductance for vacuum desorptior 
of the bed. To simulate space vacuum, the system consisted of a Veeco 2-in. oil- 
diffusion pump with a liquid nitrogen trap, plus associated valves and roughing pump. 
Vacuum desoiption pressure was monitored with an NRC Model 530 Alphatron gauge , 
whose minimum pressure sensitivity is 1 x 10 -5 torr. 

To increase data collection efficiency, the vacuum desorption system was modified in 
Phase 2 of the present work. The 5/8-in. OD copper tubing and 2-in. Veeco diffusion 
pump were replaced by 1-in. OD glass tubing and a 4-in. oil diffusion pump (see 
Fig. 7-2). The system also contained two large glass cold traps, connected in 
parallel, to allow continuous collection of all desorption products. All the pertinent 
bed dimensions, such as length, diameter, and particle size, are summarized in 
Table 7-1. 


During a vacuum desoiption experiment, pressure at the high vacuum (downstream) 
end of the bed is measured by a Hastings DV-6 gauge (down to 1 micron) and by a 
Veeco RG-3 ionization gauge below 1 micron. At the high pressure end of the bed 
(ipstream end), pressure is monitored by a Hastings DV-6 gauge. 

7.1.1 Single-Component Desoiption 

The first step in the experimental procedure, leading to vacuum desorption, is 
saturation of the bed with some contaminant to a known amount (see section 2. 1). At 
the conclusion of the flow adsorption, the bed is isolated by closing two high-vacuum 
valves upstream and downstream of the bed. 

In Phase 1 of the work, the vacuum desorption was performed by pumping, first 
with the roughing pump and then with the oil diffusion pump , according to standard 
vacuum system procedure. Desorption time was measured from the time at which 
the roughing valve was first opened. The desorption step was terminated by closing 
the valve to the diffusion pump. To determine the amount of material vacuum desorbed 
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Fig. 7-2 High Vacuum Desorption System 



Table 7-1 

DESORPTION HALF-TIMES 


Run 

No. 

■ 

1 H 



t 

1/2 



mi 

T 

(°K) 

L 

(cm) 

Species 

Exp. 

(min) 

Calc 

(min) 

D 

(cm) 

Remarks 

F5-1 to -6 

0.05 

298 

10.9 

C 3 H 8 

575 

430 

0.493 


F-9-1, 
Fll-1 to -3 

.05 

298 

4.8 

D ~ C 4 H 10 

3150 

1170 

.493 


F14-1 

.2 

298 

4.6 

n - C 4 H 10 

270 

226 

.493 

Quartz Chips 

F13-2 W 

.2 

298 

4.45 

"- C 4 H 10 

520 

446 

1. 65 


F15-l< a) 

.2 

298 

4.5 

n - C 4 H 10 

520 

446 

1. 65 


F15-2 (a) 

.2 

298 

4.5 

n - C 4 H 10 

520 

446 

1.65 


F16-1 

.2 

298 

8.0 

n - C 4 H 10 

1650 

1435 

1. 65 

Long Bed 

F17-1 

.2 

379 

4.4 

C 7 H 8 

1000 

1240 

1.65 


F17-2 

.2 

380 

4.4 

a “ C 4 H 10 

32 

31.4 

1. 65 


F17-3 

.2 

380 

4.4 

C 7 H 8 

1060 

1240 

1. 65 


F17-3 

.2 

380 

4.4 

n-C 4 H 10 

9 

8. 5 

1. 65 


F17-4 

.2 

340 

4.4 

n - C 4 H 10 

6 

3.5 

1.65 


F17-4 

.2 

340 

4.4 

- C 7 H 8 

4200 

3100 

1.65 



(a) All three runs averaged. 
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some time t , the bed was completely flow desorbed as described in section 2.1.1, 
Mi the amount vacuum desorbed was determined by taking the difference between flow 
Sorption and flow desorption. This same procedure had to he repeated several tunes, 
„ ,th a different value for time t , to establish a full curve for the amount desorbed 
v* ith time (see Fig. 7-3). To eliminate this tedious repetition procedure, it was 
decided to modify the vacuum desorption system for the Phase 2 work, so that one 
complete desorption curve was obtained with each run. The experiment^ procedure 
4or preparation of the bed and the amount adsorbed were the same in both cases 
The main difference occurred during vacuum desorption. The rate of vacuum de orp- 
r; on was continuously measured by condensing any desorbed contaminant in one o 

-f-he twoliquid-nitrogen traps (see Fig. 7-2). After a measured time interval, e 

noli trap was isolated by two stopcocks, and the other trap, in parallel, was ope 
+0 the vacuum lines. After isolation of the first trap, it was warmed up and any 
Condensed material was transferred into the calibrated volumetric system by use 
Liquid nitrogen, and isolated (see Appendix O). 

if the contaminant was, say, n-butane, the volumetric system was allowed to warm up 
*nd the pressure measured by the manometer (see Fig. 7-2). In cases w e * e ™ J 
pressure gases were used, such as toluene, and the contaminant was a hqmd at room 
temperature , a special small-bore stopcock and associated tube were used. In t 
procedure, the condensed material was isolated from the system by closing the sma 
Stopcock; then the tube and stopcock were removed from the system and weighed on an 
analytical balance. Since the tare weight was predetermined, the amount of maten 
Condensed could be calculated (see Appendix P). 


*7.1.2 Two-Component Desorption 


All the work on mixed desorption from packed beds was done in Phase 2 
of the new desorption system made this possible. 


The versatility 


The general procedure for mixed adsorption was 
component, flow-adsorption vacuum desorption. 


similar to that used earlier for single - 
A mixed stream of 26. 9 ppm n-butane 
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and 26.8 ppm toluene in 50% oxygen/50% nitrogen was admitted to the adsorption bed. 

The average flow rate through the bed was 936 ml RTP/min, and was very close to the 
flow rate used in previous single adsorption experiments with butane and toluene. 

The general procedure for obtaining vacuum desorption data was similar to that used | 
in previous desorptions; namely, the contaminants desorbed from the bed were I 

condensed in a large liquid nitrogen trap. To obtain a desoiption point, the trap was 1 

isolated from the main desorption stream and warmed up to room temperature. The 1 

trap was then opened to a manometer system and all the contaminants were recon- 
densed into a small trap with a stopcock. The trap was slowly warmed up to -22 °C 
(carbon tetrachloride slush) and held at this temperature until the butane had expanded 
to its maximum volume (pressure). The above-mentioned temperature was chosen 
because at that temperature the vapor pressure of n-butane is about 400 torr while 
that of toluene is less than 10 torr. Thus, a fairly efficient one-stage flash evapora- 
tion is achieved. Good mass-balance results obtained with toluene indicated that 
reasonable separation was achieved by this procedure. Finally, the stopcock on the 
collecting tube was closed and the tube with the toluene was weighed. 

7.1.3 Three -Component Desorption 

To analyze for the composition of the vacuum desorbed material, the following pro- 
cedure was adopted. As before, the sample was collected in a liquid nitrogen-cooled 
cold trap in the pumping line. This main trap was warmed to room temperature and 
was then opened to an auxiliary removable cold trap. After condensing in the 
removal trap, the sample was transferred to a previously evacuated 5-liter calibrated 
flask. All the condensed gases were expanded in the known volume and the total 
pressure brought up to 1 atm with helium. After allowing sufficient time for homo- 
geneous mixing of all the gases, 1-ml gas samples were taken and injected into the 
chromatograph. Using this procedure, it was found that all the Freon-12 was vacuum 
desorbed in just a few minutes and all the n-butane was desorbed within 30 min. 
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7.2 RESULTS ON SINGLE CONTAMINANTS 

All the data on vacuum desorption are summarized in Table 7-1. 

The total desorption curve at 25 °C was described by the flow experiments F5-1, 

-2, -3, -4, -5, and -6, and the data are shown in Fig. 7-3. An initial theoretical 
model for the desorption process based on Knudsen pore diffusion, which was proposed 
in the Phase 1 report (Ref. 23) , 

oo 

q^/q^ = \ (2n + l) -2 exp - [ k't(2n + 1) 2 | C 7 - 1 ) 

* n = 0 


2 2 

a constant = 7r /4 t£' 
amount remaining adsorbed 
amount originally adsorbed 
an integer 
time (sec) 

an effective diffusivity (cm /sec) 
reciprocal of a desorption rate constant (sec) 
effective bed length = £(1 - e) (cm) 
void function 

was used to predict a desorption rate curve. A comparison of this curve with experi- 
mental data is shown in Fig. 7-4. 

An initial set of experiments was conducted on the vacuum desorption of n-butane from 
a 30 x 40-mesh (0.05 cm) charcoal bed at 25°C. This work was done in Phase 1, and 
the data are shown in Fig. 7-5. To verify the effect of particle size as predicted by 
the theory of desorption described in section 7. 5 , more desorptions were run using 
n-butane adsorbed on an 8 x 12-mesh (0. 2 cm) charcoal bed. The consistency of 


where 

k* = 



n = 
t 


D k - 

T — 

V = 
€ = 
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Fig. 7 3 Vacuum Desorption of Propane From Packed Bed of BD Activated Carbon 
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Fig. 7-5 Vacuum Desorption of n-Butane From Packed Bed of BD Activated Carbon 


the results was checked by performing the desorption experiment three times. The 
results are shown in Fig. 7-6, and it can be seen that the results were quite repro- 
ducible. In agreement with the theory, the larger particles desorbed more rapidly. 
Sample calculations are shown in Appendix P. 

To test the theory of desorption from packed beds further , other experiments were 
performed. In run F14-1, the bed was diluted to one-half by volume with nonadsorbing 
quartz chips. This tested the term S Q in the desorption equation (Eq. 7.8). Since 
S is the amount initially adsorbed per unit volume, if the bed is diluted,the term 
should decrease by the same factor, and the total half-time for desorption should 
also decrease. All of this was found to be true quantitatively in experiment F14-1 
(see Fig. 7-7). 

In another test of the desorption equation, the bed length was increased from 4. 5 
to 8. 0 cm in Run F16-1. This would test the dependence of bed length to desorption. 

From theory it was predicted that the desorption half-time should increase with the 
square of the length. This was confirmed by experiment F16-1 (see Fig. 7-8). Further- 
more, once a desorption curve is established from experiment, the other curves 
could be predicted as shown in Fig. 7-8. 

Another test of the desorption theory was the variation of bed temperature. In flow 
run F17-2, the bed temperature was increased to 107 °C. Since the diffusion coef- 
ficient, Dg. , in the desorption equation is a function of the square root of absolute 
temperature, and the adsorption isotherm for n-butane on charcoal was determined 
independently on the microbalance , a prediction of desorption half-time could be made. 
The results and theory are in quite good agreement (see Fig. 7-9 and Table 7-1). 

In further investigating the desorption characteristics of packed beds, it was of inter- 
est to look at desorption of a heavy contaminant such as toluene. Furthermore, it 
was decided to do it at 107 °C to keep the amount adsorbed to a reasonable value. At 
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Fig. 7-8 Vacuum Desorption of n-Butane From Packed Beds of BD Activated Carbon 
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Fig. 7-9 Vacuum Desorption of n-Butane From a Packed Bed of BD Activated Carbon 


room temperature, the total adsorption and desorption cycle would have taken weeks. 
This experiment was used as a control or background data to determine the effects in the 
following experiments, when the coadsorption of n-butane and toluene was investigated. 
That was the purpose of flow run F17-1, and the desorption data are given in Fig. 7-10. 
The desorption was stopped after 2750 min, when only about 30% of the original toluene 
remained on the bed. There was little point in continuing the experiment because the 
main objective was to experimentally define the desorption half-time (50% remaining), 
and this was accomplished. 

7.3 RESULTS ON MULTIPLE CONTAMINANTS 

7.3.1 Multiple-Contaminant Desorption: Phase 1 

The first multiple-contaminant desorption studies were conducted in Phase 1 of the 
Contract. They were the vacuum desorption of the propylene /viny lid ene chloride pair 
and the vinylidene chloride/toluene pair. The results were inconclusive due to the very 
fast equilibrium obtained in the case of propylene and due to suspected leaks in the 
system in the experiments with toluene. In addition, some difficulty was experienced 
in getting satisfactory analysis of vinylidene chloride using the gas chromatograph. 

Thus, it was decided to perform similar experiments with other contaminant pairs 
more amenable to gas chromatographic analysis. 

7.3.2 Toluene/n-Butane Desorption at 107°C 

The first meaningful results on the simultaneous desorption of .coadsorbed species 
were achieved in. run F17-3. The desorption characteristics of each of the individual 
contaminants were examined at 107°C in flow runs F17-1, -2, and -3. The effects of 
coadsorption on vacuum desorption were also observed. It was found, as expected, 
that there was blockage of n-butane by toluene. There was about an 80% decrease in 
the amount of n-butane adsorbed in run F17-3 compared with run F17-2. The only 
difference was the presence of toluene at one-half its MAC concentration. Due to the 
lower initial amount of butane adsorbed, the desorption half-time for n-butane was 
shorter in F17-3 than in F17-2, but the value was still predictable (see Fig. 7-11 and 
and Table 7-1). 
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Fig. 7-10 Vacuum Desorption of Toluene From a Packed Bed of BD Activated Carbon 
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The presence of n-butane did not hinder the adsorption of toluene in run F17-3. The 
;ame amount of toluene was adsorbed in runs F17-1 and -3, and the desorption of 
oluene was not affected by the presence of n-butane (see Fig. 7-12). Indications are 
hat each component is desorbed independently. 

r.3.3 Three Components: Freon-12, n-Butane, and Toluene at 67 °C 

The vacuum desorption rate of a ternary mixture of Freon-12, n-butane, and toluene 
from a bed of BD charcoal has been measured at 67 °C. This temperature was chosen 
oecause at higher temperature the adsorption of Freon-* 12 and of n-butane at MAC 
becomes negligibly small (due to blockage by toluene) while at lower temperatures 
the rate of toluene desorption becomes too slow to measure in any convenient experi- 
mental time. At 67 °C the rate of Freon-12 desorption was too fast to measure (time 
for 50% desorption was less than 1 min) , while the half-times for n-butane and toluene 
were 6 min and 4200 min, respectively (see Fig. 7-13 and 7-14). These values are in 
agreement with the theory presented earlier, although agreement is less exact than 
found in other runs, 

7.4 DETERMINATION OF Dk , KNUDSEN DIFFUSIVITY , BY MEASURING 
PUMPING SPEED THROUGH A PACKED BED 

7.4.1 Experimental Procedure 

To determine the diffusion coefficient of n-butane through a packed bed of charcoal, two 
experimental procedures were available: 

(1) Since butane is adsorbed on charcoal, in order to measure °K directly for 
butane, the charcoal bed would have to be discarded and a nonadsorbing bed 
(say, q u artz chip) of equal mesh size would have to be substituted. 

(2) If the charcoal bed were used, then a nonadsorbed gas (helium) would have 
to be used. The diffusion coefficient for butane in the bed could then be 
calculated (from the helium value) , since Knudsen diffusion depends only on 
the inverse square root of mass when converting from one gas to another. 
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FRACTION REMAINING ADSORBED 



Fig. 7-12 Vacuum Desorption of Toluene From a Packed Bed of BD Activated Carbon 



FRACTION REMAINING ADSORBED 



Fig. 7-13 


Vacuum Desorption of Toluene From a Packed Bed of BD Activated Carbon During Concurrent 
Desorption of n- Butane 
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Since it appeared important to use exactly the same bed for the pumping experiments 
as for the desorption experiments, it was decided to use alternative (2), above: i. e. , 
to use the same charcoal beds as were used in the butane desorption experiments and 
to use helium, and as a double check, nitrogen. To obtain an evacuation time constant 
of reasonable magnitude, a 5-liter bulb was added at the end of the adsorption bed 
furthest from the vacuum pumps (see Fig. 7-15). 

In a typical experiment the charcoal bed is pumped to hard vacuum. The bed and 5- 

liter volume are isolated from the pumps and filled with helium to a pressure in the 
-3 

micron (10 torr) range. The main valve between bed and pumps is then opened and 
the decrease in pressure upstream of the bed is recorded as a function of time. Typical 
data are shown in Fig. 7-16. 

Experiments were carried out with helium at two temperatures (25 and 300 °C), for 
two particle sizes (0. 20 and 0. 05 cm) and for two bed lengths (4.4 and 8. 0 cm). The 
nitrogen data are shown in Fig. 7-17. Pressures were measured with a Hastings 
DV-6 vacuum gage which was calibrated with a Datametrics Barocel pressure trans- 
ducer. After correction with the pressure transducer, plots of log p versus time were 

_3 

linear over almost the whole pressure range (see Fig. 7-16 and 7-17). Above 10 
torr, some curvature of the plots would be expected due to the onset of viscous flow, 
but this is not observed in the present data. As expected from theory, the runs at 
higher temperatures are slower than those at room temperature (i.e. , t-jy 2 is larger), 
since the 5-liter volume is kept at room temperature. When proper correction for 
this is made, the measured diffusion coefficients are proportional to the square root 
of temperature, as predicted by theory. 

7.4.2 Theory and Results 

At pressures in the micron range, gases will flow through packed beds by Knudsen flow, 
which is identical to Knudsen diffusion. By measuring the pumping speed through a 
packed bed, may be measured directly. If a large volume, V, is connected to 
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Fig. 7-16 Diffusion Studies of Helium Through Packed Beds of BD Activated Carbon 
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Fig. 7-17 Diffusion Studies of Nitrogen Through Packed Beds of BD Activated Carbon 








a fast vacuum pump by a packed bed of length, L , with cross-sectional area, A 
then the rate of evacuation of the volume, V , is given by 


where c is the concentration of gas in the large volume. The left side of Eq. 7.2 is the 
rate of loss of molecules from V , and the right side is the rate of flow of molecules 
through the bed by Knudsen flow. It is assumed that the pumps are last enough so that 
the downstream pressure is effectively zero (i.e. , much less than the upstream pres- 
sure). Thus, the concentration gradient in the bed, Ac/Ax , is assumed linear and 
equal to c/L which c is the upstream concentration in the volume V. It is assumed 

in Eq. 7. 2 that the pressure is low enough so that Knudsen and not streamline flow 
occurs in the bed. 

It follows immediately from Eq. 7.2 that the time constant from evacuation of the 
large volume is 


VL (7.3) 

AD K 


The half-time for evacuation is In 2 times this, or 


t l/2 


0.693 


VL 

AD K 


(7.4) 


Thus, D„ may be computed from Eq. 7.4, so that 
K. 


D 


K 


0.693 


VL 

At l/2 


(7.5) 
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Frorn the linear plots of Fig. 7-16 and 7-17, half-times were taken and diffusion 
coefficients were calculated from the above equation. The results are collected in 
Table 7-2. The data are correlated by calculating a "tortuosity factor, " y. , from 
each D k , defined by the equation (see section 7.5) 1 ’ 


(7.6, 

where the term in parentheses is the Knudsen diffusion coefficient for a cylindrical 
capillaty of diameter d p . The factor e takes into account that only a fraction e of 
the bed is open for diffusion, while y takes into account all other factors. 

The tortuosity factors are also shown in Table 7-2. The constancy of y. for a given 
particle size (but changing T and molecular weight, shows that the (T/M)V2 i aw is 
obeyed quite accurately. However, the smaller particles have a smaller tortuosity 
factor (6. 3 versus about 9. 0) which is hard to explain. This might be due to changing 
particle shapes as the charcoal is ground smaller (the granules start as quite flat 
plate like structures before grinding) or to random errors in packing a particular 
bed. The onset of viscous flow would be in the opposite direction , to make the large 
particles have an apparent faster pumping speed, and, hence, a lower tortuosity 
factor, , than the small particles. 


7. 5 THEORY AND DISCUSSION OF RATE OF VACUUM DESORPTION FROM A BED 

To correlate and predict the experimental results presented in sections 7. 1 through 7 4 

t ree theories of bed desorption must be considered. Two of these , fortunately , can 
be eliminated. J 


The first and simplest theory is that diffusion of adsorbate through the bed is very 
fest, so that all granules in the bed are surrounded by high vacuum. AH granules 
would then desorb at the same rate as in the single -particle (microbalance) experiments 
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Table 7-2 

DIFFUSION PARAMETERS FOR BD CHARCOAL BEDS 



Dk (cmV sec) 

Yi 



d = 0. 05 cm 
P 

d = 0. 2 cm 
P 

d = 0.'05 cm 
P 

d =0.2 cm ' 
P 

Helium, 25 °C 
300 °C 

Nitrogen, 25°C(calc.)l 
150°C 
300°C 

128 

193 

49. 1 
60. 2 
73.5 

338 

592 

144.5 

177 

220 

6.24 
6.32 

6.62 

6.60 

6.25 

9.76 

8.22 

8.98 

9.05 

8.72 

n-Butane, 25 °C 

Calc, from He Dk 
C alc, from N£ 

33.6 
34. 1 

88.8 

100 

6.25 

6.55 

10.0 

8.93 


The time for half-desorption for the whole bed would then be of the order of seeonds , 
similar to the single-particle experiments. This mechanism may be eliminated 
immediately since the time for 50% desorption is on the order of days , not secon s. 
Thus, the rate of diffusion of adsorbate through the bed must be much lower than 
desorption from a single particle. Two mechanisms of diffusion through the bed 
must next be considered. 

in the Phase 1 Report (Ref. 23) for this contract, it was hypothesised that the main 
mechanism of transport in the bed was diffusion through the pore structure of the 
charcoal (see pp. 4-24 through 4-27 of Ref. 23 and also Fig. 7 4 of this report . 

This hypothesis, however, leads to an incorrect dependence of desorption rate on 
charcoal granule sines and so may be discarded. The experimental evidence is 
clearly that desorption rate is proportional to granule size, whereas a pore diffusion 
mechanism would predict either no dependence on particle size or inverse depen ence 
The most likely mechanism of adsorbate flow through the bed appears to be free- 
molecular Knudsen flow, or diffusion, through the void spaces of the bed, since the 
pressure is so low that the mean free path is larger than the channel dimensions , 
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which are about the same size as the charcoal granules. (The mean free path of n- 
butane at 10-micron pressure is about 0. 5 cm , which is larger than the granule size 
used here. ) If it is assumed, then, that Knudsen gas flow in the intergranular spaces 
is the main mechanism of butane transport in the bed, then the equation of diffusion 
accompanied by desorption to be solved is 


2 

t ~\ 9 £ 9S 

K 9x 2 " df (7.7) 

where 

c = gas concentration (mol/cm 3 ) 

S = adsorbed concentration per unit bed volume (mol/cm 3 ) 
x = distance from high-vacuum end of bed (cm) 
t = time (sec) 

D k = Knudsen diffusion coefficient in bed as defined by the equation (cm 2 /sec) 


where 

F x = flow rate per unit bed cross section at x (mol/cm 2 sec) 

InEq. 7. 7 the term 0c/8t has been neglected, since S is about 10 5 times greater 
than c . Equation 7. 7 is an equation of mass balance which states that the rate at 
which molecules are added to the flowing stream of adsorbate at a given point in the 
bed must equal the rate of desorption at this point. 

In Appendix Q, it is shown that the solution of Eq. 7.7 gives for the time of half- 
desorption 



K 


T 2 S n 

E O 

D c 
K o 


(7.9) 
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where 


L = bed length (cm) 

2 

D = gas-phase diffusion coefficient of adsorbate in the bed voids (cm / sec) 

K 3 

S = initial adsorbed concentration of adsorbate (mol/cm ) 

o 3 

c q = gas-phase concentration of adsorbate in equilibrium with S Q (mol/cm ) 

K = constant of the order of 0. 33 which depends on the shape of the adsorption 
isotherm 

A correlation described in Appendix Q for K gives the equation 

(K)" 1 = 5.16 - 3. 5(NLF) - 1.66(NLF) 10 (7-10) 


where 

(NLF) = simple function describing the degree of nonlinearity of the adsorption 
isotherm 

It may be defined as 


(NLF) = 1 - 2(c/c 0 ) l/2 


(7.11) 


where 

(c/c 0 )jy 2 = value of c/c q , which makes S/S Q = 0.50 

To determine (c/c D )jy 2 , the adsorption isotherm is plotted in reduced units, 

S/S versus c/c , where S is the initial adsorbed concentration of adsorbate 
o o o 

and c q is the gas-phase concentration in equilibrium with this. 

The value of D„ , the free -molecular-flow diffusion coefficient of adsorbate in the 
bed voids, has been measured experimentally by measuring the pumping speed through 
the bed (section 7 . 4) and this gives 
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(7.12) 


D k = ^.(9.7 XIO 3 )^!^ 

where 

e = void fraction in the bed 
T = temperature (°K) 

M = molecular weight (g/mol) 
dp = particle size (cm) 

7 1 = tortuosity factor equal to about 7. 0 (see Table 7-2 for details) 

Thus, from Eq. 7. 9 to 7. 12 the time for 50% desorption can be computed for any given 
adsorbate and at any temperature. 

Although Eq. 7 . 9 gives only the half-time for desorption, the entire course of desorp- 
tion (i. e. , a plot of percent desorbed versus a reduced time parameter, t/t^^) 
appears to follow the same pattern for both n-butane and toluene. In theory, the shape 
of this curve should depend somewhat on the shape of the adsorption isotherm , but in 
practice the shape seems to be similar for butane and toluene and thus one solution to 
Eq. 7.7 is common for both sorbates. If this is generally true, the common solution 
which is given in Table 7-3 and plotted in Fig. 7-18 is a general solution. 

Table 7-3 

PERCENT DESORBED VERSUS REDUCED TIME, t/t^ 


% Desorbed 

t/tl/2 

0 

0 

25 

0.25 

50 

1.0 

70 

2.7 

80 

4.0 

90 

8.0 


172 




Thus, from Eq. 7.9 plus Fig. 7-18, the complete course of desorption can be com- 
puted. This has been done for curve 2 in Fig. 7-8 which depicts the average desorption 
rate for three n-butane runs. The tortuosity measured by the helium transport 
experiment, a K-value of 1/3, and the value of S Q /c o from the experimental isotherm 
were used to calculate the half-time. Then the generalized rate curve of Fig. 7-18 
was used to calculate the theoretical curve (dotted line in Fig. 7-8). The agreement 
is well within experimental error for the entire rate curve. 

In summary , the preceding theory when applied to the experimental data shows excel- 
lent agreement. (The data are summarized in Table 7-1.) The theory correctly pre- 
dicts that: 

(1) Desorption rate should increase with increasing particle size. Thus, the 
above table shows that 0. 20-cm particles desorb much faster than 0. 05-cm 
granules . 

(2) The dependence of desorption rate on bed length is correctly predicted as 
an inverse-square dependence. Compare Run F15-2 with F16-1 in Fig. 

7-8. 

(3) The absolute values of the times for 50% desorption are predicted almost 
within the experimental error, the total rate curves are reasonably well 
reproduced, and the effects of bed length and adsorption capacity are 
quantitatively predicted from the detailed rate curve from one set of 
conditions . 
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Section 8 

APPLICATION TO DESIGN PROCEDURE 


The ultimate purpose in obtaining the fundamental adsorption information on activated 
carbon has been to provide the necessary physical-chemical data to form a rational 
basis for a practical design procedure for sorption filters which can be used to control 
trace contaminants for extended missions. In this section, the information developed 
under the current contract will be examined from the design point of view. The pur- 
pose is twofold. First, the development of workable design procedures gives a 
meaningful evaluation of the completeness of current progress in adsorption theory 
and practice and permits identification of the most critical areas for future research. 
Second, the development of the design procedures is of intrinsic value for advanced 
planning and design of future life support systems and as a basis for construction of 
prototype sorbent filters for engineering evaluation. 

The first important question which arises in such a design is that of regenerative 
versus nonregenerative sorbent systems. Consider, for example, a 9-man, 2-year 
mission. The nonregenerative design would require either (1) that enough charcoal 
be launched to last for the 2-year mission or (2) that the charcoal be supplied at 
intervals, say, every 90 days. In either case, the total amount of charcoal is the 
same and a crude estimate of the amount can be made without recourse to the specific 
design. To make such an estimate, it is necessary to assume a reasonable maximum 
production rate for contaminants. One suggestion for this has been given as 50 g/day 
nonbiological and 20 g/day biological. For a 2-year mission, this would result in a 
total contaminant production of 51 Kg (112 lb) . Since it may reasonably be expected 
that the charcoal will adsorb on the average about 1% by weight, then the total charcoal 
required would be about 5000 Kg (11,000 lb). If the charcoal were replenished on a 
90-day basis, the 90-day load would be about 630 Kg (1400 lb). 
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From such estimates, it is clear that only a regenerative system could be seriously 
considered. In a regenerative system, there are four general areas of concern: 

(1) adsorption capacity of the sorbent compared with production rates of trace con- 
taminants, (2) adsorption dynamics, (3) regeneration, and (4) interfacing with the 
total life support system. 

8.1 ADSORPTION DYNAMICS 

One complication which is not always recognized in the design of adsorptive removal 
systems is that adsorption of heavy contaminants even in relatively small concentra- 
tions can completely block the adsorption of light contaminants. Thus, if the hydro- 
dynamics of the adsorber are such that the heavy components break through early in 
the adsorption cycle, then the surface can no longer adsorb lighter components. In 
fact, the light components which have already been adsorbed begin to desorb as the 
heavy components break through. Of course, this can happen when the removal 
efficiency for the heavy components is still favorable for their removal. Thus, it is 
a requirement of an operable adsorber that the concentration profiles for all the heavy 
contaminants be sharp enough to leave an appreciable portion of the bed unblocked so 
that the lighter contaminants can be retained. This is best accomplished through 
careful design of bed geometry and selection of particle sizes and linear velocities 
consistent with the bed geometry. In general, this displacement phenomenon requires 
larger length-to-diameter ratios than have previously been assumed in some designs 
This, of course, complicates the selection of a design with a low pressure-drop 

limitation but suggests substantial benefits in minimizing the total gas flow rate through 
the adsorber* 


A great deal of theoretical work and some experimental work directed toward estab- 
lishing the theory in various forms has been done on the problem of adsorption 
dynamics in packed beds by workers in ion exchange, gas chromatography, and 
chemical reactor design, as well as those primarily interested in adsorption dynamics 
and equilibria. Robell and Merrill (Ref. 24) have reviewed this work and examined the 
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various approaches with the objective of recommending procedures which give 
acceptable engineering predictions of adsorption dynamics. What is suggested is the 
use of the number of transfer units, , as the figure-of-merit of the adsorption 

bed and the Pec let number 

d u 

Pe = -E_ (8.1) 

ve D f 


where 

Pe = Peclet number 
dp = particle diameter (cm) 
u = linear velocity (cm /sec) 

Df = fluid-phase diffusivity of contaminant (cm /sec) 


: as the system correlation parameter. These have been related via a semiempirical 
; correlation by Vermeulen (Ref. 25), in which the quantity N TlJ d p /bL is plotted as a 
function of the Peclet number on log-log coordinates. A version of this correlation is 
: given in Fig. 8-1, in which the quantities b and L are an isotherm shape factor and 
the bed length (cm), respectively, while the term represents contaminant 

diffusivity (cm^/sec) in the sorbent pores. 


In the work of Robe 11 and Merrill (Ref. 24), it is shown that the correlation in Fig. 8 1 
is satisfactory for engineering purposes if the number of transfer units is taken as 



where 

t = breakthrough time (sec) 

b 

w = breakthrough width (sec) 
as shown in Fig. 8-2. 

In order for the breakthrough curve to be reasonably sharp (i.e. , sharp enough to 
eliminate poisoning of the sorbent for the lightly adsorbed contaminants), it is necessary 

to maintain a design such that 

N T 2= 100 
TU 
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Fig. 8-1 Adsorption Dynamics Correlation 





Fig. 8-2 Definition of t and w 

This can usually be obtained for Peclet numbers between 10 1 and 10 where the curve 
in Fig. 8-1 passes through its maximum. The inequality, Eq. 8-3, then is the first 
criterion for an acceptable design. 

8.2 DESIGN PROCEDURE 

A workable design procedure will be described in this section which can be used to 
give a near-optimum design of a trace contaminant sorbent bed for a large number of 
trace contaminants. It is essentially the procedure used to produce a design for the 
NASA/MSC 1973 Basic Subsystems Module (Ref. 26). Because additional information 
has been obtained since that design, particularly with respect to mixed adsorption and 
the dynamics of vacuum desorption from sorbent beds, the design approach will be 
outlined here with attention given to how the newly acquired physical-chemical infor- 
mation should be integrated into the design. 

The first step in the design is to define a minimum flow rate for acceptable operation. 
This is done by writing the mass balance for a contaminant in the cabin: 


m = Fc c t) 


(8.4) 
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where 


m = production rate of contaminant (mg/ sec) 

F = volumetric flow rate into the adsorber (1/sec) 

C c ~ c ^ ara cteristic concentration of contaminant (mg/1) 
t? = removal efficiency 

Using the equation above , the minimum concentrations for a given flow rate can be 
calculated by assuming 7? = 1.0 . The minimum flow rate is set by the condition that 
c c for all contaminants must be lower than the SMAC value. Violation of this criterion 
defines flow rate-limited contaminants, and once the minimum flow rate is determined, 
the cabin concentrations for all the other contaminants may be calculated from Eq. 8.4. 
From these concentrations, it is possible to calculate the adsorption potential, A , 
for each of the contaminants. 

If the design were to be done by machine calculation, then it would be possible to carry 
each component (for the BSM design, 150 trace contaminants were considered) indi- 
vidually through the calculation. For a hand-calculated design, it is necessary to 
make some simplifications at this point for calculational facility. 

One such approach used commonly in separation design calculations for complex 
hydrocarbon mixtures is the definition of a small number of pseudo-components usually 
grouped according to boiling-point range. In the case of adsorption, a reasonable 
criterion for grouping would be according to adsorption potential, A . The trace con- 
taminants can be so grouped into 10 or fewer pseudo-components according to their 
adsorption potential, with appropriately averaged molar volumes, adsorption potentials, 
boiling points, and molecular weights. Using these pseudo-components in their re- 
spective ratios and the generalized potential plot, a relative bed profile can be con- 
structed giving the contaminant concentrations q as a function of bed length L for 
each pseudo-component, assuming completely sharp breakthrough. (See Fig. 8-3 for 
the plot derived for the BSM design.) From such a generalized plot, the breakthrough 
characteristics of any combination of charcoal load and regeneration time could be 
determined. The arrow on the abscissa in Fig. 8-3 shows the design point selected for 
the BSM design, with bed length equal to 10 1 times the length L g needed for adsorption 
of all 10 pseudo-components in a given time. 
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Next, a bed geometry consistent with the minimum flow rate and allowable pressure 
drops is selected to give sharp breakthrough curves (i. e. , in accordance with Eq. 8.3). 
The cycle time is now set by the regeneration which is a function of desorption tem- 
perature and bed geometry. Having determined a cycle time, the amount of charcoal 
is determined from the design point on the generalized concentration profile. Now the 
selected geometry must be recalculated to be consistent with the charcoal loading, and 
the procedure iterated to closure. This procedure can be done as a function of desorp- 
tion temperature to form the basis of an optimization between heating and cooling 
weight penalties and charcoal loading. 

Once the bed has been sized in this way, the actual number of transfer units for each 
component can be determined and the actual bed profiles calculated using the multi- 
component adsorption correlations developed herein (see Section 4) . From such a 
multicomponent bed profile, the actual removal efficiencies for each pseudo-component 
can be calculated. From these and the design flow rate, the actual cabin concentra- 
tions can be estimated. If these are substantially different from those found initially, the 
whole design procedure must be iterated to closure. 

Once the base design is determined in this way, it is possible to complete parametric 
calculations for interfacing with other modules of the life support system. If such 
interfacing suggests major changes from the design base case, then a redesign should 
be completed to the new conditions. 

8.3 CRITICAL DESIGN CONTAMINANTS 

One of the most enlightening features of design calculations even as crude as those 
used in BSM design is the identification of the critical design parameters. Specif- 
ically, for the BSM design a number of contaminants were identified as possible 
critical contaminants, i.e. , those that might control the system design. The first 
are those which are flow rate-limited, meaning that the cabin concentration at 
100 percent removal efficiency is greater than the MAC value due to the production 
rate being too high for the given air flow rate. Next are those not removed because 
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they are too weakly adsorbed or because their production rate is too high for an 
acceptable design. For the BSM case there were two of the former and a number o 
the latter. 


Flow Rate-Limited Contaminants 

Phenol 
Pyruvic Acid 


Unremoved Contaminants 


Acetylene 
Carbon Monoxide 
"Freon-23" 
Hydrogen 
Methane 
Nitrous Oxide 
Sulfur Dioxide 


Nitrogen Dioxide 

Formaldehyde 

Chlorine 

Hydrogen Chloride 
Hydrogen Sulfide 
Nitric Oxide 
Ammonia 


Usually, the charcoal filter is integrated with some kind of catalytic oxidizer and basic 
absorber (e.g. , lithium hydroxide). Under such circumstances, only ammonia would 
not be removed by the total system. Thus, the production of ammonia which is mostly 
biological must either be substantially curtailed, controlled at the source, or removed 
in some other way than by sorption. 


Furthermore, the BSM design has 
near the design point. They are: 

Acetaldehyde 
Acetonitrile 
Acrolein 
Carbon Disulfide 
Carbonyl Sulfide 
Chlorofluoromethane 
Chloropropane 
Ethane 
"Freon-21" 


which break through 


"Freon-23" 
"Freon-125" 
Methylacetylene 
Methyl Alcohol 
Methylhydrazine 
Methyl Mercaptan 
Nitrogen Tetroxide 
Propylene 
Vinyl Chloride 


identified those contaminants 


The production rates and allowable concentrations of these contaminants dominate the 
design. Effort directed toward more careful definition of both of these for these con 
taminants will result directly in more reliable design. 
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Under this contract and previous NASA contracts, this laboratory has investigated 
nearly half of these critical components. Nevertheless, it would be well worth the 
effort to study this group in more detail, particularly with respect to the blockage 
effects for mixed adsorption. Also, no work has yet been done to identify those con- 
taminants which are the critical ones for vacuum desorption and such work is essential 
because of the controlling nature of the pump-out times. 
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Section 9 

ADSORPTION FROM HUMID AIR 

All the experiments covered thus far in this report have dealt with adsorption of pure 
contaminants, mixtures of pure contaminants, or contaminants in dry air (50% oxygen). 

In actual spacecraft operation, however, water vapor is always present in significant 
quantities. In order to complete an accurate study of adsorption processes in space- 
craft applications , it is of extreme importance to know the effects of water vapor upon 
the adsorption of contaminants. Only with such knowledge will it be possible to decide 
upon the optimum location of the sorbent bed in the spacecraft air purification cycle. 

To this end, a series of experiments has been performed with two typical contaminant 
types, one polar and one nonpolar, and at various relative humidities. 

9. 1 EXPERIMENTAL TECHNIQUE AND RESULTS 

The flow system described in Sections 2 and 3 was modified to include a gas-stream 
humidifier. Several methods of humidification were considered, and it was concluded 
that an efficient and convenient design was one in which the gas is passed through two 
water bubblers connected in series (see Fig. 9-1). 

The following matters were next dealt with: 

(1) Determination of a desirable method for attaining various partial pressures 
of water effluent from the bubblers 

(2) Consideration of the fact that one of the contaminant gases proposed (ethanol) 
is quite soluble in water 

It was decided to introduce the ethanol into the gas stream after humidification. 
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Fig. 9-1 Flow System Humidifier 



9.1.1 Design and Fabrication 


The humidifier system was designed and built as shown in Fig. 9-1. The two 
humidifier flasks and the ethanol solution flask were immersed in a constant- 
temperature water bath at 25 ± 0 . 02°C. 

9.1.2 Method of Establishing Relative Humidity 

Two possible methods were considered. One was to lower the temperature of the 
water reservoirs and thereby decrease the equilibrium vapor pressure of water. 

The main difficulty with this procedure was that, since the flow runs in general 
were quite long, it would have been difficult to maintain accurately a constant 
subambient temperature of, say, 10°C. The second method, which was the one 
adopted, was based on the fact that the addition of concentrated sulfuric acid to water 
could be used to establish any desired partial vapor pressure. A summary of pertinent 
data is given in Table 9-1. 


Table 9-1 

SUMMARY OF DATA USED TO ESTABLISH DESIRED HUMIDITIES 


Percent 
Relative 
Humidity 
at 25°C 

ph 2 o 

(torr) 

Weight Percent 
H 2 SO 4 
(Ref. 27) 

Weight H 2 SO 4 
Used (Ref. 28) 

(gA) 

100 

24.0 

0.0 

0.0 

95 

22.8 

10.0 

106.6 

70 

16.8 

31.5 

387.8 

45 

10.8 

46.0 

624.2 


One problem with this technique was that there was a continuous depletion of water, 
lost by evaporation into the gas stream. As water was removed, the concentration of 
sulfuric acid increased and the partial pressure of water decreased. This problem 
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was solved by choosing an initial volume large enough so that the loss of small 
amounts of water did not significantly affect the diluted acid composition. In the 
present case , it was estimated that a volume of 1 liter would experience a loss of 
water of 0. 1%/hr during an average run. 

To verify the predicted partial pressures, the system described above was pumped 
until pressure readings downstream from the second bulb were constant. This 
pressure was then taken as the partial water vapor pressure at the temperature of 
the bath surrounding the bulbs. 

In the flow experiments described in Sections 2 and 3, premixed gases containing 50% 
oxygen and 50% nitrogen and one or more trace contaminants at or near their MAC 
levels were used. In the experiments described here, contaminant introduction had 
to be modified due to solubility considerations. It was found from the literature that 
the solubility of n-butane in water is 0.025 ml/ml HgO (Ref. 29). Calculations 
indicated that a 1-liter charge of water solution would be saturated with n-butane 
from the premixed cylinder in approximately 2 hr at the flow rates anticipated. 
Equilibrium in this manner was performed before diverting the gas through the 
sorbent bed . 

In the case of ethanol as a contaminant, a premixed gas could not be used due to the 
mutual solubility of ethanol and water. To circumvent the problem of solubility, the 
system was designed so that the gas stream not containing ethanol was saturated to 
the desired humidity. Then ethanol vapor was introduced into the gas stream through 
a Granville -Phillips leak valve to the desired contaminant concentration (see Fig. 9-1). 

Using the above-described flow system modifications, the humid air experiments were 
accomplished. The results are summarized in Table 9-2. 
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Table 9-2 

SUMMARY OF ADSORPTION DATA FOR n-BUTANE, ETHANOL, AND WATER ON 

BD ACTIVATED CARBON 


Run No. 

Species 

Equilibrium 

Concentration 

c eq 

(ppm) 

Equilibrium 

Pressure 

P 

Vapor 

Pressure 

po 



(torr) 

(torr) 

F18-1 

n-C^HjQ 

h 2 o 

26.4 

21.2 

25.2 

F19-1 

11-C4HXO 

h 2 0 

24.0 

21.2 

28.3 

F19-2 

C2H50H 

H20 

60,0 

0 

19.6 

F19-3 

C 2 H 5 OH 

h 2 o 

35.5 

0 

25.2 

F19-4 

n-C4H 10 

h 2 o 

25.5 

21.2 

17.9 

F19-5 

C 2 H50H 

h 2 o 

20.0 

21.2 

26.7 

F19-6 

C2H5OH 

h 2 o 

28.0 

21.2 

19.6 

F19-7 

n-C4Hio 

H 2 0 

25.6 

15.0 

26.7 

F19-8 

n-C4Hio 

h 2 o 

25.8 

15.0 

18.7 

F19-9 

c 2 h 5 oh’ 

h 2 o 

27.0 

15.0 

26.7 

F19-10 

C2H5OH 

h 2 o 

16.2 

15.0 

19.6 

F19-11 

n-C4H 10 

h 2 o 

25.1 

10.0 

25.2 

F19-12 

n-C4Hio 

h 2 o 

26.2 

10.0 

20.5 

F19-13 

C 2 HsOH 

h 2 o 

- 24 

10.0 

25.2 

F19-14 

c 2 h 5 oh 

h 2 o 

20.5 

10.0 

19.6 


Amount 
Adsorbed 
q a # 

(ml liq/g) 


Amount 
Flow 
Desorbed 

(ml liq/g) 


1.9 x 1CT 3 

1. 14 x 10" 3 
3.58 x 10“1 

1.1 x 10' 2 


2.55 x 10" 2 

9.26 x 10" 3 
None detected 

5.56 x 10“ 3 

2.0 x 10" 3 
None detected 

7.43 x 10" 3 

8.41 x 10 -3 
None detected 

1.32 x 10“ 2 
1.39 x 10 -3 
2.50 x 10' 2 
8.12 x 10 -3 
1.8 x 10" 2 

2.6 x 10' 3 


1.3 x 10- 3 


~ 3 x 10" 1 
9. 8 x 10~ 4 

1.2 x 10‘ 2 

2. 2 x 10-2 
9.75 x 10“ 2 
6.47 x 10 -3 

1. 7 x 10-3 
7.34 x 10“ 3 
8.70 x 10" 3 
1.46 x 10“2 
2. 16 x 10~ 3 
2.58 x 10“2 

1.08 x IQ' 2 


Adsorption 

Potential 

A 

(mol°K/mlliq) 


16.1 

1.19 

16.5 

2.01 

23.5 
00 

17.5 

OO 

20.4 

16.7 

19.1 

1.6 

25.2 

17.5 

16.5 

4.0 

20.1 

19.8 

18.6 

4.0 

26.4 

20.1 

16.5 
6.4 

20.3 

23.7 

18.5 
6.4 

25.9 

23.3 


Temperature | 
T 

(°C) 


26 

28 

66 

26 

64 
27 
66 
27 

65 
27 

66 
26 
67 
26 
66 
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9.2 DISCUSSION 


Study of the Table 9-3 and comparison with previous data taken under dry conditions 
reveal the important role of relative humidity in causing blockage. Some typical 

values of reduction in contaminant adsorption capacity due to the presence of water 
vapor are shown in Table 9-3. 


Table 9-3 


CONTAMINANT BLOCKAGE DUE TO WATER VAPOR 


Run 

No. 

Contaminant 

Temp. 

(°C) 

Concentration 

(ppm) 

Relative 
Humidity 
at 25°C 

Blockage 

(%> 

F19-11 

n- Butane 

26 

25 

42 

28 

F19-7 

n- Butane 

27 

26 

63 

79 

F18-1 

n-Butane 

26 

26 

89 

95 

F19-12 

n- Butane 

67 

26 

42 

0 

F19-8 

n-Butane 

65 

26 

63 

0 

F19-4 

n-Butane 

64 

26 

89 

0 

F19-13 

Ethanol 

26 

24 

42 

r 5 

F19-9 

Ethanol 

27 

27 

63 

30 

F19-5 

Ethanol 

27 

20 

89 

65 

F19-14 

Ethanol 

66 

20 

42 

56 

F19-10 

Ethanol 

66 

16 

63 

73 

F19-6 

Ethanol 

66 

28 

89 

71 


An overwhelming influence of relative humidity is evident, and clearly there is an 
urgent need to predict such blockage phenomena as part of any realistic design 
procedure. To accomplish such prediction, use is made of the modified potential 
theory as expressed in Eq. 4.50, 

T 

A T = A + v~ logX l (4.50) 

m 


190 




Values of A are shown in Table 9-2, along with adsorbed quantities of n-butane and 
sthanol. It was not possible in these mixed adsorption experiments to measure 
precisely the amount of water actually adsorbed because of limitations in the detection 
system used. The adsorption isotherm for water on BD has been measured gravi- 
metrically, however (see Fig. 2-11), and it will be assumed that the amount of water 
adsorbed is not affected by the presence of the solute (i.e. , n-butane and ethanol). 
This assumption is likely to be inappropriate only at very high values of X. . Thus , 
it is possible to calculate the mole fraction. 


X, = 



(9.1) 


where 


V. 


q s 

q w 

^rns 

mw 


amount of solute adsorbed, obtained from humid air adsorption 
amount of water adsorbed, obtained from pure water isotherm 
molar volume of solute 
molar volume of water 


and 


= q s + q w 


(9.2) 


If this model for unmixed adsorption is valid, then it may be expected that a plot of 
q from Eq. 9.2 versus A T from Eq. 4.50 will closely fit the generalized potential 
plot. Data from such calculations are summarized in Table 9-4 and plotted in Fig. 9-2. 
The line in Fig. 9-2 is the generalized potential plot and the data are in excellent 
agreement except at the highest values of q T where it is likely that q w is somewhat 
less than the value obtained from the pure water isotherm. The uncorrected adsorption 
potential is plotted in Fig. 9-3 to show the great improvement obtained in using the 
modified A as in Fig. 9-2. Thus, it is possible to estimate humid air adsorption 
from the generalized potential plot and the water isotherm. Because this has been 
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Eig. 9 2 Potential Plot for Adsorption From Dry and Humid Gas Mixtures on BD 
Activated Carbon 
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AMOUNT ADSORBED (ml liq/g) 



A = ADSORPTION POTENTIAL (mol °K/ml liq) - T/V^ log 


Fig. 9-3 Potential Plot for Adsorption From Dry and Humid 
Gas Mixtures on BD Activated Carbon 
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Table 9-4 

modified potential-theory quantities for humid 

ADSORPTION EXPERIMENTS 


Solute 

% 

% 

q T 

X 

A^p 

A 

Run 

No. 

n- Butane 
n- Butane 
n- Butane 
n- Butane 
n- Butane 

9.3 x 10“3 
8.12 X 10-3 
8.41 x 10~ 3 
2.5 x 10-2 
7.43 x 10 -3 

1.5 x 10-2 

7.0 x 10-3 
1.4 x 10-2 
3. 2 x 10-2 

2. 0 x 10-1 

2.44 x 10-2 
1.51 x 10-2 
2.2 x 10-2 

5.7 x 10-2 

2. 07 x 10-1 

1. 04 x 10"! 
1.78 x 10“1 
1. 01 x 10"! 
1.28 x 10-1 
6.88 x 10-3 

17.0 
17. 7 
16.5 
13. 7 
10.3 

20.4 
20.3 
20.1 

16.5 
16.5 

F19-4 

F19-12 

F19-8 

F19-11 

F19-7 

Ethanol 

Ethanol 

Ethanol 

Ethanol 

Ethanol 

Ethanol 

■2.6 x 10 -3 
1.39 x 10-3 
2.0 x 10~ 3 
1.8 x 10-2 
1.32 x 10~ 2 

5. 56 x 10 -3 


7.0 x 10-3 

1.0 X 10-2 

1.2 x 10-2 

3.2 x 10 -2 
1.6 x 10-1 

4. 2 x 10 -1 

9.6 x 10“ 3 
1.1 x lO- 2 
1.4 x lO -2 
5.0 x 10“1 
1. 73 x 10 -1 
4.26 x 10"! 

9.7 x 10 -2 
3.87 x lO -2 
0.61 x 10-2 
1.41 x 10~1 
2.34 x lO' 2 
3.81 x lO -3 

20.4 
18.7 
17.9 

14.4 
10. 7 

7.4 

25.9 

26.4 
25.2 

18.5 

18.6 
19.1 

"F19-14 

F19-10 

F19-6 

F19-13 

F19-9 

F19-5 


successful for n-butane, a nonpolar molecule, and 

might even exhibit hydrogen bonding with water, it 
generally. 


ethanol, a polar molecule which 
is presumed that it may be applied 
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Section 10 

CONCLUSIONS AND RECOMMENDATIONS 


10.1 CONCLUSIONS 

• A successful correlation for the estimation of sorptive capacity has been 
determined. This correlation, which has been extended to adsorption 
potential theory parameters, applies to both pure and mixed contaminants, 
under both dry and humid conditions, and at various temperatures. The 
capacity correlation is a vital design tool needed to assess the effect of 
different operating conditions on equilibrium performance. 

• The rate of vacuum desorption from single sorbent particles and from sorbent 
beds has been characterized. Mechanisms and theories proposed have been 
experimentally proven to be adequate to describe the phenomena over the 
ranges tested. 

It has been shown in single-particle studies that the rate-determining step is 
desorption from the external particle surface. Thus, the rate is inversely 
proportional to the particle diameter. 

For vacuum desorption from sorbent beds, the rate-determining step is 
Knudsen diffusion through interparticle voids. The Knudsen diffusivity has 
been experimentally measured. Although the single-particle studies showed 
that smaller particles desorb faster, the theory correctly predicts the 
surprising result that the opposite is true for beds: beds composed of larger 
particles desorb faster than those containing smaller particles. The desorp- 
tion rate dependency on bed geometry, adsorption isotherm, and initial gas- 
phase concentration is also contained in the theory, and the theoretical 
predictions agree with experimental measurements. 

Neither experimental vacuum desorption rate data nor theoretical descrip- 
tions were available prior to this contract. In the course of this work, sub- 
stantial information on both fronts has been generated, sufficiently to begin 
design work on a prototype unit. 
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• The information gained on this program was instrumental in enabling the 
development of a quantitative design methodology for vacuum-regenerative 
sorbent systems. This design procedure is the first one known to consider 
realistically the phenomenon of adsorbate interference. 

10.2 RECOMMENDATIONS 

• Further adsorption capacity measurements with contaminant mixtures are 
recommended in order to provide further confirmation of the generality of 
the modified adsorption potential theory. 

• The design, construction, operation, and test of a prototype regenerative 
sorption system is recommended. 

• Implicit in a design procedure, the following steps should be done after 
theoretical determination of the "key components" (those contaminants which 
limit the design): For the key components only, (1) obtain accurate produc- 
tion rates, and (2) verify experimentally the capacity and the kinetics. 

• Existing sorbent system designs should be updated, based on the advances 
reported here. A computer program for the design procedure should be 
developed as part of this task. 

• The successful correlation of sorption capacities under humid conditions 
indicates the theory may be applicable to waste^water purification by adsorp- 
tion. Appropriate tests should be performed. 

• Efforts to establish more carefully the production rates of the design- 
controlling contaminants should be initiated. In such tests , some attempts 
should be made to evaluate the attenuation of the rates with times of the 
order of mission times. The effect of temperature on contaminant 
production rates from materials within the spacecabin should also 

be investigated. 
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Appendix A 

SAMPLE CALCULATIONS FOR FLOW EXPERIMENTS 


A. 1 TYPICAL CONDITIONS DURING CALIBRATION 


(a) Flow rate = 250 ml RTP/min 
< b > P system = °- 5atm 

(c) Sample loop volume (Beckman Valve) = 2. 0 ml 

(d) Range x attenuation of the electrometer = 10 x 8 

(e) Contaminant peak area for 2. 0 ml volume at 0. 5 atm total pressure 
2. 5 mV-sec as determined by the disk integrator 

(f) Electrometer transconductance = 4 picocoulomb s/mV- sec (= 4 nmho) 

(g) Premixed gas contaminant concentration = 50 ppm 

A. 2 SAMPLE CALCULATIONS OF DETECTOR OUTPUT 


2. 5 mV-sec x 10 x 8 x 4 pC /mV-sec 
Detector output - o. 5 atmx2 ml 


800 pC 
" ml RTP 


The total response of the detector to 50 ppm contaminant is thus 800 pC/ml RTP 
which is the response to a known amount of contaminant equal to 5 x 10 ml RTP. 

A. 3 SAMPLE CALCULATIONS OF EFFLUENT CONCENTRATION 

p = 0. 5 atm 

system 

Loop volume = 2.0 ml 
Range x attenuation = 10 x 2 
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Peak area 


1. 5 mV-sec 


Detector response = ^ m V sec x IQ x 2 x 4pC/mV-sec 

(for unknown cone. ) ^ a ^ m x 2 ml 

= 120 pC/ml RTP 


Comparison of unknown to known detector response 


Unknown Concentration 


unknown response 

known response X standard concentration 


" 800 X 50 PP m 

= 7.5 ppm 


A. 4 SAMPLE CALCULATION OF AMOUNT ADSORBED a 

’ ^a 


A typical breakthrough curve for propylene at C = 15 ppm 



The parameters established from the breakthrough curve are: 

(1) C = 15 ppm 

(2) t g = 30 min (t g = time such that area A = area B) 

(3) Flow = 500 ml RTP/min 
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/4^ p = 0. 52 atm (see below for calculation) 
v ; bed 


Parameters established from other sources: 

(1) Weight of bed: W = 0. 500 g 

(2) Molar volume of the liquid at its normal boiling point: V m = 66. 6 ml liq/mol 

(3) Gas constant at RTP = 24 400 ml RTP/mol 




C eff (PPm) 


10 


~ " /ml RTP\ 

W(g) x gas const. [ mol ) 


15 x 30 x 66. 6 x 500 x 10 
0. 500 x 24 400 


= 4. 92 x 10 


ml liq 
g (sorbent) 


Since there was a pressure drop across the bed, an estimate was made of the pressure 
in the middle of the bed. A typical calculation is given. 



P 


4 


P 

P 

AP 


1 

2 

3 


P 


4 


= pressure upstream of bed 
= pressure downstream from bed 
= pressure drop across quartz chips (const. ) 
= pressure in the middle of the bed 


(Pi 


P 2 + 


AP 3> - P 2 
2 


= 427 torr 
= 380 torr 

= 17 torr (predetermined value) 
. 380 + < 427 - 17 2 > - 38 . ° - 
= 395 torr = = 0.52 atm 
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A. 5 SAMPLE CALCULATION OF THE ADSORPTION POTENTIAL, A 
Typical calculation of the A-value for propylene at 15 ppm and room temperature 


where 



P° 

P 

A 


vapor pressure of the pure contaminant at T°K = 11. 2 atm @ 298°K 

partial pressure of the contaminant in the system = C , x P 

inlet system 

298°K 11. 2 atm 

66.6 ml liq/mol g 10 , 6. n __ , 

(15/10 ) x 0.52 atm 

4 - 47 Sfliq log 1 ' 45 X 1()6 


27. 52 


mol °K 
ml liq 
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Appendix B 

SAMPLE WEIGHT AND PRESSURE RECORDS FROM 
GRAVIMETRIC EXPERIMENTS 


A sample of a weight record from typical gravimetric experiments, M20-65 and 66, 
is presented in Fig. B-l. This to a portion of a strip-chart recording of the electrical 
output of the electrobalance unit referred to in section 2. 1. 2. On it are included 
electrobalance settings, temperature, and other pertinent notes concerning the nature 

of the experiment. 

A sample total and partial pressure record for the same experiments is giver , i in Fig. 
B-2, which is the recording from a two-pen strip-chart recorder attached to both the 
total pressure (Barocel) gage and the partial pressure (QKGA) gage described in 
section 3 1 2. Total pressure may be read directly from the recording as indicate , 
while partial pressure of each contaminant may be calculated from the mass spectrom 
eter peaks on the chart. Normally, this was only estimated to ensure that the more 
strongly adsorbed contaminant was not being displaced by the second, less strongly 
adsorbed one. The calculation is similar to that described in section 5. 2. 2, excep 
that the peak M/e = 91 was used for toluene, and the ratio S 43 /S 91 was taken o 
0.9. Partial pressure of toluene would then be given by the relation 


where 

Pt 

h 91’ h 43 
P 


toluene partial pressure (torr) 

peak heights for M/e = 91 and 43, respectively ( P A) 
total pressure (Barocel reading) (torr) 
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Fig. B-l Weight Record for Run No. M20-65 
and M20-66 
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Run No. 
F8-1 

Temp. 

T 

. CC) 

j Sorbent 
30 x 40 BD 

Contaminant 

Concen- 

tration 

c eff 

(ppm) 

Partial 

Pressure 

P 

(torr) 

Total 

Pressure 

P 

(atm) 

Adsorption 
1 Potential 

A 

(mol *K/znl liq) 

Amount Adsorbed 

Vacuum 

Desorptio 

Amount 
n Flow 

Desorbed 
ddf 

(ml liq/g) 

Amount 

Otherwise 

Removed 

Purpose of Experiment 

Remarks 

da 

(ml liq/g) 

(ml S^P/g) 

Time 

(min) 

F9-1 

-2 

-a 

-4 

25 

25 

27 

25 

30 x 40 BD 
30 x 40 BD 
30 x 40 BD 
30 x 40 BD 

(1) *-C 4 H l0 

(2) O.H. 

(D -C 4 H 1q 

< 2 > C 3 H 6 
C 3 H 6 

n ' C 4 H 10 

20.4 

.54 

1.08 

13.2 

13.0 

7.9 x 10~ 3 

2.09 x io -4 
4.52 x nf 4 
5.51 x 10" 3 
5.43 x i<f 3 

0. 51 
.51 
.55 
.55 
.55 

16.6 

34.0 

20.4 

27.6 

27.8 

3,4 x 10 -2 
1.67 X io" 5 
8.0 V 10~ 3 
1.03 X 10“ 3 
1.18 x io“ 3 

4 

7.9 X 10° 
5.62 x io" 3 
1.86 x io° 
3.47 X io' 1 
3.97 x 10' 1 

4 305 

5 760 

80 

1. 4 x 10~ 2 
0 

4.5 X IO -2 
0 

5.9 X 10“ 4 

2.0 x io" 2 
3.5 x 10~ 2 
5.9 x 10 -4 

Desorption rate study and block- 
age of C 3 H 6 byn-C 4 H 10 
Desorption rate study and biock- 

of ^3 H 6 ^ low conc - n- ^4 H io 
Adsorption of C 3 Il g and desorption 
rate study 

i — 

i 

(Single desorption point) 


25 

| 30 x 40 BD |a~C 4 H 10 

— 



- 

* 

- 






Fll-l 

-2 

-3 

-A 

-5 

—6 

-7 

1 

24 

25 
25 
27 

27 

25 

25 

30 X 40 BD 
30 x 40 BD 
30 x 40 BD 
30 x 40 BD 

30 X 40 BD 
30 x 40 BD 

30 x 40 BD 

n ~ C 4«10 

D - C 4 H 10 

°~ C 4 H 10 

<i>-c 4 h 10 

< 2 > C 3 H 6 
C 3 H 6 
(1) C_H fi 
<2) tHC 4 H 1Q 

(D-c 4 h 10 

< 2 > C 3 H 6 

21.4 
17.8 

18.5 
0.9 

1580 

1600 

13.2 

20.2 

19.5 
13.4 

9.02 x l(f 3 
7.45 x 10 -3 
7.74 x 10“ 3 
3.8 X 10~ 4 
6.61 X io' 1 
6.33 x io' 1 

5.32 x l(f 3 
8. 13 X 10“ 3 
8.00 x io' 3 
5.50 x io' 3 

0.56 

.55 

.55 

.55 

.55 

.52 

.53 

.53 

.54 

.54 

16.4 

16.7 
16.6 
20.9 

18.5 

18.6 

27.8 
16.55 
16.6 
27.8 

2.98 x 10* 2 
2.41 X 10“ 2 
2.62 X io~ 2 
5.7 x IO -3 
4.0 x io~ 2 
4. 5 x 10' 2 

1.35 X l(f 3 
! 2. 71 x l(f 2 
2.58 x 10~ 2 
4.13 x 10~ 4 

6.91 X 10° 
5.59 x 10° 
6.08 X 10° 
6.32 X 10° 
1.35 X 10 1 
1,51 x io 1 

4.54 x Kf 1 
6.30 x io° 
6.00 x io° 
1.39 x 10' 1 

6 431 
1 050 
16 767 

1.01 X i<f 2 
1.86 X 1Q~ 2 
0.33 X ia“ 2 
5,3 x 10“ 3 
3.8 x io" 2 
4.6 x 10’ 2 

4.5 x io -4 
2.81 x uf 2 

1.97 X io“ 2 
0.55 x ia“ 2 
2.29 X IO -2 

9.0 x io' 4 

Desorption rate study 

Desorption rate study 

Desorption rate study 

Blockage of high come. CgHg by 

low cone. n-C.H,. 

4 10 

Adsorption of CgHg 

Displacement of C-H. bv n-C,H A 
o 6 “ 4 10 

Blockage of OjHg by n-C 4 H 10 

q dfspl.ce<i* 117> ' 10 ' 3mll Vg 

F12-1 

-2 

25 

25 

30 x 40 BD 
30 x 40 BD 

< l > C 3 H 6 
(2) C 2 H 4 

C 2 H 4 

13.8 

27.2 

26.2 i 

5.35 x l(f 3 
1.05 x l(f 2 
9,95 x io" 3 

0.51 

.51 

.50 

27.8 

40.3 

40.3 

1.49 x 10' 3 
4,97 x io" 5 
4.88 x 10 -5 

5.01 x iff 1 
2.25 x 10 -2 
2.22 x io -2 

- 

1.52 x IO -2 
4. 12 x 10‘ 2 
4. 64 x 10~ 5 

- 

Blockage of C 2 H 4 by C^Hg 
Adsorption of C 2 H 4 


F13-1 

-2 

25 

25 

8 x 12 BD 
8 x 12 BD 

- c « H l( > 

r c i H w 

19.5 

25.0 

8.15 x 10“ 3 
1.00 x IO* 2 

0.55 

.53 

16,2 

16.3 

2.63 x 10‘ 2 
3.53 x 10‘ 2 

6. 11 x 10° 
8.20 x 10° 

4 560 

3.3 x 10“ 3 1 

| 

Desorption Rate Study 
Desorption Rate Study 

Bed leak 

F14-1 

25 

8 x 12 BD 

- c ."i. | 

25.3 

1.00 x 10 -2 

0.52 

16.4 

3.28 x 10‘ 2 

7.62 x 10° 

4 320 

2.5 x 10 -2 

3. 6 x 10~ 2 

Desorption Rate Study 

Quartz chips in bed 

F15-1 

-2 

-3 

25 

24 

105 

i X 12 BD 
8 x 12 BD 
8 x 12 BD i 

a - C 4«10 

Q - C 4 H !0 

°- C 4 H 10 

25.3 

26.4 
27.0 

1.01 x l<f 2 

1.05 X 10~ 2 
1.05 x 10* 2 

0.52 

.52 

.51 

16.3 

16.2 

24.2 

3.40 X 10" 2 
3.56 x IO -2 
1.05 x io" 3 

7.90 x 10° 
8.26 x 10° 
2.44 x 10 _l 

5 750 
5 580 
11 250 

2.9 x 10“ 2 
2.9 X l(f 3 

3.5 x 10 -2 

3.6 x 10" 2 
1.2 x IO -3 

Desorption Rate Study 
Desorption Rate Study i 

Desorption Rate Study 

Confirms F13-2 data 

Confirms F13-2 data 
Bed problem 

| F16-1 

25 

8 x 12 BD i 


26.5 

1.03 x 10 -2 

0.51 

16.2 

3.64 x io" 2 

8.45 X 10° 

10 200 

5. 1 x 10“ 3 

3.6 x IO -2 

Desorption Rate Study 

Long bed (9 g) 




Run No. 

Temp. 

T 

(*C) 

Sorbent 

Contaminant 

Concen- 

tration 

c eff 

(ppm) 

Partial 

Pressure 

P 

(torr) 

Total 

Pressure 

P 

(atm) 

Adsorption 
potential ' 

A 

mol "K/ml liq) 

F17-1 

-2 

-3 

107 

107 

107 

8 x 12 BD 
8 x 12 BD 

C 7 H 8 
“' C 4 H L0 
(1) b-C 4 H 10 
< 2 > C^f, 

27.3 

27.0 

26.9 

26.8 

1.05 x 10‘ 2 
1.04 x 10' 2 
1.04 x 10" 2 
1.04 x 10" 2 

0.51 

.51 

.51 

.51 

15.0 

24.2 

24.8 

15.4 

-4 

67 

8 x 12 BD 

(1) CjH, 

B«-c 4 » 10 

(3) CC1 2 F 2 

26.5 
46.0 

38.5 

1.01 X 10 
1.75 x IQ' 2 
1.46 x 10“ 2 

.50 

.50 

.50 

12.2 

19.5 

25.0 

FI 8-1 

28 

30 x 40 BD 

(1) -c 4 h 10 

(2) H O 

26.4 

1.16 X 10~ 2 

0.58 

16.1 





2.12 x 10 1 


1. 19 

FI 9-1 

28 

30 x 40 BD 

iM. 

(1)H 2 o 


2.12 x 10 1 

0.505 

2.01 

-2 

66 

30 x 40 BD 

(2) n-C 4 H 10 

c 2 hoh 

24.0 

60 

9.12 x 10 3 
2.3 x 10“ 2 

.50 

.506 

16.5 

23.5 

-3 

26 

30 x 40 BD 

<w>« 

35.5 

1.34 X 10 2 

.497 

17.5 

-4 

-5 

64 

27 

30 x 40 BD 
1 30 x 40 BD 

(1) n-C 4 H 1Q 

(2) i^O 

(i) c 2 h 5 oh 

25.5 

20 

1.00 x 10 2 
2.12 x 10 1 
7.6 x 10 -2 

.514 

.502 

20.4 

16.7 

19.1 

-6 

66 

30 x 40 BD 

( 2 ) 1^0 
<« C 2 H 5 OH 

28 

2.12 x 10 l 
1. 06 x 10“ 2 

.50 

1.6 

25.2 




(2) I^O 


2.12 x 10 l 
-3 


17.5 

-7 

-8 

27 

65 

30 x 40 BD 
30 X 40 BD 

(1) n-C 4 H lQ 

(2) H 2 0 

(1) n-C 4 H 10 

25.6 

25.8 

9.7 X 10 
1.50 x 10 1 
1.00 x 10" 2 

.50 
. 508 

16. 5 
4.0 
20.1 




(2) H 2 0 


1.50 x 10 1 


19.8 

-9 

27 

30 X 40 BD 

(1) c 2 h 5 °h 

(2) i^O 

27.0 

1.03 x 10' 2 
1.50 x 10' 1 

.501 

18.6 

4.0 

-10 

66 

30 x 40 BD 

w C 2 H 5 OH 

16.2 

6.2 x 10 * 

.505 

26.4 




(2) 1^0 


1.50 x 10 1 
-3 


20. 1 

-11 

26 

30 x 40 BD 

i (1) n-C 4 H l0 

25.1 

9.4 x 10 d 

.49 

16.5 




(2) H 2 0 


1.00 x 10 1 


6.4 

-12 

67 

30 x 40 BE 

1 (1) "“ C 4 H 10 

26.2 

1.00 x 10" 2 
1 

.50 

20.3 






1.00 X 10 

-3 


23.7 

-13 

26 

30 x 40 BI 

) (1) c 2 h 5 oh 

24 

9.1 x 10 3 

.50 

18. 5 




(2)^0 


1.00 x 10 1 
-3 


6.4 

-14 

66 

30 x 40 BI 

) (1) C 2 H 5 OH 

20.5 

7.8 x 10 

1 .50 

! 25.9 




(2) H z O 

. 

1.00 x 10 1 


23.3 


Amount Adsorbed 


q a <u 

(ml liq/g) (ml STP/g) 


Vacuum Amount 

Desorption Flow 

Time Desorbed 

l d q df / 

(min) 

2 800 1.8 X io“ 2 

1 000 o 

5 250 o 

5 250 l.X X 10~ 2 

10 250 5.2 X 10" 2 

10 250 0 


Amount 

Otherwise 

Removed 

q* - qdf 

(ml liq/g) I 


Purpose of Experiment 


3.1 x 10 - ' 
5. 8 x 10“' 

5.2 X 10“ ! 

9.0 x 10“' 

2.1 x 10“ 

9. 5 x 10" 


Desorption Rate Study 
Desorption Rate Study 
Desorption Rate Study 
Desorption Rate Study 
Desorption Rate Study 
Desorption Rate Study 
Desorption Rate Study 


Displacement and vacuum 
desorption of two contaminants 
Blockage and vacuum desorp- 
tion of three contaminants 





TnlH»l 

Weight 

(Degassed) 


Partial 

Contain- Pressure 


C 3 H 8 

50. 318 C 3 H 8 

49.477 I 

49.091 C3H3 

50.000 None 


CgH5CH3 7.0 

C6H5CH3 
C 6 HsCH 3 3. 5 

C 6 HsCH3 
c 6 h 5 ch 3 
C 6 H 5 CH3 3 - s 

c 6 h 5 ch 3 - 

Q5K5CH3 15 


Adsorption . . 

Weight Adsorption Amount 

Increase Potential Adsorbed 

aw* ^ Qa 

(mol*K/ml Uq) (ml liq/g) 


f 300 

n~C4H 10 6. 7 


3.16 x 10-1 
2. 12 x 10-I 
2. 23 x lO'l 
2. 25 x io _1 
2.30 x 10” 1 


9.9 x 10-2 
3, 44 x 10-2 
2.82 x 10-Z 
1. 51 x 10-2 
5.38 x 10-2 

1.78 x 10-2 
3.25 X 10-2 

3. 89 X 10- 3 
3. 05 x 10~2 
3. 66 X 10" 3 
4.74 x 10-3 


3. 90 x 10 -2 
1. 21 x 10-2 

3.31 X 10-2 
3.34 X 10-2 
3.33 x io-2 
3.30 x 10-2 


Desorption Details 


Al uminum 


Flow nm in O 2 at 0. 5 atm 

Flow run in O 2 /N 2 at 0. 5 atm 

Differential flow runs in O 2 /N 2 at 
0. 5 atm 

Vacuum background wt_ increase 
Instrument drift 

System Check 

Propane adsorption and desorption 


1 Nickel screen pans on all Minimize flow effects on desorption d 
experiments henceforth 


Check on magnitude of pressure 
effects of unadsorbed gas with 
similar molecular weight 


Wt. variation with pressure; buoyancy 


Excessive base line drift 
Excessive base line drift 

~0. 3%/hr wt. increase 

+5pg/hr 

Leaky 


Desorption data useless due 
to Strong flow effects 


I Av - 3. 22 x 10 -3 mi liq/g 
Desorption curve appears 
normal after few seconds 


Av qa = 2. 66 x 1 0-3 ml liq/g 


Leaky manifold - data useless 

Leaky manifold - data useless 
System leak 

Microbalance slit out of adj. 


>50% desorbed 
Manifold leak — data useless 
Data not used 
Data not used 


Adjust furnaces for minimum therm ogravimetrtc effects 

[ (Av AW a and 


M13-30D = long desorption 
M13-57D j long desorption 


M14-21D ■ long desorption 


M14-61D* long desorption 
M14-83D = long desorption 


High chart speed 
High chart speed 
High chart speed 
High chart speed 


1 27 *C isotherm 

Find T such that q at 200p = q at 20p A 27*C 
62 *C isotherm 

C 3 H 6 adsorption on "dirty" charcoal 

C 3 Hg adsorption on fresh charcoal 
47 'C isotherm 
Repeat M13-58 to -60 

Repeat M14-1 to -3 after pump oil removal 
Brief 27 *C Isotherm after pump oil removal 
Brief 62 *C isotherm after pump oil removal 
37 ‘C isotherm 
77 'C isotherm 

System check - 00 data takrm 

Long desorption after isosteric 
adsorption: variation with 
temperature 


(q&A at p = 20j*Hg) 

T too high 
(q&A atp - 20pHg) 

No significant difference 
No significant difference 
(q&A atp - 20pHg) 

No significant difference 
No significant difference 
No significant difference 
No significant difference 
(q i A at p » 20pHg) 

(q k A at p = 20pHg)(sample loss) 

(sample loss) 

(sample loss) 


Appendix D 

GRAVIMETRIC SYSTEM EXPERIMENTAL SUMMARY 




-10 

-11 

-12 

-13 

-14 

-15 

-16 

-17 

-18 

-19 

-20 

-21 

-22 to -24 
-25 
-26. -27 
-28 to -34 
-35 to -37 


Temp. 

T 

(*C) 

Sorbent 

Initial 

Weight 

(Degassed) 

wo 

<g> 

Contam- 

inant 

Partial 

Pressure 

Pa 

(pHg) 

Adsorption 

Weight 

Increase 

AW a 

(mg) 

Adsorption 

Potential 

A 

mol*K/ml liq) 

Amount 

Adsorbed 

(ml^q/g) 

26.6 

8 x 12 BD 

50.365 

n-CdHin 

13.2 

1.009 

16.1 

3.33 x 10-2 

26.2 

| 


50.347 



14.7 

1.010 

15.9 

3.33 x 10-2 

86.8 

1 


50.331 



143 

1.005 

17.9 

3.31 X 10-2 

67.0 

1 


50.285 



102 

1.002 

16.7 

3.31 xiO-2 

46.8 

♦ 


50.300 

n-C4Hio 

48 

1.015 

16.3 

3.35 x 10-2 

26.4 

8xl 

2 BD 

50. 308 

c 3 h 8 


9.6 

.102 

23.7 

3.43 x 10-2 

27.4 

30 x 40 BD 

49.750 

C 3 H8 


9.3 

0.109 

23.8 

3.76 x 10' 3 

27.4 

30 x 40 BD 

49.759 

C 3 H 8 


8.3 

.099 

24.1 

3.36 x lO" 3 

27.3 

50 x 70 BD 

' 49.970 

n-C4Hio 

7.6 

0.997 

16.9 

3.31 x 10-2 

87.6 

16 x 20 BD 

50. 191 

0-C4H1O 

220 

1.004 

17.2 

3.32 x lO -2 

87.3 



50.182 



235 

1.001 

17.1 

3.31 x 10-2 

87.3 



50. 185 



220 

1.000 

17.2 

3.31 x 10-2 

87.7 



50.219 



9.8 

.130 

22.3 

4.30 x lO" 3 

87.1 



50. 233 



9.5 

. 133 

22.3 

4.39 X lO- 3 

87.1 



50.248 



7.5 

.130 | 

22.6 

4. 30 x lO- 3 

86.9 



50.245 



>430 

2. 598 


8. 59 x 10"2 

87.2 



50.233 



22 

.264 

20.9 

8.74 x lO -3 

87.0 



50. 241 



25 

.265 

20.8 

8.78 X lO" 3 

86.3 



50. 232 



25 

.264 

20.8 

8,74 x lO" 3 


-38 

66.8 

-39 to -41 

66.8 

-42 to -47 

66.8 

-48 

.8 

-49 to -52 

106.8 

-53 to -59 

106.8 

■ -60 to -61 

106.6 

-62 to -68 

66.7 

-69 to -75 

27.0 

-76 to -78 

26.8 

-79 A/D 

26.8 

-80 to -82 

26.8 

-83 to -86 

106.8 

-87 to -90 

106.8 

-91 to -95 

106.9 

M19-1 to -A 

26.8 

-5 to -7 

26.8 

-8 to -9 

106.8 

-10 to -15 

107 

-16 to -20 

106.6 

-21, -22 

106.6 

-23 

107 

-24 

27 

-25 

27 

-26 

26.8 

-27 

66.8 

-28 

107 

-29 

107 

-30 

26.8 


16 * 20 BD 
30 x 40 BD 


n-C 4 H l( ) 

C6H5CH3 


C6H5CH3 0.1 to 1. 
11-C4H10 1 to 2. 

Q-C4H10 30/8.4 

n-C4Hio 7.1 to 10 
C 6 H 5 CH 3 0. 26 to 2 
n-C4Hio 1.35 to 1 
o-C 4 Hio 0. 1 to 30 

0-C4H10 0.8 to 30 

0-C4H10 16 to 100 

n-C4Hio 1-0, 3.0 
11-C4H10 1 to 1000 

C6H5CH3 0.2 to 10 
C 6 H 5 CH 3 48 to 100 

11-C4H10 20 

C 6 H 5 CH 3 40 total 
ln-C4Hio 20 
lc 6 H 5 CH 3 40 total 
C6H5CH3 -0.1 
\ Q -C 4 H T o 48 
n-C4Hio 9.0 
n-C4Hio 98 

n-C4Hio 710 
h-C^Hiq 8 . 45 


1 to 100 1 

- Various 

- 

200 to 1000 

Various 

- 

100 | 

9.4 

2.69 x lO" 1 

200 to 1000 

- Various 


1 to 50 

- Various 

— 

91 

10.24 9.5 

2.64 x 10" 1 

1 to 10 

— Various 

- 

10 to 100( 

— Various 

— 

.070 | 

| 10.23 9.0 

2.64 x 10-1 

0,08 to 1.1 

4 - Various 

- 

0.1 to 1.46 

- 

_ 

1 to 2. 2 I 


— 

30/8.4 | 


- 

7.1 to 1000 

1 — ! 

— 

0 . 26 to 20 

1 

— 

1.35 to 10C 

to 

_ 

0. 1 to 300 

- Various 


0.8 to 30 

Various 

- 

16 to 1000 

- Various 

- 

1.0, 3.0 

- Various 

— 

1 to 1000 

- Various 

- 

0.2 to 100 

— 

— 

48 to 1000 

_ _ _ _ 

- 


20 

3.098 total 
1.447 
11. 69 total 
1.900 
3. 112total 
.976 
.971 


Desorption Details 


High chart speed'] 


High chart speed J 

High chart speed 
High chart speed 

High chart speed 
High chart speed 


3.33 x 10-2 

3.30 x 10'2 

3.31 x 10-2 
3.26 x 10-2 

3.31 X 10-2 

3.32 x 10-2 
3.32 x 10-2 
1.78 x 10-1 


. High chart speed 
High chart speed 


3.90 x 10- 3 Mass spec analysis 
4. 97 x 10-2 Mass spec analysis 

5, 15 x 10~ 2 Mass spec analysis 

3.43 x 10-2 

3.41 x 10-2 

3.40 x 10-2 

3.42 X 10-2 


Long desorption after isosteric adsorption: 
variation with temperature 
Also variation with particle size from 
previous group of experiments 

Long desorption after isosteric adsorption: 
variation with particle size 

Variation with particle size 

Long desorption after isosteric adsorption: 
reproducibility, variation with temperature, 
and variation with particle size 


Huge the rmogravi metric effects 


Isotherm points 
Long desorption 

Calibration of Pi rani gage after adj 
Isotherm points 


Isotherm points 
Long desorption 
Isotherm points 
Isotherm points 
Long desorption 
Isotherm points 


Isotherm points 
27 “C isotherm points 
27* C isotherm points 
107*0 isotherm points 
107 °C isotherm points 
107 *C isotherm points 
107 °C isotherm points 

2-component desorption 
2 -component desorption 

2-component desorption 
Long desorption, repeat M13 -90 
Long desorption, repeat M13-88 

Long desorption 

Long desorption, repeat M19-26 


Data unsatisfactory - huge AW/At 


Data unsatisfactory 


Aborted - pan stuck to tube 

First expt w/ oon-magnetic 
Long furnaces 


Used later data instead 


Not done rigorously for isotherm use 
First Baroce) isotherm data 



212 


Temp. 

T 

rc> 


TniHal 

Weight 

(Degassed; 

W° 

(g> 

Contain- 

in±nt 

Partial 

Pressure 

Pa 

0*Hg) 

T 

Adsorption 

Weight 

Increase 

AW a 

(mg) 

Adsorption 

Potential 

A 

(morK/ml Uq) 

| Amount 
Adsorbed 
qa 

[ (ml Uq/g) 

50.618 

fCeHsCHi 

0.8 

6. 588 
1.477 



ln-C4Hio 

380 


4.84 x 10-3 

50.758 

cci 2 f 2 

I 

0.1 to 1000 

— 

Various 


50. 735 

i.l to 30 

— 

i 

_ 

50.752 

1 

100 to 725 

_ 

f 

_ 

50.760 

T 

50 to 1000 

- 

T 

„ 

50.807 

cci 2 f 2 

1000. 21 

- 

Various 

_ 

50. 819 

C6H5CH3 

<0.1 

1.043 

>14 

2.64 x io-2 

50.819 

n-C4Hio 

1.5, 5.8 

- 

Various 

_ 

50.710 

C6H5CH3 

<0.1 

1.044 

>17 

2.64 x io-2 

50.710 

n-C4Hio 

0. 9 to 983 

- 

Various 

_ 

50.709 

c 6 h 5 ch 10 

-.08 

1.043 

-18 

2.64 x io-2 

50. 709 

n-C4Hio 

2.9 to 1000 

- 

Various 

_ 

50. 707 

n-C^Hjo 

2.6 to 453 

- 

Various 

_ 

50. 707 

c 6 h 5 ch 3 

6.8 

10.432 

9.3 

2.64 x 10-1 

50. 732 

c 6 h 5 ch 3 

7.2 

10.416 

9.3 

2.64 x 10-1 

50. 732 

n-C 4 H 10 

22. 2 to 997 

- 

Various 

— 

50.703 

c 6 h 5 ch 3 

91 

10.424 

9.5 

2.64 x 10“1 

50.703 

n-C 4 Hio 

10.5 to 914 

— 

Various 


50,708 

n-C 4 Hio 

0.73 to 980 

- 

Various 


50.774 

C6H5CH10 

0.55 

2.087 

16.2 

5.28 x 10-2 

50.774 

50.768 

n “ C 4 H 10 

c 6 H 5 ar 10 

3.75 to 992 
<.01 

2.089 

Various 

16.0 

5. 28 x 10-2 

50.768 

n-C 4 Hio 

1.8 to 985 

— 

Various 


50. 707 

n-C4Hio 

0.64, 1.20 

— 

Various 

_ 

50.739 

□-C4H10 

5. 9 to 2xl 0* 

- 

Various 

_ 

50.718 

n-C^io 

199 to 2*10* 


Various 


50.014 

n-C4Hio 

10 

-0.002 



50. 014 

I 

LOO 

- .001 



50. 014 

i 

1 000 

.000 



50.014 

t 

10 000 

+ .004 



50. 012 

C 2 H 5 OH 

10 

- .003 



50.012 

100 

- .001 



50,012 

I 

1 000 

- .001 



50.012 

t 

10 000 

.000 



“ 

C6H5CH3 





— 

n-C^io 






C6H5CHJ 






n-C4Hio 






C 2 HsOH 





~ 

C2H50H 






C2H5OH 





50.862 

H20 

5 000 

0. 853 

9.34 

1.75 x lO- 2 

50.862 

1 

10 000 

4.121 

4.60 

8.46 xlO-2 

50.70 


— 



50.70 

i 

15 000 

21.29 

1.84 

4.37 x 10-1 

50. 70 

T 

20 000 

27.82 

0.05 

5. 71 x 10-1 

51.05 

h 2 o j 

5 000 

0.58 

13.8 

1.19X10-1 

51.05 

1 

10 000 

1.17 

8.90 

2. 39 x 10-2 

51.05 

i 

15 000 

2. 63 

6. 05 

5. 38 x 10-2 

51.05 

♦ 

20 000 

9. 16 

4.02 

1.87 x 10-1 

50.993 

h 2 o 

5 000 

0.259 

29.1 

5.31 x 10-3 

50.993 

1 

10 000 

.345 

23.7 

7. 02 x 10-3 

50.993 


15 000 

.427 

20.5 

8.76 x 10-3 

50.993 

i 

20 000 

.497 

18.2 

1.02 x 10-2 

50.993 

♦ 

25 000 

.595 

16.5 

1.22 x 10-2 

51.10 

H z O 

5 000 

0.36 

21.7 

7.3 x 10-3 

51.10 


10 000 

.535 

16.2 

1. 10 x 10-2 

51.10 


15 000 

.725 

13.1 

1.48 x io-2 

51.10 

51.10 


20 000 

1.025 

10.8 

2.10 xl0~ 2 


Desorption Details 


Purpose 


He marks 


M20-1 

-2 to -8 
-9 to -12 
-13 to -15 
-16. -17 
-18 A/D 
-19 

-20. -21 

-22 

-23 to -30 
-32 

-32 to -38 
-39 to -41 
-42 
-43 

-44 to -49 
-50 
-51 to -56 
-57 to -84 
-65 
-66 to -73 
-74 
-75 to -81 
-82. -83 
-84 to -90 
-91 to -95 


27.1 

66.9 

26.8 

26.8 

26.8 

66.8 

26.8 

26.7 

26.8 
26.8 


27.0 

26.8 


66.8 

66.8 


26.8 

26.8 

26.8 


8 x 12 BD 


8 x 12 BD 


2 -component desorption 

67 "C Freon-12 isotherm points 
27 “C Freon-12 isotherm points 
27 *C Freon-12 isotherm points 
27 *C Freon-12 isotherm points 
67 "C Freon-12 isotherm points 
Toluene praadsorption - blockage 
of n -butane adsorption isotherm 
Toluene praadsorption - blockage 
of n- butane adsorption isotherm 
Toluene preadsorption - blockage 
of n -butane adsorption isotherm 
Few 27 *C isotherm points 
Toluene preadsorption 
Toluene preadsorption — blockage 
of n -butane adsorption isotherm 
Toluene preadsorption — blockage 
of n-butane adsorption isotherm 
67 “C n-butane isotherm points 
Toluene preadsorption - blockage 
of n-butane adsorption isotherm 
Toluene preadsorption - blockage 
of n-butane adsorption isotherm 
27 *C ti - butane isotherm points 
27 ‘C n-butane isotherm points 
67 ®C n-butane isotherm points 


Poor reproducibility 


Poor reproducibility 
Mass spec analysis shows 
Freon-12 contamination 
Mass spec analysis shows 
no displacement of toluene 
Mass spec analysis shows 
no displacement of toluene 

Aborted - mass spec failure 
(mass spec not used) 

Mass spec analysis shows 
large displacement of toluene 

Mass spec analysis shows 
no displacement of toluene 
Mass spec analysis shows 
negligible displacement of toluene 


M21-1 

-2 

-3 

-4 

-5 

-6 

-7 

-8 

-9 

-10 

-11 

-12 

-13 

-14 

-15 


"Blank" run in n-butane 

"Blank" run in n-butane 

"Blank” run in n-butane 

"Blank" run in n-butane 

"Blank" run in ethanol 

"Blank" run in ethanol 

"Blank" run in ethanol 

"Blank" run in ethanol 

Maas spec ionization cross-section 

determinations 


Negligible weight change 


Negligible weight change 


M22-1 
-2 
-3. -4 
-5 
-6 
-7 
-8 
-9 
-10 
-11 
-12 
-13 
-14 
-15 
-16 
-17 
-18 
-19 


31.8 


30 x40 BD 


30 x 40 BD 


22 ‘C water isotherm point 


22 *C water 
32 "C water 
32 °C water 
32 *C water 
32 "C water 
67 -C water 


isotherm point 
isotherm point 
isotherm point 
isotherm point 
isotherm point 
isotherm point 


67 °C water 
47 *0 water 
47 “C water 
47 *C water 
47 *C water 


isotherm point 
isotherm point 
isotherm point 
isotherm point 
isotherm point 


Aborted — ins uffi cient tim e 


Aborted - pressure won't go above 
20 torr 




Appendix E 

SUMMARY OF ISOTHERM POINTS FROM DESORPTION RUNS ON 

GRAVIMETRIC SYSTEM 


Table E-l 

SUMMARY OF ISOTHERM POINTS FROM DESORPTION RUNS ON GRAVIMETRIC 

SYSTEM - n-BUTANE 


Run No. 

Particle 

Diameter 

(mm) fC) 

Temperature 

Equilibrium 

Pressure 

Amount Adsorbed 

Adsorption 
Potential 
(mol °K/ml liq) 

(°K) 

Pirani 

(torr) 

Actual 

(torr) 

(ml STP/g) 

(ml liq/g) 

M14-87 

0.5 

46 

319 

0.102 

0. 128^ 

7.70 

3.37 x 10' 2 

16.8 

M14-88 

0.5 

67 

340 

. 297 

.093 

7.76 

3.34 x 10 

16.9 

M14-89 

0.5 

87 

360 

.585 

.30 

7.75 

3.33 x 10 

16.8 

M14-90 

0.5 

26 

299 

.019 

. 0067^ 

7.68 

3.31 x 10" 2 

16.9 

Ml 5-1 

2.0 

27 

300 

0.0345 

0.0132^* 

7.73 

3,33 x 10" 2 

16.1 

Ml 5-2 

2.0 

26 

299 

.0385 

.0147 

7.74 

3.33 x 10 

15.9 

Ml 5-3 

2.0 

87 

360 

. 550 

.143 

7.70 

3.31 x 10 

17.9 

M15-4 

2.0 

67 

340 

.385 

,102 

7.69 

3,31 x 10 

16.7 

Ml 5- 5 

2.0 

47 

320 

.165 

.048 

7.78 

3.34 x 10“ 2 

16.3 

M17-1 

0.25 

27 

300 

0.0225 

0.0076^ 

7.70 

3,31 x 10~ 2 

16.9 

M18-1 

1.0 

88 

361 

0.90 

0.220* C * 

7.71 

3.32 x 10’ 2 

17.2 

M18-2 

1.0 

87 

360 

.95 

.235 

7.70 

3.31 x 10 

2 

17. 1 

M18-4 

1.0 

87 

360 

.90 

.220 

7.69 

3.30 x 10 

_Q 

17.2 

M18-5 

1.0 

88 

361 

.026 

.0098 

1.00 

4.30 x 10 

-3 

22.3 

M18-6 

1.0 

87 

360 

.025 

.0095 

1.02 

4.39 x 10 

22.3 

Ml 8-8 

1,0 

87 

360 

.020 

.0075 

1.00 

4.30 x 10 

22,6 

M18-9 

1.0 

87 

360 

2 


19.97 

-3 


M18-10 

1.0 

87 

360 

.085 

.022 

2.03 

8.72 x 10 

_Q 

20.9 

M18-11 

1.0 

87 

360 

.095 

.025 

2,04 

8.75 x 10 

20.8 

M18-12 

1.0 

87 

360 

.095 

.025 

2.03 

8.74 x 10 

-2 

20.8 

M18-14 

1.0 

26 

299 

.0445 

.0171 

7.74 

3.33 x 10 

_2 

15.7 

M18-15 

1.0 

27 

300 

.0255 

.0096 

7.68 

3.30 x 10 

_o 

16,5 

M18-16 

1.0 

87 

360 

.70 

.175 

7.70 

3.31 x 10 

-2 

17.6 

M18-18 

1.0 

88 

361 

.65 

.165 

7.70 

3.31 x 10 

-2 

17,7 

M18-19 

1.0 

88 

361 

.90 

.225 

7.71 

3.32 x 10 

_ 9 

17.2 

M18-20 

1.0 

88 

361 

.80 

.20 

7.73 

3.33 x 10 

17.4 

M19-23 

0.5 

107 

380 

- 

0.020 (a) (b) (c) (d) 

0.905 

3.90 x 10” 3 

_2 

23.0 

M19-24 

0.5 

27 

300 

- 

.020 

11.53 

4.97 x 10 

-2 

15.5 

M19-26 

0.5 

27 

300 

- 

.0090 

7.96 

3.42 x 10 

-2 

16.6 

M19-27 

0.5 

67 

340 

- 

.098 

7.92 

3.41 x 10 

-2 

16. 8 

M19-29 

0.5 

107 

380 

- 

.710 

7.89 

3.39 x 10 

_2 

16.8 

M19-30 

0.5 

27 

300 

- 

. 00845 

7.94 

3.41 x 10 

16.7 


(a) McLeod gage reading. 

(b) From calibration versus McLeod gage. 

(c) From calibration versus Barocel. 

(d) Barocel reading. 
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Table E-2 

SUMMARY OF ISOTHERM POINTS FROM DESORPTION RUNS ON GRAVIMETRIC SYSTEM - TOLUENE 


Run No. 

Particle 

Diameter 

(mm) 

Temp< 

srature 

Equilibrium 

Pressure 

Amount Adsorbed 

Adsorption 

(°C) 

(°K) 

Pirani 

(torr) 

Actual 

(torr) 

(ml STP/g) 

1 ~ 

(ml liq/g) 

Potential 
(mol °K/ml liq) 

M9-1 

0.5 

27 

300 

0.0020 

0.0074* a) 

51.6 

2.72 x 10 -1 

9.2 

Ml 0-2 

0.5 

27 

300 

0.0010 

0.0035^ 

47.2 

2.50 x 10 _1 

9.6 

Ml 3-1 

0.5 

87 

360 

0.0010 

0.0035 (a > 

18.7 

-2 

9.9 x 10 

14.8 

M13-2 

0.5 

i 87 

360 

.0010 

. 0035 (a) 

19.0 

1.00 x 10 -1 

14.8 

M13-3 

0.5 

87 

360 

.0010 

.0035* a) 

18.7 

9.9 x 10 -2 

14.8 

M18-21 

1.0 

27 

300 

0.0028 

o.ooio (a) 

33.8 

1.78 x 10 -1 

11.4 

Ml 8-25 

1.0 

27 

300 

.014 

.0050 

50.1 

2.65 X 10 -1 

9.7 

M18-48 

1.0 

67 

340 


.091^ 

50.1 

2.65 X 10 -1 

9.5 

M18-60 

1.0 

107 

380 

- 

1.070 

50.1 

2.65 X 10 -1 

9.0 

M19-25 

0.5 

27 

300 

— 

0.0001^ 

0.75 

5.15 X 10" 2 

14 

M2 0-1 

2.0 

27 

300 


0.0008^ 

31.6 

1.67 x 10 -1 

11.7 

M2 0-19 

2.0 

27 

300 


.0001 

5.00 

2.62 x 10" 2 

14 

M20-22 

2.0 

27 

300 

— 

.00001 

5.01 

2.64 x 10“ 2 

17 

M2 0-31 

2.0 

67 

340 

- 

.00008 

5.01 

2.64 x 10' 2 

18 

M2 0-42 

2.0 

27 

300 

— 

.0068 

50.1 

2.64 x 10 -1 

9.3 

M2 0-43 

2.0 

27 

300 

- 

.0072 

50.0 

2.63 x io -1 

9.3 

M20-50 

2.0 

67 

340 

— 

.091 

50.1 

2.64 x io -1 

9.5 

M2 0-65 

2.0 

67 

340 

— 

. 00055 

10.0 

5.28 x 10~ 2 

16.2 | 

M20-74 

2.0 | 

27 

300 

— 

.0001 

10.0 

5.28 x 10 -2 

16.0 


(a) From calibration versus Barocel. 

(b) Barcel reading. 
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Table E-3 

SUMMARY OF ISOTHERM POINTS FROM DESORPTION RUNS ON GRAVIMETRIC SYSTEM - PROPANE 


Run 

No. 

Particle 

Diameter 

(mm) 

Temperature 

Equilibrium 

Pressure 

Amount Adsorbed 

Adsorption 
Potential 
(mol°K/ml liq) 

(°C) 

(“K) 

Pirani 

(torr) 

Actual 

(torr) 

(ml STP/g) 

(ml liq/g) 

M2-1 

0.5 

27 

300 

— 

0.0133( a ) 

0.69 

2.30 x 10~ 3 

23.1 

M3-1 -2 

0.5 

27 

300 


0.0133( a ) 

1.24 

4.12 x 10' 3 

23.1 

M7-1 

0.5 

27 

300 

0.0133 

0.0055( b ) 

0.86 

2.85 x 10 -3 

24.7 

M7-2 

0.5 

27 

300 

.0133 1 

. 0055( b ) 

1.16 

3.87 X 10 -3 

24.7 

M7-3 

0.5 

27 

300 

.0133 

. 0055( b ) 

1.05 

3.48 x 10 -3 

24.7 

M7-4 

0.5 

27 | 

300 

.0133 

. 0055(h) 

.96 

3.18 x 10" 3 

24.7 

M7-5 

0.5 

27 

300 

. 0133 

. 0055(h) 

.93 

3.08 x 10“ 3 

24.7 

M7-6 

0.5 

27 

300 

.0133 

. 0055(h) 

.98 

3.25 x 10" 3 

24.7 

M7-7 

0.5 

27 

300 

. 0133 

.0055(h) 

.85 

2.82 X 10“ 3 

24.7 

1 

M8-1 

0.5 

27 

300 

0.0133 

0.0055(h) 

0.80 

2.67 x 10" 3 

24.7 

M8-2 

0.5 

27 

300 

.0133 

. 0055(h) 

.81 

2. 70 x 10“ 3 

24. 7 

M8-3 

0.5 

27 

300 

.0133 

.0055(h) 

.80 

2.67 x 10 -3 

24.7 

M8-4 

0.5 

27 

300 

.0133 

. 0055(h) 

.82 

2. 73 x 10 -3 

24.7 

M8-5 

0.5 

27 

300 

.0133 

. 0055(h) 

.76 

2.53 x 10“ 3 

24.7 

M15-6 

2.0 

27 

300 

0.023 

0.0096(h) 

1.03 

3.43 x 10 -3 

23.7 

M16-1 

0.5 

27 

300 

0.0225 

0.0093(h) 

1.11 

3.70 x 10 -3 

23.8 

M16-2 

0.5 

28 

301 

.020 

.0083 

1.01 

3.36 x 10 -3 

24.1 


(a) Based on analyzed concentration in premixed gas from Matheson. 

(b) From calibration versus Barocel. 




Appendix F 

ADSORPTION ISOTHERM DATA FROM FLOW EXPERIMENTS 


Hun No. 

Sorbent 

Contaminant 

Temperatun 

T 

fC) 

Partial 
2 Pressure 

P 

(torr) 

Amount Adsorbed 

Adsorption 

% 

(ml STP/g) 

q a 

(ml liq/g) 

A 

(mol °K/ml Hal 

■ 

16 x 20BD 
16 x 20BD 
16 x 20BD 
16 x 20BD 
16 x 20BD 

C 2»4 

C 2 H 4 

C 2 H 4 

C 2 H 4 

C 2 H 4 

24 

24 

24 

24 

24 

2.28 x 10~ 2 
2.16 x 10" 2 
2.21 x 10’ 2 
2.16 x 10“ 2 
1.89 x 10" 2 

3.75 x lo’ 2 
3.55 x io -2 
3,53 x 10“ 2 
3.55 x 10~ 2 
3.73 x lo“ 2 

8.28 x io’ 5 
7.84 x io" 5 
7.80 x 10“ 5 
7.83 x 10' 5 
8.21 x 10~ 5 


F3-1 

16 x 20 I3jc 

C 2 H 4 

24 

2.25 x 10~ 2 

1. 12 x 10 _1 

2.47 x 10~ 4 

kuh 

-2 

16 x 20 13X 

C 2 H 4 

24 

3.06 x lo" 2 

1.13 x io’ 1 

2.48 x 10~ 4 

37.4 

F5-1 

30 x 40BD 

C 3 H 8 

■ 

1.38 x 10” 2 

1.24 x 10° 

4.12 x IO -3 

22.7 

-2 

30 x 40BD 

C 3»8 


1.29 x 10~ 2 

1. 14 x 10° 

3.79 x 10~ 3 

22.9 

-3 

30 x 40BD 

C 3 H 8 


1.35 x 10‘ 2 

1.10 x 10° 

3,67 x io -3 

22. 8 

-4 

30 x 40BD 

C 3 H 8 

1 

1.22 x io“ 2 

1.02 x 10° 

3.41 x io“ 3 

23. 0 

-5 

30 x 40BD 

C 3 H 8 

1 

1.26 x 10“ 2 

1.08 x io° 

3.60 x 10“ 3 

22. 8 

-6 

30 x 40BD 



1.27 x 10“ 2 

1.05 x 10° 

3.49 x io" 3 

hum 

F6-1 

30 x 40BD 


25 

8.45 x io’ 3 

6.6 x io" 1 

1.96 x 10~ 3 

26. 8 

-2 

30 x 40BD 


46 

8.28 x 10 -3 

1.98 x io' 1 

5.9 x 10 -4 

29.8 

F7-2 

30 x 40BD 

CH 3 CC1 3 

102 

1.01 x lo' 2 

6.16 x io" 1 

2.20 x io’ 3 

27.0 

-3 

30 x 40BD 

C 7 H 8 

102 

9.80 X 10~ 3 

9.85 x io 1 

5.2 x IO -2 

15.2 

-4 

30 x 40BD 

CH 3 CC1 3 

102 

9.97 x 10“ 3 

7.45 x io“ 2 

2.66 x io” 4 

26. B 

-5 

30 x 40BD 

CH 3 CC1 3 

101 

8.1 x 10“ 3 

5.3 x 10 _1 

1.89 x 10‘ 3 

27.3 

-7 

30 x 40BD 

C 7 H 8 

102 

2.55 x 10~ 3 

4.06 x 10° 

2.15 x 10~ 2 

16.9 

F9-1 

30 x 40BD 

n ~ C 4 H 10 

25 

7. 9 x 10 -3 

7.88 x 10° 

3.4 x io" 2 

16.6 

-2 

30 x 40BD 

n - C 4 H 10 

25 

4.52 x lo” 4 

1.86 x xo° 

8.0 x 10~ 3 

20.4 

-3 

30 x 40BD 

C 3 H 6 

27 

5.43 x 10" 3 

3.97 x 10* 1 

1.18 x io" 3 

27.8 

Fll-1 

30 x 40BD 

n - C 4 H 10 

24 

9.02 x 10" 3 

6.91 x 10° 

2.98 x 10~ 2 

16.4 

-2 

30 x 40BD. 

n ' C 4 H 10 

25 

7.45 x 10~ 3 

5.59 x io° 

2,41 x 10 -2 

16.7 

-3 

30 x 40BD 

n - C 4 H 10 

25 

7.74 x io” 3 

6.08 x io° 

2.62 x 10‘ 2 

16, 6 

-4 

30 x 40BD 

n - C 4 H 10 

27 

3.76 x 10 -4 

1.32 x 10° 

5.7 x 10“ 3 

20.9 

>5 

30 x 40BD 

C 3 H 6 

25 

6,33 x 10~ X 

1.51 x 10 1 

4.5 x io” 2 

18.6 

-6 

30 X 40BD 

C 3 H 6 

25 

5.32 x 10“ 3 

4.54 x io* 1 

1.35 x 10~ 3 

27.8 

-7 

30 x 40BD 

"- C 4 H 10 

25 

8.00 x 10~ 3 

6.00 x io 1 

2.58 x 10‘ 2 

16.6 

FI 2-1 

30 x 40BD 

C 3 H 6 

25 

5.35 x 10" 3 

5.01 x io -1 

1.49 x io” 3 

27. 8 

-2 

30 x 40BD 

C 2 H 4 

25 

9.95 x 10" 3 

2.22 x 10~ 2 

4.88 x 10~ 5 

40.3 

F13-1 

8 x 12BD 

"- C 4 H 10 

25 

8, 15 x 10“ 3 

6.11 x 10° 

2.63 x 10" 2 

16. 2 

-2 

8 x 12BD 

n ' C 4 H 10 

25 

1.0 x 10" 2 

8.20 x 10° 

3.53 x 10" 2 

16.3 

F14-1 

8 x 12BD 

n ' C 4 H 10 

25 

1.0 x 10~ Z 

7.62 x 10° 

3.28 x 10~ 2 

16.4 

F15-1 

8 X-12BD 

n - C 4»10 

25 

1.01 x 10" 2 

7.90 x 10° 

3.40 x 10“ 2 

■MB 

-2 

8 x 12BD 

n - C 4 H 10 

25 

1.05 x lo“ 2 

8.26 x 10° 

3.66 x 10‘ 2 

■pvB 

*3 

8 x 12BD 

n - C 4 H 10 

105 

1.05 x 10“ 2 

2.44 x l(f 1 

1.05 x 10‘ 3 


F16-1 

8 x 12BD 

n_C 4 H 10 

25 

1.03 x 10“ 2 

8.45 x 10° 

3.64 x i(f 2 


F17-1 

8 x 12BD 

C 7»8 

107 

1.05 x 10~ 2 

1.10 x IO -1 

5.8 x 10“ 2 


-2 

8 x 12BD 

n - C 4 H !0 

107 

1.04 x 10" 2 ! 

5.04 x IO" 1 

2.17 x 10‘ 3 

24.2 

-4 

8 * 12BD 

C 7»8 

67 

1.05 x io‘ 2 | 

3.68 x IQ’ 1 

1.94 x io" 1 

12,2 
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Appendix G 


DATA FOR HIGH-PRESSURE VOLUMETRIC ADSORPTION EXPERIMENTS WITH 
n-BUTANE ON 30 x 40 BD ACTIVATED CARBON 


Table G-l 

HIGH-PRESSURE VOLUMETRIC-SYSTEM n-BUTANE 
ISOTHERM DATA FOR 30 x 40 BD ACTIVATED CARBON 
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Table G-2 

SELECTED POTENTIAL- PLOT 


Run 

No. 

Temperature 

Equilibrium 

Pr^ccn ro 

Adsorption 
Potential 
(mol°K/ml liq) 

Amount 
Adsorbed 
(ml liq/g) 

(°C) 

(°K) 

(atm) 

Vl-7 

V2-1 

V2-2 

0.0 

273.2 

V 

1 

1.82 

8.1 

7.03 

4. 70 x 10~1 
1.40 x 10 -1 
2.74 x 10-1 

V3-8 

V4-1 

V4-2 

21.5 

294. 7 

1 

I 3.60 

9.4 

7.5 

I 4.29 x 10 -1 
1.48 x 10" 1 
2.62 x 10' 1 

V5-1 

V5-3 

V5-11 

51.0 

324.2 

0.0053 

.0209 

.232 

10.1 

8.04 

4.52 

1.57 x 10"1 
2.35 x 10“1 
3.64 x 10-1 
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Appendix H 

EQUATIONS PERTAINING TO MASS SPECTROMETER SAMPLING TECHNIQUE 

FOR GRAVIMETRIC EXPERIMENTS 


The mass spectrometer section of the gravimetric system pictured in Fig. 2-4 and 
referred to in Section 3. 1. 2 may be represented by the line drawing in Fig. H-l. 

In this figure, the following symbols are used to represent system parameters: 


V 1 

= volume of the microbalance jar section 

w 

V 2 

= volume of mass spectrometer section 

u> 

V 3 

= volume of getter-ion pump section 

w 

Pi 

= partial pressure of a particular species in 
microbalance jar section 

(torr) 

P2 

= partial pressure of a particular species in 
mass spectrometer section 

(torr) 

P3 

= partial pressure of a particular species in 
getter-ion pump section 

(torr) 

Qi 

- throughput of a particular species past the 
limiting conductance between microbalance 
jar and mass spectrometer sections 

(torr-f/sec) 

% 

- throughput of a particular species past the 

limiting conductance between mass spectrometer 
and getter-ion pump sections 

(torr- 9 . /sec) 

f i 

= limiting conductance for a particular species 
between microbalance jar and mass spectrometer 
sections 

(£ /sec) 

F 2 

= limiting conductance for a particular species 
between mass spectrometer and getter-ion pump 
sections 

(£/sec) 

S 

= pumping speed of getter-ion pump for a 
particular species 

(f /sec) 

T 

= response time of mass spectrometer in sampling 
from microbalance jar section 

(sec) 
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Fig. H-l Mass Spectrometer Sampling System for Gravimetric 
Experiments 
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The throughput past F equals that past F„ and is equal to pumping speed times 

1 

pressure in the getter-ion pump: 


Ql - Q 2 - Sp 3 


(H. 1) 


‘Since conductance may be defined (Dushman, Ref. 2, p. 82) as throughput divided by 
p ressure drop, the following equations hold: 


Q 1 = F 1 (P 1 ' p 2> 


(H.2) 


% = F 2 (p 2 ' p 3 ) 


(H. 3) 


S ubstituting Eq. H.2 and H.3 into Eq. H.7 , one obtains, respectively. 


F l p l ' F l p 2 = Sp 3 


(H.4) 


and 


F 2 P 2 = (S + F,) p 


2' ^3 


(H.5) 


J-fEq. H.4 is multiplied by (S + F 2 )/S, 


S + F r 


S + F r 


S F 1 P 1 


S F 1 p 2 + P 3 


(H. 6) 


S ubtracting Eq. H. 6 from Eq. H.4, one obtains 


1 + X F 1 + F 2 


Vo = 1 + 


it) F 1»1 


(H.7) 
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If Eq. H. 7 is divided by the large term on the left-hand side, 

(H. 8) j 


If the limiting conductance F 2 is significantly lower than the pumping speed S for 
the species in question, Eq. (H. 8) simplifies to 



Since factors affecting the values of conductance for each species would be the same 
for both F 1 and F 2 , these factors would cancel in Eq. H.9. Thus, partial 
pressure in the mass spectrometer is dependent only on the partial pressure in the 
microbalance jar and on the physical dimensions of the limiting conductances F 1 
and F 2 . This means that no species discrimination will occur as long as the limiting 
conductance F 2 is much less than the pumping speed S for the species in question. 

Response time of the mass spectrometer section may also be determined from the 
above parameters. If response time is defined as the sampling (mass spectrometer 
section) volume times pressure divided by the throughput. 


y 2 

T " Q 


(H. 10) 


Substituting Eq. H. 2 in the above equation one obtains 


T = 


V 2 

(Pi - P 2 ) 


(H. 11) 
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If Eq. H.9 is then substituted into the above equation, 



which simplifies to 


T 



(H. 12) 


(H. 13) 


Again, if the limiting conductance F 2 is much smaller than the pumping speed S , 
the above equation simplifies to 


T 



(H. 14) 


For the system described in Section 3. 1. 2, in which the mass spectrometer section 
volume V 2 is less than 300 ml and the limiting conductance F g is 1. 5 l/sec, the 
response time r is less than 0. 3/1. 5 = 0. 2 sec. 
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Appendix I 

CALCULATION OF THERMODYNAMIC QUANTITIES USED IN DISCUSSION 

IN SECTION 4.2 


Calculations of thermodynamic properties of n-butane and toluene from Eq. 4. 35 and 

4.40 are shown in Tables 1-1 and 1-2, respectively. The term 0 in Eq. 4.35 may be 
defined (see Eq. 4.33) by 


<p = R/m 


where 


R = gas constant 1.987 (cal/mol°K) 

m - slope of the tangent to the potential-plot curve (ml liq/mol 0 K) 


( 1 . 1 ) 


The isosteric heat of adsorption calculated here is defined by Eq. 4.35. The values 

for the heat of vaporization, AH y , used in these calculations are 5.5 kcal/mol for 
n-butane and 9. 1 kcal/mol for toluene. 
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Table 1-1 


CALCULATION OF THERMODYNAMIC QUANTITIES FOR n-BUTANE ADSORPTION 
AT 300°K BASED ON THE POTENTIAL-PLOT CORRELATION (Fig. 4-2) 


Adsorption pressure, p (torr) 

1.3 x 10-3 

1.7 x 10~ 2 

8.8 x 10" 2 

1.3 x 10° 

Amount adsorbed, q (ml liq/g) 

1.0 x 10" 2 

1.0 x 10" 1 

1.7 x 10' 1 

3.5 x 10-1 

Intercept at A = 0 of tangent to 
potential-plot curve at q , q° 
(ml liq/g) 

8.8 x 10° 

3.0 x 10° 

1.7 x 10° 

1.06 x 10° 

q/q° 

1.14 x 10“ 3 

3.33 x 10” 2 

1.0 x lO" 1 

3.30 x 10" ] 

ln(q/q°) 

-6.78 

-3.40 

-2.30 

-1.11 

1 /slope of tangent to potential-plot 
curve at q , l/m (mol°K/ml liq) 

-7.1 

-8.6 

-9.5 

14.2 

(V m R/m)'ln(q/q°) (kcal/mol) 

9.2 

5.6 

4.2 

3.0 

Isosteric heat of adsorption, 
<q S t)th (kcal/mol) 

14.7 

11.1 

9.7 

8.5 

Fractional coverage , 9 

1.81 x 10 -2 

1.81 x 10" 1 

3.1 x 10 _1 

6.3 x 10" 1 

Entropy of adsorption, 
(qst)th/T (cal/mol°K) 

49.0 

37.0 

32.3 

28.3 

-R In 9 (cal/mol°K) 

9.66 

3.395 

2.32 

A 

0.92 

Adsorption pressure, p (atm) 

1. 7 x 10 -6 

2.2 x 10 -5 

1.16 x 10”* 

1.7 x 10~' i 

-R In p (cal/mol°K) 

Theoretical adsorption entropy, 

26.4 

21.3 

18.0 

12.7 

-AS a (cal/mol°K) 

53.0 

41.6 

37.2 

30.5 
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Table 1-2 


CALCULATION OF THERMODYNAMIC QUANTITIES FOR TOLUENE ADSORPTION 
AT 300°K BASED ON THE POTENTIAL- PLOT CORRELATION (Fig. 4-3) 


Adsorption pressure , p (torr) 

— 

ca. 1 x 10~ 5 

3.5 x 10~ 4 

4.5 x 10 ~ 2 

Amount adsorbed, q (ml liq/g) 

1.83 x 10“ 2 

4.98 x 10-2 

1.35 x 10" 1 

3.68 x 10" 1 

Intercept at A = 0 of tangent to 
potential-plot curve at q, q° 
(ml liq/g) 

4.90 X 10° 

4.0 x 10° 

2.0 x 10° 

1.0 x 100 

q/q° 

3. 74 x 10 -3 

1.24 x 10 -2 

6.75 x lO -2 

3.68 x 10" 1 

In (q/q°) 

-5.59 

-4.38 

-2.69 

-1.00 

1/slope of tangent to potential- 
plot curve at q, 1/m 
(mol°K/ml liq) 

-8.54 

-8.90 

-10.9 

-16.5 

(V m R/m) In (q/q°) (kcal/mol) 

11.20 

9.17 

6.90 

3.88 

Isosteric heat of adsorption, 
(<lst>th( kc al/mol) 

20.3 

18.3 

16.0 

13.0 

Fractional coverage, 9 

3.17 x 10" 2 

8.63 x 10“ 2 

2.34 x 10” 1 

6.38 x 10 _1 

Entropy of adsorption, 
( q st>th/ T (cal/moFK) 

67.7 

61.0 

53.3 

43.3 

-R In 6 (cal/mol°K) 

6.85 

4.87 

2.88 

0.89 

Adsorption pressure, p (atm) 

— 

ca. 1. 3 x 10~ 8 

4.6 x 10"7 

5.9 x 10" 5 

-R In P (cal/mol°K) 
Theoretical adsorption entropy, 

— 

ca. 36. 0 

29.0 

19.4 

-AS a (cal/mol°K) 


58.0 

49.0 

37.3 
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Appendix J 

DERIVATION OF EQUATION 4.45 


In this section Eq. 4. 45 is derived by a kinetic argument, and it is shown that the 
Langmuir theory of mixed adsorption and the Freundlich theory of mixed adsorption 
(as derived by Glueckauf , Ref. 20) are special cases of this equation. 

The following simple kinetic argument leads to Eq. 4.45: 

Consider an adsorbing surface covered to a fraction e 1 by component 1 and to 6 2 by 
component 2 , with equilibrium gas phase pressures p 1 and p 2 . For each compon- 
ent at equilibrium , the rate of evaporation must equal the rate of adsorption. The rate 
of adsorption is set equal to the collision rate times a sticking coefficient: 

dq/dt = ^P 1 /y27rmkT |o! ( J - 1 ) 


where 


q 

t 

P 

m 

k 

T 

a 


amount adsorbed 
time 

pressure 
molecular mass 

Boltzmann constant = 1.380 x 10 

temperature 

sticking coefficient 


2 

(molecule/cm ) 
(sec) 

(dyne/cm ) 

(g/molecule) 

(erg/molecule°K) 

(°K) 


The rate of desorption is set equal to a frequency factor k Q times an Arrhenius activa 
tion energy term times the fraction of surface covered: 


-dq/dt = 


k Q exp 


Q(0 t )' 

RT 


(J.2) 
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where 


Q(0-p) molar heat of adsorption at total coverage = 0^, (Kcal/mol) 
® \ — fraction of surface covered by component 1 

Equating these rates gives a Langmuir-type equation: 


2ir mkT 


= k Q exp - 


Q(0 t ) 


where the subscript m ' on P m j denotes that it is the equilibrium pressure in a 
mixture. 


There are several assumptions in Eq. J. 3. Q(0 T ) is written to indicate that the heat 
of adsorption of component 1 depends only on the total surface coverage 6 = 9 + e 

We picture here a heat of adsorption which decreases with surface coverage and de- 
pends only on the total amount adsorbed, roughly independent of the composition of the 
adsorbed layer. Tacitly a uniform surface is also assumed since the rate of evapora- 
tion is taken to be proportional to ^ and not to a sum over a distribution of heteroge- 
neous sites. 


Equation J. 3 applies to adsorption from a mixture. For adsorption from the pure 
component 1 the pressure of pure component 1, p Ql , which gives exactly the same 
total surface coverage 6^ will be given by an analogous equation , 



= k Q exp 


Q(0 t )‘ 

RT °T 


(J.4) 


where we assume that a and k Q are roughly the same for a surface covered by pure 
component 1 as for the mixture. In Eq. J.4 the subscript "o" denotes a pure 
component. 
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Now, dividing Eq. J. 3 by Eq. J.4 gives immediately 


p ml = 9 _1_ 
Pol d T 


(J.5) 


Since 9^/6^ is the mole fraction of component 1 in the adsorbed mixtures (assuming 
equal areas per molecule), then rearrangement of Eq. J. 5 gives 


which is the desired equation. 

To show that the above equation holds for mixed Langmuir adsorption it is to be noted 
that the normal Langmuir equation for a single gas , 


6 


ap x 

1 + ap x 


may be rewritten 


P 


I g 

a 1 - 9 


(J.7) 


(J • 8) 


Thus, the pressure of the pure gas required to give a particular adsorption <? T is 


1 9 T 

PqI^T* a 1 - 0 T 


(J • 9) 


Here "a” is the ratio of rate constants for adsorption to those for desorption in the 
original Langmuir derivation. 
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Now for mixed adsorption it is basic to the Langmuir theory that the pressure of com 

ponent 1 necessary to give an adsorption ^ in the presence of an adsorption 0 of 
the second component is ^ 


ml 


a 1 - *1 - *2 


a 1 - 0. 


(J • 10) 


where 


0 T 9 1 + 6 2 


Now if Eq. J. 9 for the pure component 1 is multiplied by X = 6 /(0 + 0 ) = g /g 
there results 1 1 1 2 IT’ 


x i p oi<V = s r^ 1 


e r 


Equating the left sides of Eq. J. 11 and J, 12 gives 


^ml ^lPol^T^ 


which is the desired equation (Eq. 4.45). 


(J.ll) 


(J.12) 


We next consider Glueckauf’s equation for mixed Freundlich-type adsorption on a 
heterogeneous surface: 


l0g X™ 1 “ m log < 0 i + d 2 } ~ const * (J. 13) 

which is Eq. 26, page 382 of Young and Crowell (Ref. 14). Here m is the Freundlich 
exponent, which both gases presumably obey. 
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Equation J. 13 can be rewritten 


P ml 


= X 1 A(0 1 + 



(J.14) 


But Ale ♦ e ) 1/m is exactly P ol (0 T ) , the equilibrium pressure the pure gas would 

( 1 2 f q — Q l\ so Eq. J.14 is simply a special case of Eq. J.12 

exert as a coverage of 0 ^ 1 u 2 ' 

(4.45) adapted for the Freundlich isotherm. 


Thus, Eq. 4.45 appears to be a rather general 
models for mixed adsorption. 


equation which fits a number of diverse 
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Appendix K 

DETAILS OF DATA OBTAINED WITH THE VOLUMETRIC DESORPTION APPARATUS 

Data on n-butane desorption rate from BD charcola at 26" C have been summarised in 
Section 1. 1. Some details are presented here. The data from some thirty deter- 
minations are shown in Table K-l. Each run is designated by a symbol such as 
-49, in which the letter designates a certain sample of charcoal. With a given letter 
(say A) a change in number refers to a major change in treatment of that sample, 
usually a repeated bakeout under high vacuum. 

All data for a given sample of carbon and a given initial loading show excellent agree- 
ment with second order desorption kinetics. The most extensive data are found in 
run A-49 m which seven points were obtained (including the first point of A-51). All 
seven points, when plotted as inverse of fraction left vs. time lie on a straight line 
going through the correct intercept at zero time. The point is 75. 6 percent desorbed 
is a bit low, but experimental error in reading the Pirani scale at this high pressure 
is large. In some cases (e. g. , D-72 and E-79) only one point at a given loading was 
token, since second order seemed well established. In these cases the time for half 
desorption (used in Fig. 5. 3. 1-1) was read from a linear second order plot going 
through the one experimental point and through the correct Intercept. 

C alibration of Pirani Oage. The Pirani gage was calibrated for n-buta„e using a 
McLeod gage. The Pirani scale reading had to be divided by 3.33 to give the true 

butane pressure. The calibration was checked almost every day. 

A dsorption Isotherm . Each time butane was admitted to freshly baked charcoal, a 
point on the adsorption isotherm was measured. Just as for the rate data, for a given 
sample of charcoal these data were reproducible, but had a scatter for different samples 
o charcoal. In particular both samples of 0. 1-cm particles adsorbed more strongly 
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Table K-l 


SUMMARY OF DATA FOR n-BUTANE DESORPTION AT 26°C 
t»t tt/^t TiMrTRTP APPARATUS 
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an both samples of 0 . 2-cm particles, so that at a given loading the equilibrium 
pressure was about a factor of two lower for the smaller particles. This odd 
P enomenon has also been observed on the microbalance. These data also check with 
e microbalance data in that, between 3 and 30 pHg pressure, the observed isotherm 
“ linear on a log-log plot (Freundich isotherm) with a slope of 0. 5. However the 
absolute pressures are a factor of about two higher in the present work <even for 0 1- 
cm particles) than in the earlier microbalance data (Fig. 2 - 8 ). This is most likely due 

o c i ra ion errors in one or both researches, the McLeod Gage used in the present 
work being of unknown accuracy. 

rupted Experiment. In run A-51 an interrupted experiment was performed to 
compare with the continuous data of A-49. In A-51 desorption was carried out for 1- 
minute intervals separated by 10-minute waiting periods. The desorbed butane was 
accumulated in the cold trap and was never readmitted to the charcoal. The total 
desorption time for A-51 is, thus, the sum of the individual desorption times, and the 
o al percent desorbed is the sum of the percents desorbed in each increment The 
interrupted experiment was somewhat faster than the continuous experiment, especially 
at high percentages desorbed. Thus, in 5.4 minutes charcoal in the interrupted experi- 
ment desorbed 74. 8 percent while it took 11. 0 minutes to desorb this percent in the 
continuous experiment. This result is unexpected, since neither pore diffusion effects 
or adiabatic cooling effects are to be expected at room temperature. Experimental 

in a h , Urely T ° rder ’ and the P r °P° Bed experiments with a thermocouple imbedded 
charcoal granule may shed further light on this problem. 
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Appendix L 

DETERMINATION OF DESORPTION RATE ORDER 


Many desorption rates have been found to follow kinetics analogous to those of chemical 
reactions in which the rate of the reaction is proportional to a function of the concentra- 
tion of a single reactant. Thus, the rate expression for desorption can be written: 

- dq/dt = kq 11 

where q = amount remaining adsorbed 
t = time 
k = rate constant 
n = rate order. 

This rearranges to 

which integrates to 
Setting boundary conditions, 

A characteristic time, r , may be defined such that, 

a = aq^. (0 < a < 1) . (L- 5 ) 

For instance, if the characteristic time were defined with a = l/2 , it would then be 
the half-life of the desorption. Since most desorption data cover a fairly limited range 
of relative amounts remaining adsorbed, usually from 10 to 100 percent, using the half- 
life would yield only three or four points for the following analysis. A better choice 
would be a = 10 -0 * 1 = 0.7943 which would give 10 points per decade. 


-q aq - kdt , (^-2) 

q 1-n /(n - 1) = kt + c (if n * 1). (c = constant) (L.3) 

q = q° at t = 0 , 

q 1-n /(n - 1) = kt + q° 1_n /(n - 1) . ( L * 4 ) 
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Then, at some time = t + r , 


(aq) 

Subtracting Eq. (L. 1) 

Then, log (a 1-n - 

Differentiating: 


1 n /(h - 1) = k(t + r) + q° 1-n /(n - 1). 

(a - l)q* n /(n - 1) = kr . 

1) + (1 - n) log q - log (n - 1) = log k + log r . 
( 1 - n)d log q - d log t . 


(L. 6) 

(L. 7) 

(L. 8) 
(L. 9) 


Thus, the slope of a plot of log t vs. log q , would have a slope = (1 - n) . (For the 

special case of n = 1 , it can be easily shown that r is independent of q, so that the 
slope would equal zero, which is still n - 1. ) 


Further information can be gained from the above plot since, at t = 1 (log t = 0) , 
Eq. (L,8) becomes 

1 n 

log k = log (a - 1) + (1 - n) log q' - log (n - 1). 

The rate constant may be estimated by the above equation, which rearranges to: 

k = (a 1-n - ljq' 1 "^ - 1) (L 

where n = rate order 

q' = intercept, t = 1 (log t = 0). 


Alternatively, at q = 1 (log q = 0) , Eq. L. 8 becomes 

log k = log (a 1 n - 1) - log r' - log (n - 1) , (L. 12) 

which on rearrangement yields 


k - (a 1 n - l)/(n - l)r’ (L.13) 

where 

r' = intercept at q = 1 (log q = 0) . 
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For a series of desorptions of slightly varying order, it would be necessary to approx - 
mate the slope of a single order for each set of data points, in order to determine the 
rate constants in this manner. Thus, the method is somewhat less precise than that of 
making a plot according to the applicable rate equation. 

The method of making such a rate order-determining plot is as follows. Values of a , 

with b varying from 0 to a number sufficiently high such that a b = 0. 1 are determined, 
preferably from mathematical tables. The value of amount adsorbed at or very close to 
the beginning of the desorption is multiplied by each value of a , giving various 
amounts remaining adsorbed that differ by a factor of a . The calculation of amount 
adsorbed is reversed so that the weight increase, and thus the total weight, of the 
desorption sample is known for each value of a" b . The time required to reach each 
of these successively smaller weights is then determined directly from sample weight 
recordings. (If the desorption is measured by other than a gravimetric method, some 
analogous method of determining each successive desorption time must be used. ) Next, 
the difference between each successive desorption time is determined. This difference, 
r , is plotted on log-log paper against the amount remaining adsorbed just prior to the 
beginning of each time period. The slope of the plot = 1 - n; the intercept at t - 1 
(log t = 0) = log q' at r = 1 . The rate order = n . Substitution of n and log q' at 
t = 1 in Eq. L. 11 then yields the rate constant. 
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Appendix M 

KINETIC THEORY OF VACUUM DESORPTION 


In this section a more detailed kinetic picture of the theory is developed, taking into 
account explicitly the shape of the adsorption isotherm and the pumping speed. 

There are clearly three rate processes which determine the overall rate of desorption 
from a single particle: 1) rate of diffusion of an adsorbed molecule in the pore struc- 
ture to the particle surface, 2) rate of the actual desorption step from the adsorbent 
surface into the gas phase, and 3) pumping speed. 


Diffusion within the pore structure can be either Knud sen gas diffusion in the sm all 
pores or surface migration of the adsorbed molecule along the pore wall. 

M. 1 CASE I: GAS DIFFUSION IN PORES 

Diffusion of gas within the pore structure when accompanied by adsorption is ruled by 
the equation (for spherical granules): 

9c as 

9t + at ( M - 1 ) 


3-c 2 8c 

8x“ 2 ' ?9lT 


where c and S are the gas-phase and ads orbed-phase concentration, respectively, 
D is the diffusion coefficient and r is the radial distance from the center of the 
granule. Since S is enormously greater than c for a good adsorbent, the term 
9c/9t can be dropped'from Eq. M. 1 . Since there are two unknowns (c and S), a 
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second equation is necessary and this is the equation of adsorption kinetics on the pore 
wall: 


as 

at 



(M. 2) 


where Langmuir type kinetics have been assumed, k^ and k^ are rate constants for 
adsorption and desorption, and S 1 is the adsorbed concentration at one monolayer. 

Equations M. 1 and M. 2 are to be solved for boundary conditions determined by the 
pumping speed. Before investigating this, however, it is noted that if adsorption is 
very fast, then adsorption equilibrium can be assumed throughout the granule, so that 
Eq. M. 2 can be replaced by the equation of the adsorption isotherm: 


S = f x (c) 


(M.3) 


Thus 9S/9t in Eq. M. 1 can be replaced by 


as 

at 


df 2 ( c ) 

dc 


9c 

at 


(M. la) 


so that Eq. M.l becomes: 


d; 



2 9c 
r 9r 


«!( c > 8C 

dc at 


(M. 4) 


Clearly the derivative df^/dc depends on the shape of the adsorption isotherm and is 
just the slope of the isotherm, dS/dc. For example, for the Langmuir isotherm it is 
easy to show that: 


df j(c) S Q R q (M.5) 

” c o lR o + (1 - R o >c/c o ] 
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where S q is the initial adsorbed concentration in the granule and c q is the gas phase 
concentration in equilibrium with it. R q is a parameter determining the fraction of 
surface covered before desorption begins. When R q = 1 , the surface is only slightly 
covered, and the isotherm is linear. When R q = 0 , the surface is fully covered and 
the isotherm is very non-linear, i.e., irreversible adsorption). Thus for the Langmuir 
isotherm, the final diffusion equation to be solved is: 

D ( 9^.9 H . 2 £c\ _ R 0c 

\ Sr 2 r dr J " c q R + (1 - R)c/c q at ( M - 6 ) 

This equation is non-linear, except when R = 1 , so that numerical integration is 
necessary. 

An alternative form to Eq. M-6 is to eliminate c rather than S in Eq. M.l using the 
the adsorption isotherm. Thus if the isotherm is inverted to give c as a function of 
S : 


C = f 2 (S) 


(M. 7) 


then 


dc _ df 2^ S ^ 3S 
8r dS 9r 


and Eq. M.l may be rewritten as: 


1 _ 

2 


_9 2 df 2 (S) as 

8r 1 D dS ar 


as 

at 


(M. 8) 


(M. 9) 
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This is exactly equivalent to diffusion of a single component in a medium in which the 
diffusion coefficient depends on concentration > 


D x (c) 


D 


df 2 ^ 

dS 


(M. 10) 


We note that if Eq. M. 9 had been written for a plane sheet, rather than a sphere, it 
would be 


D 


_ 9 _ 

ax 


df 2 (S) 

"~dS~ 


as 

8x 


as 

at 


(M. 11) 


The importance of Eq. M. 11 is that Crank (Ref. 30, Ch. 12) gives numerical solutions 
for a number of problems with variable diffusion coefficients. Each of Crank’s func- 
tions which give this variation may be converted into an equivalent isotherm, and thus 
much numerical work may be avoided. 


The boundary conditions for which Eq. M. 4 or M. 9 must be solved depend on the 

pumping speed and warrant some investigation. First of all, if pumping speed were 
infinitely fast, then the boundary condition would be that at the pellet surface the con- 
centration of gas (and a sorbed phase) is zero. For a finite pumping speed , the 
boundary condition depends on the particular experimental set-up. A rough and ready 
intuitive boundary condition would be that the outward flux of gas from the granules 
equals the pumping speed times the concentration at the surface. This clearly is sort 
of a steady state approximation which neglects the initial transient fall in pressure 
when the pumps are first turned on. 

A mathematically and experimentally exact solution of this troublesome boundary 
problem is as follows: Experimentally, the charcoal granule is separated from the 
pumps by a tube of radius ^ (^ = about 1. 8 cm for the microbalance). If at distance 
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L, from the granules, there is located an ion gage (or other fast response manometer), 
then at zero time, just before the pumps are turned on, this gage will read c , the 
initial pressure in equilibrium with the adsorbent. The pumps are now turned°on and 
the ion-gage reading is recorded as a function of time, and will rapidly decrease by 
some law, say. 



(M. 12) 


where c f is the ultimate pressure and b is an empirical time constant depending on 
the pumping speed. We now know empirically the pressure history at this one point 
in the apparatus. (We note that in Eq. M. 12 c refers to the partial pressure of 
adsorbate, and not to the total pressure including water vapor etc. ). To relate this 
known pressure history at a distance L from the adsorbate to that at the adsorbate 
one must solve the law of flow in the connecting tube. At low pressure around one 
micron this will be Knudsen flow (or diffusion), and hence one must solve the diffusion 
equation in the connecting tube simultaneously with Eq. M.4. That is, the pumping 
or diffusion equation in the connecting tube is 

p. 9_ c gc 

D K 0Z 2"8t < M - 13 ) 

where is the Knudsen diffusion coefficient in the tube and Z is distance from the 
adsorbate. Equation M. 13 must be solved with the boundary conditions: 

a) At t = 0 , c = c q for all Z 

b) At x — L, c is a known empirical function of time, for example given by 
Eq. M. 12. 

c) At x = 0 the flux must equal the rate of desorption. The rate of desorption 
for one granule is the outward flux through the pellet surface, 

- 4 " R2D (fr') s 
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(here R is the granule radius) while the flux through the connecting tube at x = 0 is 


- 7rr 


2 

1 




(M. 15) 


so that a boundary condition for both Eq. M.4 and Eq. M. 13 is 


+ N4ttR 2 



2 

7rr l D K 



(M. 16) 


where N is the number of adsorbent granules used in this experiment. 

Thus the exact solution of the pumping speed problem is to solve Eq. M. 4 and Eq. 

M. 13 simultaneously with Eq. M. 16 as a common boundary condition. This is likely 
quite possible but it is probably not necessary for a very good approximation, due to 
the very small time constant involved in Eq. M. 13. The time constant for Eq. M. 13 
is clearly L 2 /D k , where L is the connecting tube length and D^. is the Knudsen 
diffusion coefficient in the connecting tube. For microbalance, the connecting tube is 
about 35 cm long and about 3. 5 cm in diameter. The Knudsen diffusion coefficient in 
a tube of this size is 


D K = 9-7 X 1(>3 X T 1\M = 3 ’ 8 X 1<)4 (M. 17) 

Thus the time constant L 2 /D is about 3 x 10 2 seconds. This means that the gas in 
the connecting tube will adjust very rapidly to changing concentrations at the boundaries 
of the tube providing the time constant of these changes is long compared to 3 x 10 
seconds. Since the time constant for desorption rate is roughly a thousand times 
slower, the connecting tube will be in a pseudo-steady state with a linear concentra- 
tion gradient. This linear gradient will be 


9c 

az 


Co - 


L 


(M. 18) 
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where Cg is the concentration at the pellet surface (Z = O) and is the concentra- 
tion at the end nearest the pumps (Z = L). Putting Eq. M. 18 into M. 17 gives the 
desorption rate: 


- d < S / S o> _ ”l D K , 

dt L (C S c f ) 


(M. 19) 


The multiplying factor 


L (M. 20) 

is exactly the pumping speed of the connecting tube of length L . Thus Eq. M. 19 is 
the theoretical justification for the intuitive guess mentioned earlier, i.e. that a prac- 
tical boundary condition might be to equate the rate of desorption to the pumping speed 
times concentration at the pellet surface. 

Thus the theoretical rate of desorption including pumping speed effects can be carried 
through m either of two stages of approximations: a) the rigorous method of solving 
Eq. M. 4 and Eq. M. 13 simultaneously, which will include all transient pumping effects 
or b) using Eq. M. 19 as the external boundary condition for the pellet, which will 
neglect the initial transient pumping effects at the very start of the experiment. 

A third and most practical method of analyzing for effects of pumping speed is to test 
a given experiment after it has been performed. Since the rate of desorption is 
measured in the experiment, this empirical quantity can be put into Eq. M. 19 and the 
concentration at the pellet surface can be computed, since c g is the only unknown in 
Eq. M. 19. For the pumping speed to be great enough so that it does not affect the 
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desorption rate , the computed value of c g must be much less than the gas concen- 
tration in equilibrium with the granule at this stage of desorption. Using Eq. M. 19, 
this criterion works out to be 


-d(S/S Q )/dt 



(M. 21) 


where 

3 

S = pumping speed of the connecting tube (cm / sec) 

p = gas pressure which is in equilibrium with the adsorbent at this stage of 
6 

desorption (atm) 

Pj = pressure at a distance L from the absorbent (atm) 

Equation M. 21 is rapid and easy to apply to a given run. Because desorption rates are 
fast, it turns out that very large pumping speeds are necessary to measure the true 
rate of adsorption. This will be discussed later. 

To summarize: For the case of gas phase diffusion being the main method of trans- 
port in the pores, and for a system with very fast pumping speed, desorption rate 
should be ruled by Eq. M. 4 or M. 9 solved with the boundary condition c = O at the 
granule surface. Whether or not pumping speed is fast enough to make this the 
experimentally determined rate should be tested by using Eq. M.21. 

M. 2 CASE II: ADSORBED PHASE DIFFUSION IN PORES 

Before discussing solutions of Eq. M.4, the second case for which diffusion in the 
adsorbed phase is the main mechanism of flow within the granule is considered. Here 
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the gas phase in the pores is neglected, and in an elementary theory diffusion within 
the granule would be ruled by the simple diffusion equation: 

„ /a 2 s . 2 as\ _ as 

S \ar 2 r9r / ” “ (M,22) 

solved with the boundary condition that the diffusive flux at the surface must equal the 
rate of desorption from the surface or: 

- D (if ) s = ^ - k A c S < M23 > 

where Freundlich desoiption kinetics have been assumed for the actual surface desorp- 
tion step. In Eq. M. 23 Ng is the number of molecules adsorbed per unit surface at 
the external surface and Cg is the gas phase concentration at the external surface. 

For very fast pumping speed, c g is zero, so the last term in Eq. M. 23 may be 
dropped. 

There are however problems connected with Eq. M.22, for nonlinear isotherms. 
Equation M. 22 assumes that the adsorbed phase diffusion rate should be proportional 
to the adsorbed concentration gradient. This is likely a poor assumption in general, 
since diffusion rates for "imperfect systems" should be proportional to the gradient 
in free energy or activity. Since this quantity is measured by the pressure of gas 
c e (S) in equilibrium with the absorbed concentration S , it appears more reasonable 
to use the basic flux equation than the diffusion rate across a plane: 

9c 

D 1 Tr" (M. 24) 


rather than 


D 


as 
s ax 


(M. 25) 
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If Eq. M. 24 rather than Eq. M. 25 is used, and if c 0 is given by some adsorption iso- 
therm, Eq. M. 7, then the diffusion equation to be solved for adsorbed-phase diffusion 
would be exactly Eq. M. 9 with a changed meaning of the diffusion coefficient. In 
Eq. M. 9 it is the Knudsen gas diffusion coefficient, while in the version of Eq. M.9 
modified for adsorbed phase diffusion it is defined by Eq. M. 24. 

Thus , although the diffusion equations for gas diffusion in the pores and for adsorbed 
phase diffusion are mathematically similar there are three important differences: 
a) for gas phase diffusion the diffusion coefficient D is the Knudsen diffusion coef- 
ficient which is proportional to pore radius. For adsorbed phase diffusion the coef- 
ficient (or Dg) is for two-dimensional migration along the surface, and hence 

will vary roughly inversely with the pore radius , since the surface area available for 
two-dimensional migration will vary inversely with pore radius. Thus we expect 
adsorbed phase (surface) diffusion to be important in very small pores such as those 
found in charcoal, b) Although the equations are formally similar, the boundary con- 
ditions are quite different, since for adsorbed phase diffusion desorption occurs only 
at the granule surface while for gas phase diffusion it can occur from the internal pore 
walls throughout the whole granule. Thus when diffusion is very fast compared to 
desorption rate, for the surface diffusion case the rate becomes equal to the rate of 
desorption from the external surface only, while in the other case it becomes equal to 
the desorption rate from the whole internal surface. In particular when adsorbed 
phase diffusion is very fast, the adsorbed concentration will be uniform throughout 
the pore structure , and equal to the average concentration S within the granule so 
that the rate of desorption for a single granule becomes, from Eq. M. 23, 

4« 2 (k E N"-k A c s ) < M - 26 > 


The adsorption rate constant per unit surface, k A , may be eliminated in favor of 
c (S) , the gas concentration in equilibrium with S , so that Eq. M.26 becomes. 

G 


Rate of Desorption 
from one granule 


47rR 2 k E Ng 



(M. 27) 
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When gas -phase pore diffusion is very fast, the external surface (47 tR^) in Eq. M. 27 is 
replaced by the total internal surface of the granule 


Rate of Desorption 
from one granule 


4jr 

3 


R 3 p S k 
y P g E S 



(M. 28) 


where 

Pp = particle density 
S = internal surface area 

B 

M. 3 SOLUTIONS FOR RATE OF DESORPTION 


In this section equations for the absolute rate of desorption will be collected, for the 
various rate-determining cases. Particular attention will be paid to the dependence 

of rate on particle size and on temperature. Freundlich kinetics and isotherms will 
be assumed. 

M. 3. 1 Pumping Speed Limiting 

The simplest case is when pumping speed is limiting. Here, by assumption, the 
desorption rate from the granule is faster than pumping speed, so that each granule 
is surrounded by its equilibrium gas pressure. Hence, the rate of desorption is the 
rate at which gas at this pressure is pumped away. For granules containing a concen- 
tration, S , of adsorbate, this is, in moles per second, 

“ W = (S p )C e (S) (M.29) 


where 

S = pumping speed (1/sec) 

c e ( s ) = equilibrium concentration of gas (mol/1) 
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The number of moles , N , of adsorbate in the charcoal is the volume of activated 
carbon pellets used times the concentration S , or 


N 



(M. 30) 


where w = weight of activated carbon used (g). 
Substituting Eq. M.30 into Eq. M. 29 gives 


P = —E- (S ) c (S) 
dt W P e 


(M. 31) 


If c g (S) is given by the Freundlich isotherm , c g - kS , or its equivalent c g 


c /c 
e o 


<s/s 0 ) 


n 


(M. 32) 


fS c are initial concentrations), then substituting Eq. M.32 in Eq. M.31 gives 

' n n 


d(S/S Q ) 

dt 



(M. 33) 


Here S/S is clearly the fraction remaining. Hence when pumping speed is limiting, 
the kinetics of desorption are n th order, where n is the exponent in the Freundlich 
isotherm , and the effective time constant is 


7 1/2 



(M. 34) 
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This expression is simply equivalent to the number of moles initially adsorbed divided 
by the number of moles pumped away per second at the initial equilibrium pressure. 
Clearly the adsorption rate depends on the number of particles (i.e. , the weight of car- 
bon) used. (This emphasizes the intuitively obvious criterion that for absorption rate 
to be independent of pumping speed, the time constant for desorption must be indepen- 
dent of the number of particles used. ) The apparent activation energy when pumping 
speed is limiting will be the heat of adsorption, since at constant initial loading, S 
the quantity, c q , in Eq. M.32 must have this exponential increase with temperature 

in order to obey the usual Claus ius-Clapeyron relation used for determining the isosteric 
heat of adsorption, 

M 3 2 Fr o m^Lrte rio DeS0rpti0n from Externa l Surface Fed by Fast Surface Migration 

The second simplest case is that of Eq. M.31, which gives -dN/dt , the number of 
molecules per second for one pellet for the above case. Since N = (4/3 7 tR^)S and 
since S = P p S g N g , 


dS/S 

o 

dt 


3k E N o /_S_\ 

P A R \ s o) 


(M. 35) 


where N q is the value of N g in equilibrium with S Q . Because according to the 
Freundlich isotherm 


j n 

'E o 


k„fC = k. c 
A o 


(M. 36) 


Eq. M. 34 may be rewritten: 


dS/S 3 k. c /_\ n 

dt R S y s ) (M. 37) 
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where is the rate 


constant for adsorption per unit surface. 


Equation M. 35 (or M. 37) predicts a rate varying inversely with particle size and with 


an 


activation energy equal to the heat of adsorption. 


M.3.3 Pore Diffusion Rate Controlling. 

For the case of pore division rate contro.iing Eq. M.4 or Eq M.9 must he solved with 
the boundary condition c = S - O at the external surface. (It is assumed that Eq. M.9 
holds for both gas diffusion and surface migration, so that we will not distinguish 
between these mechanisms in this section. ) A considerable amount of ‘“-can 
be deduced without even solving the equations , as follows. Both Eq. M an . 
involve the derivatives which are the slope of the isotherm or the slope of i 
isotherm . Because we are only concerned with that portion of the isotherm up c Q 
(the initial equilibrium gas concentration) it is convenient to write the isot erm in 

reduced form. 


S 

s' 


= f 3^ c/c 0 ) 


(M. 38) 


, „„ f or the Freundlich iso- 

where f 3 (c/c Q ) is some simple function such a ( / Q ) , Q and must 

therm The function f 3 (c/c ) clearly must equal zero when c/ Q ' 
equal unity when c/c c - 1 . When the slope of the unreduced isotherm dS/dc 

taken to insert, say, in Eq. M.4, it will always have the form of Eq. M.5 such a 


dS/dc = — x a function of c/c q 
c 

o 


(M. 39) 


« follows that for any given adsorption isotherm . the diffusion equation (Eq. M. 5 or 
M.9) can be written in reduced form , involving only c/c c and its derivatives , p us 
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single reduced-time parameter, 



(M. 40) 


The fraction desorbed as a function of time will depend only on this parameter, but will 
be a different function of this parameter for each shape of adsorption isotherm . It 
follows immediately that the time constant (say the time to reach 50% desorption) will 
be proportional to the square of the particle size, and that for gas diffusion at least, 
the activation energy will be the heat of desorption, by the usual Clausius -Clapeyron 
argument concerning the ratio c q /S o . It may be noted that the effective value of the 
diffusion coefficient, according to Eq. M.40 is lowered by the huge factor of c /S 
(about 10 to 10 for a good absorbent) since 10 to 10 more molecules must diffuse 
out than in the absence of adsorption. (It may be further noted that if, for the surface 
migration case, Eq. M. 22 is used then this huge factor does not arise, but if Eq. M. 22 
is used for non-linear isotherms , then Dg will be so concentration and temperature 
dependent as to be useless. ) 


In summary, the essential purpose of detailed solution will be to obtain absolute plots 
of the fraction desorbed vs. the parameter r of Eq. M.39. Detailed solutions are 
available for the following isotherms: 


M.3.3.1 Linear Isotherm: For the linear isotherm (i.e. , R = 1 in Eq. M. 6) the resulting 

equation is linear and has the classical solution given by Crank (Ref. 30, p. 90, Fig. 6.4, 

curve marked "Zero"). Crank plots the fraction desorbed vs. the square root of a 

parameter which is identical with our r , but in his nomenclature lacks the c /S 

o o 

factor. Crank’s "a" is our radius of the sphere , R . Rate of desorption is almost 
linear with Vr up to 50% desorption and then bends over. The value of r at 50% de- 
sorption is (0. 175) 2 = 0. 0306 so that the time for half desorption (sec) is: 

TJ 2 S 

fc l/2 = °- 0306 d'^ (M. 41) 

o 
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Although recently Eq. M.4 has been numerically solved for adsorption for the Langmuir 
isotherm (see Weisz, Ref. 31; also, see Hall, Ref. 32), no solutions for desorption for 
the Langmuir isotherm seem to exist. (For nonlinear isotherms, the results for adsorp- 
tion are not applicable to desorption, due to the changed boundary conditions which do 
not transform into the other problem. ) The only other known solutions applicable to 
desorption appear to be those given by Crank (Ref. 30) for the problem of diffusion 
coefficients which depend on concentration, which may be transformed into the present 
problem. 


Crank (Ref. 30, ch. 9 and 12) gives solutions for the boundary conditions of desorption 
of the equation; 



(M. 42) 


which is for desorption from a plane sheet or slab in which the diffusion coefficient de- 
pends on concentration through the function , f 4 (S/S Q ) . Crank gives solutions for vari- 
ous elementry forms of the function, f 4 (S/S Q ) , such as exp (kS/S Q ) . To apply to 
j the present problem it is first necessary to find out what isotherms the function f 4 

' corresponds to, and secondly to convert to the geometry of desorption from a sphere. 

The first problem can be solved exactly and the second to a high degree of approximation 

. The first problem is exactly solved by the theorem: 

To each Crank function, f 4 (S/S Q ) , there corresponds a reduced inverted isotherm: 



(M. 43) 
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where y - S/S Q . The reduced inverted isotherm, Eq. M.43, has the correct proper- 
ties that c/c q = 0 where S/S Q = 0 , and that c/c q = 1 when S/S Q = 1 . Further- 
more, c/c q has the important property that 


or 


9 c/c 

o 

d~S/s'~ 

o 


f 4 ( E/S o> 

f f 4 <y) <Jy 


dc 

dS 


% J _ 4 (S/S Q ) 

S 1 

° / f 4 (y) dy 


(M. 44) 


If Eq. M. 44 is now inserted into Eq. M.ll [remembering that f 2 (S) in Eq. M. 11 equals 
c in Eq. M.44, because of Eq. M.7] , then it follows that 


D 



< S/S o> 


9S 

9x 


!] 


as 

at 


(M. 45) 


Comparing Eq. M. 45 with Eq. M.42 shows that the solution for desorption from a plane 
sheet using the (inverted) isotherm , Eg. M. 43 , is exactly the same as Crank’s solution 
for the function, f 4 , providing that for Crank’s D q there is substituted 

c 

T) = D — 

o SI (M . 46) 

0 f f 4 <y) d y 


We now apply Eq. M.43 and Eq. M.46 to some results of Crank. 
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M. 3. 3. 2 Logarithmic Isotherm 

A function, f 4 (S/S Q ) , for which Crank gives solutions of Eq. M.40 is 

f 4 ( s / s 0 ) = ex P ( k S / S Q ) 

According to Eq. M.43, the inverted isoterm this corresponds to is 

exp (k S/S q ) - 1 
c/c o exp (k) - 1 

Solving this for S/S Q gives 

S/S q - ~ log [1 + (exp k - 1) C/C q ] 


(M. 47) 


(M. 48) 


(M. 49) 


which may be called the logarithmic isotherm for present purposes. 

Crank (Ref. 30, p. 279, Fig. 12.17) gives solutions for two values of k (k = In 25 = 
3.22 and k = In 10 = 2.303) plotted as fraction desorbed vs. 



(M. 50) 


where l is the half thickness of the slab. Since for this case the integral in Eq. M.46 
is 


£ (exp k - 1) 


(M. 51) 


then Crank's parameter corresponds to plotting fraction desorbed vs. 
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!i/2 


n C o k t_ 

S q (exp k - l) f 2 


(M. 52) 


Thus, for the two adsorption isotherms defined by Eq. M.48 (or M.49) with k = 3. 22 
and 2. 30, we have solutions for fraction desorbed as a function of the parameter given 
by Eq. M. 52, but for a slab of thickness, 2JL , not for a sphere of radius, R . 

To deduce the time to desorb from a sphere from the time to desorb from a plane 
sheet or slab, we introduce the following approximate theorem: If a fraction, 
is desorbed from a plane sheet in time , t^ , then in the same time , t , it will desorb 
from a sphere a fraction, F 2 , given by 

F 2 = 1 ' " F x) 3 (M. 53) 

This approximate theorem is based on the geometry of a sphere and the physical fact 
that the time to diffuse a distance, x , does not depend on the shape of the granule. 
Thus, if molecules diffuse a radial distance of 20% into a sphere, they have occupied 
50% of the volume. Equation M. 53 is simply a geometrical illustration that for equal 
diffusion distances (which should occur in equal times) the fractional volumes outside 
these distances are given by Eq. M. 53. For the only case available, namely the linear 
isotherm , Eq. M. 53 is accurate to a few percent. That is, the rate of diffusion into 
or out of a sphere can be calculated to a few percent from the rate into a slab by 
Eq. M. 53 , as is seen by comparing Fig. 4.6 with Fig. 6.4 of Crank (Ref. 30). 

For the moment we are interested in only one main number from Eq. M. 53, namely 
that when 50% is desorbed from a sphere only 20. 5% is desorbed from a slab. From 
Crank's data when 20. 5% is desorbed from the slab, this Vt is 0. 080 and 0.11, for 
k = 3.22 and k = 2.303, respectively. Hence, for a sphere, when 50% is desorbed, 
t is the square of these quantities or 6.4 x io ~ 3 and 1.21 x io~ 2 . Since t is defined 
by the bracket of Eq. M. 50, we have the relation 
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(M. 54) 



2 

exp k - 1 R 
T 1/2 k D 


S 

o 

c 

o 


where tjy 2 is the time for 50% desorption from a sphere, and Tjy 2 is tlie value of 
Crank's parameter at 20. 5% desorption from a slab. When k = In 25 = 3. 22, then 
(exp k - l)/k = 7.46, and Tjy 2 = 6.4 x 10 . So for k = 3.22: 


t 


1/2 



(M. 55) 


When k = In 10 = 2.303 , then (exp k - l)/k equals 3.91 and: 


t^ = 4.73 x 10 


-2 


„2 S 
R o 

D c 


oJ 


(M. 56) 


Equations M.40, M. 55, and M. 56 now all give the half-desorption times for three dif- 
ferent adsorption isotherms. It is seen that as the isotherm gets more non-linear the 
time for half-desorption gets longer, since the numerical coefficients in front of the 
three equations increase with degree of non-linearity. 

So far we have not discussed the shape of the curve, percent desorbed vs. time, since 
interest has been focussed on the absolute value of the half-desorption times as a 
function of the important parameters . We are now in a position to make our impor- 
tant deduction concerning the shape of this curve: All of Crank's curves for percent 
desorbed from a slab vs. VF are linear up to 50% or so desorbed, independent of the 
shape of the adsorption isotherm. This is a basic characteristic of diffusion into a 
slab. Referring to Eq. M. 53, this means that 


F 


1 



(M. 57) 


where A = the desorption curve shape parameter, 
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holds up to Fj = 0. 50 or so, so that 

Fg = fraction desorbed from a sphere 

= 1 - (1 - At 1 /^ 3 (M. 58) 

should hold up to very high percentages desorbed from a sphere. Thus the simple de- 
duction can be made that if pore diffusion rate is ruling, the fraction F desorbed 
from a spherical granule should be ruled by: 


F 2 = 1 - (l - At 1/2 ) (M. 59) 

from zero to 80 or more percent desorbed. Equation M. 59 is very easy to test experi- 
mentally. It may be rearranged to give: 


1 - (1 - F 2 ) 1/3 = At 1/2 (M. 60) 

1/3 

so that a plot of log [1 - (1 - F 2 ) ] vs. log t should be a straight line of slope 0. 50. 

In preliminary plots this law seems to hold very well, especially at the lower tempera- 
tures. At the highest temperatures the slope tends to be less than 0. 50. It should not 
be minimized that the desorption curves which are of complex shape can be fitted by a 
single parameter A in Eq. M. 60. 
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Appendix N 

DERIVATION OF EQUATION 6. 9 


In this section Eq. 6.9 is derived from first principles. The rate of desorption from 
a single granule (assumed spherical) may be formulated as follows: If the concentration 
of adsorbate in the granule is S , then the total number of moles in one granule is the 
volume of the granule times S or 


n 



(N.l) 


where 


n 


"P 

S 


number of moles in one granule (mol) 

granule diameter ( cm ) 

concentration of adsorbate in the granule (mol/cm 3 ) 


The rate of desorption is clearly the time derivative of this , obtained by writing 
- dS/dt for S in the above expression. 


If the rate of desorption equals the rate of evaporation from the external granule 

2 

surface, then this rate is the external surface area (7r d p ) times the rate per unit 
area. The rate equation is thus 



6 dt 


,2 

7rd r. 
P A 


where 

r A 

t 


desorption rate per unit area 
time 


(mol/cm 2 -sec) 

(sec) 


(N.2) 
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When the granule is in adsorption equilibrium, the rate of adsorption per unit area is 
the rate that molecules strike the surface times a sticking coefficient, a , or 

r A = I C e (S) “ ( N - 3 ) 

where 

c e (S) = vapor concentration in equilibrium with the adsorbed 

concentration S (mol/cm 3 ) 

v = average molecular velocity (cm /sec) 

At equilibrium the rate of adsorption equals the rate of desorption. It is now assumed 
that the rate of desorption into high vacuum is the same as the rate in the presence of 
the equilibrium vapor (i.e. , that the presence of a low pressure of vapor does not 
affect the rate that molecules leave the surface). Thus, Eq. N.3 is inserted for the 
bracket of Eq. N. 2. This gives 


dS 3 1- 

" dt = 2d" VQ!C e (S) < N ' 4 > 

P 

If S is now eliminated in favor of the more practical variable, q (whose units are 
ml STP/g), and if c g (S) is eliminated in favor of the more practical variable of 
pressure in torr, then Eq. 6.9 immediately results. 

Langmuir (cf. Dushman, Ref. 21, p. 18) has used a similar argument to estimate the 
vapor pressure of metals from the rate of evaporation of hot wires. Here the argument 
is inverted, so that the rate of evaporation may be estimated from the equilibrium 
pressure. 


260 



Appendix O 

CALIBRATION OF MANOMETER SYSTEM 


The manometer system, shown in Fig. 0-1 and Fig. 7-2, was used in the vacuum 
desorption system (see section 7.2.1). The desorbed material, collected in the main 
traps, was transferred through stopcock No. 3 and condensed with liquid nitrogen in 
the collecting tube (see Fig. 0-1). The stopcock was then closed, and the condensed 
material expanded in the volume V 2 . This procedure demands an accurate value 
for volume V so calculations can be made in regard to the amount desorbed. 

Z 

The calibration of the manometer volumes was performed on sides A and B to obtain 
accurate volumes (see Fig. 0-1). 

The calibration procedure was as follows: 

(1) The calibrated volume was mounted and evacuated through stopcock No. 2 
on a calibrated volumetric system . 

(2) After evacuation, n-butane gas was admitted into the volumetric system 
and the calibrated volume. The total pressure was read on a Wallace 
and Tier nan FA- 14 5 gage. 

(3) Stopcock No. 2 was then closed and the calibrated volume was transferred 
to the manometer system as shown in Fig. 0-1. 

(4) The manometer system up to stopcock No. 2 was evacuated through the 
main vacuum system and isolated through stopcock No. 3. 

(5) Collecting tube A was immersed in liquid nitrogen and stopcock No. 2 was 
opened . 

(6) All of the n-butane in the calibrated volume was condensed in the collecting 
tube and stopcock No. 3 was closed. 

(7) The liquid-nitrogen trap was removed and the system allowed to come to 
equilibrium at room temperature. 

(8) The pressure of the closed system was read with the mercury manometer. 
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Not to scale. 

Manometer and collecting tube made 
from 4 mm ID precision glass tubing. 
& = high vacuum glass stopcocks. 


Fig. 0-1 Manometer System for the High Vacuum Desorption System 
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Since the total volume of the system varied depending on the total pressure, 
it was useful to use the "zero-point" volume. This is the volume of the 
system when both sides of the mercury column are equal as would be the 
case when both sides are under vacuum or exposed to the atmosphere. 

(9) Given the pressure in the closed manometer system, the volume can be 
calculated from the simple (PV) gas law relationship. Finally, a AV is 
subtracted from the calculated volume (see Fig. 0-1) to obtain the "zero- 
point" volume. 

(10) The same procedure of calibration was used on both sides of the manometer 
system . 

For Run No. F17-1 and thereafter there was a modification of the manometer system. 
Two microstopcocks were glass-blown into the collecting tube with a ball joint in 
between. (See Fig. 0-1.) This allowed one to collect all the sample in the lower 
part of the tube and then disconnect the section for weighing on a balance. 
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Appendix P 

SAMPLE CALCULATION OF VACUUM DESORPTION DATA 


As described in previous sections, the desorbed contaminant was collected in a liquid 
nitrogen trap (see section 7.2.1). At various times t the trap was isolated and the 
condensed material was transferred to the manometer system (see Fig. 0-1, 
Appendix O). Using n-butane as an example, the following calculation was performed 
to determine the amount of material desorbed in time, t 


where 


*de 


= V 


P(V 


m RT x 760 W 


<lde 

v m 

P 

V 

AV 

R 

T 

W 


amount of n-butane desorbed in time, t 
molar volume for n-butane = 96.4 
pressure in the manometer system 
"zero -point" volume 

added volume due to manometer depression due to 
pressure P . Calculated from tube dimensions 
gas constant = 83.0 
manometer system temperature 
weight of activated carbon in bed 


(ml liq) 

(ml liq/mol) 
(torr) 

(ml) 

(ml) 

(ml atm/mol°K) 
(°K) 

(g) 


For example, if the total manometer pressure is 300 torr (manometer depression = 
15 cm), the zero-point volume is 10.95 ml, and the temperature is 297. 2°K, then 


AV = 
V + AV = 

qde = 


<lde = 


7r (0.2) 2 (15) = 1.884 ml (tube ID = 0.4 cm) 
10.95 + 1.88 = 12.83 ml 


(96.4 ml liq/mol) 


(12.83 ml)(300 torr) 

(83.0 ml-atm/mol°K)(297. 2°K)(760 torr/atm)(l. 00 g) 


(96. 4)(12. 83)(300) ml liq/g 
(83. 0)(297. 2)(760)(1 . 00) 


1.98 x 10-2 m i n q /g 
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The same calculation is made for all samples taken at various times t , therefore, 
the total amount vacuum desorbed is 


total 

q de 


I q de 

t=l 


The total amount of n-butane desorbed from the bed is 


= q, 


..total 

de 


Mf 


where 

a = amount desorbed by flow desorption at the conclusion of vacuum 
n df . 

desorption (ml liq/g) 
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Appendix Q 

DERIVATION OF SOLUTIONS FOR EQUATION 7. 7 


In this section solutions of Eq. 7.7 are derived. In particular it is shown that Eq. 7.7 
can be converted to a diffusion equation with variable diffusion coefficients which has 
already been solved by Crank (Ref. 30, Ch. 9 and 12, esp. p. 278 et seq.). 

Equation 7.7 is first rewritten: 


r> JL 9c _ as 

K Ox 3x " at (Q- 2 ) 

A second equation is clearly needed since the above equation has two dependent 
variables. Since by hypothesis desorption from a single particle is faster than 
diffusion through the bed, the bed must be in local adsorption equilibrium. That is, 

S , the adsorbed concentration, varies through the bed (i.e. , is a function of x), but 
each particle is surrounded by a gas pressure of adsorbate which is equilibrium with 
its particular value of S . Thus, the second equation is the equation of the adsorption 
isotherm, and the isotherm equation may be used to eliminate c from Eq. Q.l in 
favor of S . Thus , because the bed is in adsorption equilibrium 

dc _ /dc\ as 

dx Ids/ 3x 2 ) 

e 

where (dc/dS)^ is the slope of the inverted adsorption isotherm (c versus S). If the 
inverted isotherm is written in reduced form , 

c/c o = ff (S /s o ) (Q.3) 
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then 


(i) =^ f 2< s/s o> (Q ' 4 > 

e o 

where f 2 (S/S ) is the slope of the reduced isotherm, Eq. Q.3. Here S q is the initial 

adsorbed concentration in the bed , and c is the vapor concentration in equilibrium 

o 

with S . 
o 

Substituting Eq. Q.2 and Q.4 into Eq. Q. 1 gives 

“Kr^lv^SIJ-f (Q - 5) 

o 

Equation Q.5 may be written in reduced form: 

A l f 2 (y) 111 = (Q - 6) 

where y = S/S q , Z = x/L , L = bed length, and , a reduced time parameter, is 
defined by 


T 


1 


D„ c 
K o 

t 2 S 
L o 


(Q.7) 


Now Eq. Q.5 (or Eq. Q. 6) is exactly the equation for diffusion without adsorption for 
the case when the diffusion coefficient varies with concentration, (see Crank, Ref. 30, 
loc . cit. , Eq. 9.5). Since in general the adsorption isotherm is non-linear, the above 
equations are non-linear, and numerical integration is necessary. However, Crank 
has carried out the necessary tedious numerical work for a variety of fg-functions . 

For Crank the f -function describes how the diffusion coefficient varies with 

Li 
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concentration. For the present work it is the slope of an adsorption isotherm. Thus, 
to every ^-function for which Crank derives solutions there corresponds for the 
present problem an adsorption isotherm which can be derived by integration of Crank's 
fg-function. The exact relationship between the reduced isotherm and the f 9 -function is 

La 


y 

J f 2 (y) dy 

c o 

C 1 (Q-8) 

0 J f 2 (y) d y 

O 

where y = S/S q . Thus, the slope of the reduced isotherm (Eq. Q.8) is proportional to 

f 2 (y) > and the integral in the denominator is a normalizing factor so that c/c is 

unity when S/S is unity. Thus, to get exact analogy with Crank's solution, his D 
u o 

(the diffusion coefficient at zero concentration) must be replaced for the present 

desorption case by an effective diffusion coefficient, 

c 

n — n —2. 1 

o KS q I (Q.9) 

° J f 2 (y) dy 

O 

To illustrate a particular example in detail, one ^-function used by Crank (Ref. 30, 
Fig. 12-17) is 


f Q (S/S ) = exp k S/S 

Li O O 


(Q- io) 


(i.e. , for his case the diffusion coefficient increases exponentially with concentration). 
For the present case this corresponds to the isotherm (obtained by putting Eq. Q. 10 into 
Eq. Q.8), 


c/c 


o 


exp k S/S - 1 
o 

exp k - 1 


(Q- 11) 
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Crank gives his solutions in the form of plots of percent desorbed versus a reduced 
time parameter which for present purposes (using Eq. Q.9) is 



D — — 
S .2 
o L 


kt 

(exp k 



(Q. 12) 


Thus, from Fig. 12-17 of Crank it may be read that, when 50% is desorbed from the 

1 /2 

bed (for the case k = 3.219), the value of the reduced t is 0.223. Hence, the 
value of t is the square of this, or 0.0497, at 50% desorption. Putting these values 
(t = 0.0497 and k = 3.219) into Eq. Q. 12 gives 


0.0497 


% J_ 3.219 

S .2 24 1/2 

o L 


(Q. 13) 


where t^ is the actual time for 50% desorption (not the reduced time). Solving 
Eq. Q. 13 for t^ gives 



0.373 


T 2 s 

L o 

D„ c 
K o 


(Q- 14) 


Now, the time for 50% desorption for other-shaped isotherms will differ from 
Eq. Q. 14 only in the numerical factor, which has been called K in Eq. 7.9 of 
section 7.5. In each case this numerical factor, K , will be a number read from 
Crank’s curves times a normalization factor obtained by integrating the ^-function. 

Thus, for each f -function used by Crank, K was computed and also the adsorption 
isotherm was calculated and plotted from Eq. Q. 8. To correlate K as a function of 
isotherm shape, the simple non-linear function (NLF) of Eq. 7.11 was calculated for 
each isotherm, and K was plotted as a function of this. The resulting curve was 
then found to be fitted by Eq. 7. 10 to a good degree of accuracy. 

Thus, from Eq. 7.9, 7.10 and 7.11 the time for 50% desorption for any-shaped 
isotherm can be calculated. 
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